From a7805df3d485263c1330b89c42fa0f33376054ce Mon Sep 17 00:00:00 2001
From: Thomas Fischer <thomas.fischer@ufz.de>
Date: Sat, 26 Dec 2020 14:31:38 +0100
Subject: [PATCH] [PL/BC/NormalTraction] Subst. MaL::Vector3 by
Eigen::Vector3d.
---
.../NormalTractionBoundaryConditionLocalAssembler.h | 6 +++---
1 file changed, 3 insertions(+), 3 deletions(-)
diff --git a/ProcessLib/BoundaryCondition/NormalTractionBoundaryConditionLocalAssembler.h b/ProcessLib/BoundaryCondition/NormalTractionBoundaryConditionLocalAssembler.h
index facd6e966b9..7453f758bb9 100644
--- a/ProcessLib/BoundaryCondition/NormalTractionBoundaryConditionLocalAssembler.h
+++ b/ProcessLib/BoundaryCondition/NormalTractionBoundaryConditionLocalAssembler.h
@@ -11,7 +11,6 @@
#pragma once
#include "MathLib/LinAlg/Eigen/EigenMapTools.h"
-#include "MathLib/Vector3.h"
#include "MeshLib/Elements/FaceRule.h"
#include "NumLib/DOF/DOFTableUtil.h"
#include "ParameterLib/Parameter.h"
@@ -87,8 +86,9 @@ public:
// TODO Extend to rotated 2d meshes and line elements.
if (e.getGeomType() == MeshLib::MeshElemType::LINE)
{
- auto v1 = MathLib::Vector3(*e.getNode(1)) -
- MathLib::Vector3(*e.getNode(0));
+ Eigen::Vector3d const v1 =
+ Eigen::Map<Eigen::Vector3d const>(e.getNode(1)->getCoords()) -
+ Eigen::Map<Eigen::Vector3d const>(e.getNode(0)->getCoords());
element_normal[0] = -v1[1];
element_normal[1] = v1[0];
element_normal.normalize();
--
GitLab