diff --git a/src/algebra/EigenvalueSolver.hpp b/src/algebra/EigenvalueSolver.hpp
index 98696152215164477c34e4743accb5ad4659e85d..de7faf8e42d03ee59a5dfe3b7f783ec99b1695b8 100644
--- a/src/algebra/EigenvalueSolver.hpp
+++ b/src/algebra/EigenvalueSolver.hpp
@@ -128,17 +128,6 @@ class EigenvalueSolver
     }
   }
 
-  template <typename MatrixType>
-  void
-  computeExpForSymmetricMatrix([[maybe_unused]] const MatrixType& A, [[maybe_unused]] SmallMatrix<double>& expA)
-  {
-#ifdef PUGS_HAS_SLEPC
-    this->computeExpForSymmetricMatrix(PETScAijMatrixEmbedder{A}, expA);
-#else    // PUGS_HAS_SLEPC
-    throw NotImplementedError("SLEPc is required to solve eigenvalue problems");
-#endif   // PUGS_HAS_SLEPC
-  }
-
   template <typename T>
   PUGS_INLINE TinyMatrix<3, 3, T> swap(TinyMatrix<3, 3, T>& matrix, size_t i, size_t j) const;