diff --git a/src/post/SaveFunction.cpp b/src/post/SaveFunction.cpp
index 1a15f40f8d7469cb06a9c2842ceb83351ccc1ef1..ddacf019d25b153f74d94d47be82e71606849059 100644
--- a/src/post/SaveFunction.cpp
+++ b/src/post/SaveFunction.cpp
@@ -27,7 +27,7 @@ namespace gmshfem::post
   }
 
   template< class T_Scalar >
-  static void s_addModelData(const int tag, const int step, const double time, const int partition, const unsigned int nbrNodesByElements, std::vector< std::size_t > &gmshElementsTags, std::vector< std::vector< double > > &gmshData, std::vector< typename MathObject< T_Scalar, Degree::Degree0 >::Object, numa::allocator< typename MathObject< T_Scalar, Degree::Degree0 >::Object > > &values)
+  void s_addModelData(const int tag, const int step, const double time, const int partition, const unsigned int nbrNodesByElements, std::vector< std::size_t > &gmshElementsTags, std::vector< std::vector< double > > &gmshData, std::vector< typename MathObject< T_Scalar, Degree::Degree0 >::Object, numa::allocator< typename MathObject< T_Scalar, Degree::Degree0 >::Object > > &values)
   {
     if(scalar::IsComplex< T_Scalar >::value) {
       // real part
@@ -78,7 +78,7 @@ namespace gmshfem::post
   }
 
   template< class T_Scalar >
-  static void s_addModelData(const int tag, const int step, const double time, const int partition, const unsigned int nbrNodesByElements, std::vector< std::size_t > &gmshElementsTags, std::vector< std::vector< double > > &gmshData, std::vector< typename MathObject< T_Scalar, Degree::Degree1 >::Object, numa::allocator< typename MathObject< T_Scalar, Degree::Degree1 >::Object > > &values)
+  void s_addModelData(const int tag, const int step, const double time, const int partition, const unsigned int nbrNodesByElements, std::vector< std::size_t > &gmshElementsTags, std::vector< std::vector< double > > &gmshData, std::vector< typename MathObject< T_Scalar, Degree::Degree1 >::Object, numa::allocator< typename MathObject< T_Scalar, Degree::Degree1 >::Object > > &values)
   {
     if(scalar::IsComplex< T_Scalar >::value) {
       // real part
@@ -138,7 +138,7 @@ namespace gmshfem::post
   }
 
   template< class T_Scalar >
-  static void s_addModelData(const int tag, const int step, const double time, const int partition, const unsigned int nbrNodesByElements, std::vector< std::size_t > &gmshElementsTags, std::vector< std::vector< double > > &gmshData, std::vector< typename MathObject< T_Scalar, Degree::Degree2 >::Object, numa::allocator< typename MathObject< T_Scalar, Degree::Degree2 >::Object > > &values)
+  void s_addModelData(const int tag, const int step, const double time, const int partition, const unsigned int nbrNodesByElements, std::vector< std::size_t > &gmshElementsTags, std::vector< std::vector< double > > &gmshData, std::vector< typename MathObject< T_Scalar, Degree::Degree2 >::Object, numa::allocator< typename MathObject< T_Scalar, Degree::Degree2 >::Object > > &values)
   {
     if(scalar::IsComplex< T_Scalar >::value) {
       // real part
diff --git a/src/post/SaveFunction.h b/src/post/SaveFunction.h
index 7c3932f4ad58fabd90092365ca5a461b9ad713e4..7ba626c49bbcc507883242385c6ee7cd22f6f740 100644
--- a/src/post/SaveFunction.h
+++ b/src/post/SaveFunction.h
@@ -17,6 +17,15 @@
 namespace gmshfem::post
 {
 
+  template< class T_Scalar>
+  void s_addModelData(const int tag, const int step, const double time, const int partition, const unsigned int nbrNodesByElements, std::vector< std::size_t > &gmshElementsTags, std::vector< std::vector< double > > &gmshData, std::vector< typename MathObject< T_Scalar, Degree::Degree0 >::Object, numa::allocator< typename MathObject< T_Scalar, Degree::Degree0 >::Object > > &values);
+
+  template< class T_Scalar>
+  void s_addModelData(const int tag, const int step, const double time, const int partition, const unsigned int nbrNodesByElements, std::vector< std::size_t > &gmshElementsTags, std::vector< std::vector< double > > &gmshData, std::vector< typename MathObject< T_Scalar, Degree::Degree1 >::Object, numa::allocator< typename MathObject< T_Scalar, Degree::Degree1 >::Object > > &values);
+
+  template< class T_Scalar>
+  void s_addModelData(const int tag, const int step, const double time, const int partition, const unsigned int nbrNodesByElements, std::vector< std::size_t > &gmshElementsTags, std::vector< std::vector< double > > &gmshData, std::vector< typename MathObject< T_Scalar, Degree::Degree2 >::Object, numa::allocator< typename MathObject< T_Scalar, Degree::Degree2 >::Object > > &values);
+
 
   template< class T_Scalar, Degree T_Degree >
   class SaveFunction final : public PostInterface