[PL/CT] Add shape matrix cache to process data

ProcessLib/ComponentTransport/ComponentTransportProcessData.h

@@ -17,6 +17,7 @@
 #include "LookupTable.h"
 #include "MaterialLib/MPL/MaterialSpatialDistributionMap.h"
 #include "MathLib/LinAlg/Eigen/EigenMapTools.h"
+#include "NumLib/Fem/ShapeMatrixCache.h"
 #include "NumLib/NumericalStability/NumericalStabilization.h"
 #include "ParameterLib/ConstantParameter.h"
 #include "ParameterLib/Parameter.h"
@@ -84,6 +85,9 @@ struct ComponentTransportProcessData
 
     bool const isothermal;
 
+    /// caches for each mesh element type the shape matrix
+    NumLib::ShapeMatrixCache shape_matrix_cache;
+
     static const int hydraulic_process_id = 0;
     // Thermal process is optional, indicated by -1. If present, it is positive.
     const int thermal_process_id = isothermal ? -1 : 1;

ProcessLib/ComponentTransport/CreateComponentTransportProcess.cpp

@@ -317,7 +317,7 @@ std::unique_ptr<Process> createComponentTransportProcess(
         mesh_space_dimension,
         *aperture_size_parameter,
         isothermal,
-    };
+        NumLib::ShapeMatrixCache{integration_order}};
 
     SecondaryVariableCollection secondary_variables;