diff --git a/qutip/nonmarkov/bofin_solvers.py b/qutip/nonmarkov/bofin_solvers.py
index 86328acd9b..cd199b8424 100644
--- a/qutip/nonmarkov/bofin_solvers.py
+++ b/qutip/nonmarkov/bofin_solvers.py
@@ -634,9 +634,6 @@ def _grad_next_fermionic(self, he_n, k):
             )
         return op
 
-
-        return op
-
     def _rhs(self, L):
         """ Make the RHS for the HEOM. """
         ops = _GatherHEOMRHS(self.ados.idx, block=L.shape[0], nhe=self._n_ados)