From 3df5356cfafc7c3ad5f84cd2ce376c460c68391a Mon Sep 17 00:00:00 2001
From: Tom Fischer <thomas.fischer@ufz.de>
Date: Fri, 5 Apr 2024 09:12:00 +0200
Subject: [PATCH] [PL/CT] Add shape matrix cache to process data
---
ProcessLib/ComponentTransport/ComponentTransportProcessData.h | 4 ++++
.../ComponentTransport/CreateComponentTransportProcess.cpp | 2 +-
2 files changed, 5 insertions(+), 1 deletion(-)
diff --git a/ProcessLib/ComponentTransport/ComponentTransportProcessData.h b/ProcessLib/ComponentTransport/ComponentTransportProcessData.h
index 2ac5512eed9..121c3e7641c 100644
--- a/ProcessLib/ComponentTransport/ComponentTransportProcessData.h
+++ b/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;
diff --git a/ProcessLib/ComponentTransport/CreateComponentTransportProcess.cpp b/ProcessLib/ComponentTransport/CreateComponentTransportProcess.cpp
index bb69b4fe94e..5710f1e4d51 100644
--- a/ProcessLib/ComponentTransport/CreateComponentTransportProcess.cpp
+++ b/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;
--
GitLab