diff --git a/Tests/NumLib/TestFunctionInterpolation.cpp b/Tests/NumLib/TestFunctionInterpolation.cpp
index 1d2335e32a4802d766ec3d4eaca6599130a05a02..dd43128f41867a126e6a72b594343461cddad83c 100644
--- a/Tests/NumLib/TestFunctionInterpolation.cpp
+++ b/Tests/NumLib/TestFunctionInterpolation.cpp
@@ -67,7 +67,8 @@ TEST(NumLibFunctionInterpolationTest, Linear1DElement)
     const unsigned integration_order = 2;
     IntegrationMethod integration_method(integration_order);
 
-    ShapeMatricesType::ShapeMatrices shape_matrix;
+    ShapeMatricesType::ShapeMatrices shape_matrix(
+        ShapeFunction::DIM, ShapeFunction::DIM, ShapeFunction::NPOINTS);
 
     finite_element.computeShapeFunctions(
             integration_method.getWeightedPoint(0).getCoords(),