Removal of now-obsolete components

Author: EdCaunt

devitoboundary/stencils.py

@@ -406,55 +406,3 @@ def _all_variants(self, deriv):
                 self._stencil_list[le][ri] = sp.simplify(stencil_entry)
 
         self._stencil_dict[str(self._bcs)+str(self._s_o)+str(deriv)+'ns'] = self._stencil_list
-
-    def subs(self, eta_l=None, eta_r=None):
-        """
-        Obtain a numpy array of the stencil coefficients given values of eta_l
-        and eta_r. This is the offset between the central stencil point and
-        the boundary on the respective side. As such eta_l should always be
-        negative and eta_r positive.
-
-        Parameters
-        ----------
-        eta_l : float
-            The offset between the center of the stencil and the left side
-            boundary, measured in grid increments. Should always be negative.
-            Default is None.
-        eta_r : float
-            The offset between the center of the stencil and the right side
-            boundary, measured in grid increments. Should always be positive.
-            Default is None
-
-        Returns
-        -------
-        coeffs : ndarray
-            An array of the stencil coefficients for these eta values.
-        """
-
-        if eta_l is None or abs(eta_l) > self._s_o/2:
-            dist_l = 0
-            sub_l = 0  # Placeholder to stop subs() encountering None
-        else:
-            if eta_l >= 0:
-                raise ValueError("eta_l must be negative. Current value is %.6f" % eta_l)
-            # Turn eta into the relevent le index in the stencil list
-            dist_l = self._s_o - np.ceil(abs(eta_l)*2).astype(np.int) + 1
-            sub_l = eta_l
-
-        if eta_r is None or eta_r > self._s_o/2:
-            dist_r = 0
-            sub_r = 0  # Placeholder to stop subs() encountering None
-        else:
-            if eta_r <= 0:
-                raise ValueError("eta_r must be positive. Current value is %.6f" % eta_r)
-            # Turn eta into the relevent ri index in the stencil list
-            dist_r = self._s_o - np.ceil(eta_r*2).astype(np.int) + 1
-            sub_r = eta_r
-
-        stencil = self._stencil_list[dist_l][dist_r].subs([(self._eta_l, sub_l),
-                                                           (self._eta_r, sub_r)])
-        coeffs = np.empty(self._s_o+1)
-        for i in range(self._s_o+1):
-            coeffs[i] = float(stencil.coeff(self._f[i-int(self._s_o/2)], 1))
-
-        return coeffs