diff --git a/Numeric/nodalBasis.h b/Numeric/nodalBasis.h
index 3570b76959484e02cd0ab034253d83ebe8eea954..9a614d1021e0b074e047b6c65a14871c43486e96 100644
--- a/Numeric/nodalBasis.h
+++ b/Numeric/nodalBasis.h
@@ -20,7 +20,7 @@ class nodalBasis {
 
   // Basis functions evaluation
   inline virtual void f(double u, double v, double w, double *sf) const {};
-  inline void f(fullMatrix<double> &coord, fullMatrix<double> &sf) const;
+  inline void f(fullMatrix<double> &coord, fullMatrix<double> &sf) const {};
 
   // Basis functions gradients evaluation
   inline virtual void df(double u, double v, double w, double grads[][3]) const {};
diff --git a/wrappers/gmshpy/gmshNumeric.i b/wrappers/gmshpy/gmshNumeric.i
index 4cb1bb74cee0787da44b36af64fb5802c523a24c..cc7712bac2f0b914c6b4df8c6e0b447b849e967e 100644
--- a/wrappers/gmshpy/gmshNumeric.i
+++ b/wrappers/gmshpy/gmshNumeric.i
@@ -10,6 +10,7 @@
   #include "GaussIntegration.h"
   #include "JacobianBasis.h"
   #include "fullMatrix.h"
+  #include "nodalBasis.h"
   #include "polynomialBasis.h"
 %}
 
@@ -18,4 +19,5 @@
 %include "fullMatrix.h"
 %template(fullMatrixDouble) fullMatrix<double>;
 %template(fullVectorDouble) fullVector<double>;
+%include "nodalBasis.h"
 %include "polynomialBasis.h"