From ade940fedded58a6326199e2743ec22d0b6245cd Mon Sep 17 00:00:00 2001
From: Christoph Lehmann <christoph.lehmann@ufz.de>
Date: Wed, 24 Aug 2016 11:41:58 +0200
Subject: [PATCH] [PL] use EigenMapTools in local assemblers.
---
ProcessLib/TES/TESLocalAssembler-impl.h | 16 +++++++---------
1 file changed, 7 insertions(+), 9 deletions(-)
diff --git a/ProcessLib/TES/TESLocalAssembler-impl.h b/ProcessLib/TES/TESLocalAssembler-impl.h
index ba18d2feed5..ddc02a240f3 100644
--- a/ProcessLib/TES/TESLocalAssembler-impl.h
+++ b/ProcessLib/TES/TESLocalAssembler-impl.h
@@ -10,6 +10,7 @@
#define PROCESS_LIB_TES_FEM_IMPL_H_
#include "MaterialLib/Adsorption/Adsorption.h"
+#include "MathLib/LinAlg/Eigen/EigenMapTools.h"
#include "NumLib/Fem/FiniteElement/TemplateIsoparametric.h"
#include "NumLib/Fem/ShapeMatrixPolicy.h"
#include "NumLib/Function/Interpolation.h"
@@ -135,15 +136,12 @@ void TESLocalAssembler<ShapeFunction_, IntegrationMethod_, GlobalDim>::assemble(
// This assertion is valid only if all nodal d.o.f. use the same shape matrices.
assert(local_matrix_size == ShapeFunction::NPOINTS * NODAL_DOF);
- local_M_data.resize(local_matrix_size * local_matrix_size);
- auto local_M = Eigen::Map<NodalMatrixType>(
- local_M_data.data(), local_matrix_size, local_matrix_size);
- local_K_data.resize(local_matrix_size * local_matrix_size);
- auto local_K = Eigen::Map<NodalMatrixType>(
- local_K_data.data(), local_matrix_size, local_matrix_size);
- local_b_data.resize(local_matrix_size);
- auto local_b =
- Eigen::Map<NodalVectorType>(local_b_data.data(), local_matrix_size);
+ auto local_M = MathLib::toZeroedMatrix<NodalMatrixType>(
+ local_M_data, local_matrix_size, local_matrix_size);
+ auto local_K = MathLib::toZeroedMatrix<NodalMatrixType>(
+ local_K_data, local_matrix_size, local_matrix_size);
+ auto local_b = MathLib::toZeroedVector<NodalVectorType>(local_b_data,
+ local_matrix_size);
unsigned const n_integration_points =
_integration_method.getNumberOfPoints();
--
GitLab