diff --git a/contrib/DiscreteIntegration/Integration3D.cpp b/contrib/DiscreteIntegration/Integration3D.cpp
index a281fd0f27805c60415ad2196cf0fd12fd71900d..a4ad62b753431f1526a7487a67d2c1cdfbacc3eb 100644
--- a/contrib/DiscreteIntegration/Integration3D.cpp
+++ b/contrib/DiscreteIntegration/Integration3D.cpp
@@ -873,6 +873,7 @@ void DI_Element::evalC (const double u, const double v, const double w, double *
   int nbV = nbVert() + nbMid();
   std::vector<double> s(nbV);
   ev[0] = 0; ev[1] = 0; ev[2] = 0;
+  getShapeFunctions (u, v, w, &s[0], order); 
   for(int i = 0; i < nbV; i++){
       ev[0] += x(i) * s[i];
       ev[1] += y(i) * s[i];