From a47bd4ff4bc47d0bc90d240b5fc0c7bb30fce16d Mon Sep 17 00:00:00 2001
From: Dmitri Naumov <github@naumov.de>
Date: Thu, 6 Jan 2022 21:35:03 +0100
Subject: [PATCH] [PL/RCT] Avoid vector sizes mismatch
Using dynamic size vector mapping for the vector in number of integration points.
---
.../RichardsComponentTransportFEM-impl.h | 3 ++-
1 file changed, 2 insertions(+), 1 deletion(-)
diff --git a/ProcessLib/RichardsComponentTransport/RichardsComponentTransportFEM-impl.h b/ProcessLib/RichardsComponentTransport/RichardsComponentTransportFEM-impl.h
index cb2aa7d26f5..c2201fc0402 100644
--- a/ProcessLib/RichardsComponentTransport/RichardsComponentTransportFEM-impl.h
+++ b/ProcessLib/RichardsComponentTransport/RichardsComponentTransportFEM-impl.h
@@ -350,7 +350,8 @@ LocalAssemblerData<ShapeFunction, IntegrationMethod, GlobalDim>::
auto const local_x = x[process_id]->get(indices);
cache.clear();
- auto cache_vec = MathLib::createZeroedVector<LocalVectorType>(
+ auto cache_vec = MathLib::createZeroedVector<
+ Eigen::Matrix<double, 1, Eigen::Dynamic, Eigen::RowMajor>>(
cache, n_integration_points);
for (unsigned ip = 0; ip < n_integration_points; ++ip)
--
GitLab