diff --git a/openmc/deplete/abc.py b/openmc/deplete/abc.py
index e8f65f887..691606681 100644
--- a/openmc/deplete/abc.py
+++ b/openmc/deplete/abc.py
@@ -133,27 +133,5 @@ class Operator(metaclass=ABCMeta):
 
         pass
 
-    @abstractmethod
-    def form_matrix(self, y, mat):
-        """Forms the f(y) matrix in y' = f(y)y.
-
-        Nominally a depletion matrix, this is abstracted on the off chance
-        that the function f has nothing to do with depletion at all.
-
-        Parameters
-        ----------
-        y : numpy.ndarray
-            An array representing y.
-        mat : int
-            Material id.
-
-        Returns
-        -------
-        scipy.sparse.csr_matrix
-            Sparse matrix representing f(y).
-        """
-
-        pass
-
     def finalize(self):
         pass
diff --git a/openmc/deplete/openmc_wrapper.py b/openmc/deplete/openmc_wrapper.py
index 7d776825b..56a248c30 100644
--- a/openmc/deplete/openmc_wrapper.py
+++ b/openmc/deplete/openmc_wrapper.py
@@ -278,24 +278,6 @@ class OpenMCOperator(Operator):
             number = density * 1.0e24
             self.number.set_atom_density(mat_id, nuclide, number)
 
-    def form_matrix(self, y, mat):
-        """Forms the depletion matrix.
-
-        Parameters
-        ----------
-        y : numpy.ndarray
-            An array representing reaction rates for this material.
-        mat : int
-            Material id.
-
-        Returns
-        -------
-        scipy.sparse.csr_matrix
-            Sparse matrix representing the depletion matrix.
-        """
-
-        return copy.deepcopy(self.chain.form_matrix(y[mat, :, :]))
-
     def initial_condition(self):
         """Performs final setup and returns initial condition.