diff --git a/ProcessLib/NeumannBcAssembler.h b/ProcessLib/NeumannBcAssembler.h
index 48ad268a67b7fc7883fcd2b5389708a91ce5d880..f731702f24fb2f4226d69c334fade0da4f51f9dd 100644
--- a/ProcessLib/NeumannBcAssembler.h
+++ b/ProcessLib/NeumannBcAssembler.h
@@ -66,7 +66,7 @@ public:
                                          ShapeFunction::NPOINTS);
             fe.computeShapeFunctions(
                     integration_method.getWeightedPoint(ip).getCoords(),
-                    _shape_matrices[ip]);
+                    _shape_matrices[ip], GlobalDim);
         }
 
         _neumann_bc_value = value_lookup(e);
diff --git a/ProcessLib/Utils/InitShapeMatrices.h b/ProcessLib/Utils/InitShapeMatrices.h
index 874effec0107580b0e90c3d34f3af5374bc0a3d7..08e1f06d5226af793ea2ab623a3de16f7c65d907 100644
--- a/ProcessLib/Utils/InitShapeMatrices.h
+++ b/ProcessLib/Utils/InitShapeMatrices.h
@@ -40,7 +40,7 @@ initShapeMatrices(MeshLib::Element const& e, unsigned integration_order)
                                      ShapeFunction::NPOINTS);
         fe.computeShapeFunctions(
                 integration_method.getWeightedPoint(ip).getCoords(),
-                shape_matrices[ip]);
+                shape_matrices[ip], GlobalDim);
     }
 
     return shape_matrices;