diff --git a/src/problem/Formulation.cpp b/src/problem/Formulation.cpp
index d9eafcafc16f019184e7abccf521f8ad8e53e37b..c8576b20ee5ca88a508a186f027ffa1c3578653c 100644
--- a/src/problem/Formulation.cpp
+++ b/src/problem/Formulation.cpp
@@ -1230,8 +1230,10 @@ namespace gmshddm
     }
 
     template< class T_Scalar >
-    static int MatVectProduct(Mat A, Vec X, Vec Y)
+    int MatVectProduct(Mat A, Vec X, Vec Y)
     {
+      PetscFunctionBegin;
+
       const int outerVerbosity = gmshfem::common::Options::instance()->verbose;
       const int innerVerbosity = (outerVerbosity != 0 ? outerVerbosity - 1 : 0);