diff --git a/NumLib/Fem/CoordinatesMapping/ShapeMatrices-impl.h b/NumLib/Fem/CoordinatesMapping/ShapeMatrices-impl.h
index eeb5ee23787a6e716fe830ace6048e1bd2186bf1..89f010bf63e536dd163bb8b6d88abd98e34c2ac4 100644
--- a/NumLib/Fem/CoordinatesMapping/ShapeMatrices-impl.h
+++ b/NumLib/Fem/CoordinatesMapping/ShapeMatrices-impl.h
@@ -64,6 +64,7 @@ inline void setZero(ShapeMatrices<T_N, T_DNDR, T_J, T_DNDX> &shape, ShapeDataFie
     setZero(shape, ShapeDataFieldType<ShapeMatrixType::DNDR>());
     setMatrixZero(shape.J);
     shape.detJ = .0;
+    shape.integralMeasure = 0.0;
 }
 
 template <class T_N, class T_DNDR, class T_J, class T_DNDX>
diff --git a/NumLib/Fem/CoordinatesMapping/ShapeMatrices.h b/NumLib/Fem/CoordinatesMapping/ShapeMatrices.h
index a1d6d5d8ca46982ce12ea9d1e58a711f0ef3037e..c1c60d3730f0c8a8f0b390c687355a92d05beb2d 100644
--- a/NumLib/Fem/CoordinatesMapping/ShapeMatrices.h
+++ b/NumLib/Fem/CoordinatesMapping/ShapeMatrices.h
@@ -54,6 +54,7 @@ struct ShapeMatrices
     double detJ;        ///< Determinant of the Jacobian
     JacobianType invJ;  ///< Inverse matrix of the Jacobian
     DxShapeType dNdx;   ///< Matrix of gradient of shape functions in physical coordinates, dN(r)/dx
+    double integralMeasure;
 
     /** Not default constructible, dimensions always must be given.
      *
@@ -80,7 +81,8 @@ struct ShapeMatrices
           J(local_dim, local_dim),
           detJ(.0),
           invJ(local_dim, local_dim),
-          dNdx(global_dim, n_nodes)
+          dNdx(global_dim, n_nodes),
+          integralMeasure(0.0)
     {
         setZero();
     }