diff --git a/ProcessLib/BoundaryCondition/NormalTractionBoundaryConditionLocalAssembler.h b/ProcessLib/BoundaryCondition/NormalTractionBoundaryConditionLocalAssembler.h
index facd6e966b9b7dbb0d4260884cd9e3e9945e5718..7453f758bb9f87406322705765cd1fb59ec697f5 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();