[T] only create mesh once

Tests/NumLib/TestExtrapolation.cpp

@@ -12,6 +12,7 @@
 #include "MathLib/LinAlg/LinAlg.h"
 
 #include "MeshLib/MeshGenerators/MeshGenerator.h"
+#include "MeshLib/IO/writeMeshToFile.h"
 
 #include "NumLib/DOF/DOFTableUtil.h"
 #include "NumLib/DOF/MatrixProviderUser.h"
@@ -31,7 +32,6 @@
 
 namespace
 {
-
 template<typename ShapeMatrices>
 void interpolateNodalValuesToIntegrationPoints(
         std::vector<double> const& local_nodal_values,
@@ -45,7 +45,7 @@ void interpolateNodalValuesToIntegrationPoints(
     }
 }
 
-}
+} // anonymous namespace
 
 class LocalAssemblerDataInterface : public NumLib::ExtrapolatableElement
 {
@@ -171,6 +171,7 @@ public:
     {
         auto const extrapolatables =
             NumLib::makeExtrapolatable(_local_assemblers, method);
+
         _extrapolator->extrapolate(extrapolatables);
         _extrapolator->calculateResiduals(extrapolatables);
 
@@ -234,19 +235,18 @@ TEST(NumLib, DISABLED_Extrapolation)
      * x.
      */
 
+    const double mesh_length = 1.0;
+    const std::size_t mesh_elements_in_each_direction = 5;
+
+    // generate mesh
+    std::unique_ptr<MeshLib::Mesh> mesh(
+                MeshLib::MeshGenerator::generateRegularHexMesh(
+                    mesh_length, mesh_elements_in_each_direction));
+
     for (unsigned integration_order : {2, 3, 4})
     {
-
         namespace LinAlg = MathLib::LinAlg;
 
-        const double mesh_length = 1.0;
-        const double mesh_elements_in_each_direction = 5.0;
-
-        // generate mesh
-        std::unique_ptr<MeshLib::Mesh> mesh(
-                    MeshLib::MeshGenerator::generateRegularHexMesh(
-                        mesh_length, mesh_elements_in_each_direction));
-
         auto const nnodes    = mesh->getNumberOfNodes();
         auto const nelements = mesh->getNumberOfElements();
         DBUG("number of nodes: %lu, number of elements: %lu", nnodes, nelements);