Skip to content
Snippets Groups Projects
Commit d9c9d51c authored by Tom Fischer's avatar Tom Fischer
Browse files

[PL/LIE] Subst. MaL::Vector3 by Eigen::Vector3d.

parent a7805df3
No related branches found
No related tags found
No related merge requests found
Hide whitespace changes
Inline Side-by-side
ProcessLib/LIE/Common/Utils.cpp
+ 3
− 3
View file @ d9c9d51c
...@@ -9,7 +9,6 @@ ...@@ -9,7 +9,6 @@
#include "Utils.h" #include "Utils.h"
#include "MathLib/Vector3.h"
#include "MeshLib/Elements/FaceRule.h" #include "MeshLib/Elements/FaceRule.h"
namespace ProcessLib namespace ProcessLib
...@@ -23,8 +22,9 @@ void computeNormalVector(MeshLib::Element const& e, unsigned const global_dim, ...@@ -23,8 +22,9 @@ void computeNormalVector(MeshLib::Element const& e, unsigned const global_dim,
if (global_dim == 2) if (global_dim == 2)
{ {
assert(e.getGeomType() == MeshLib::MeshElemType::LINE); assert(e.getGeomType() == MeshLib::MeshElemType::LINE);
auto v1 = Eigen::Vector3d const v1 =
MathLib::Vector3(*e.getNode(1)) - MathLib::Vector3(*e.getNode(0)); 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[0] = -v1[1];
element_normal[1] = v1[0]; element_normal[1] = v1[0];
element_normal[2] = 0; // not used in 2d but needed for normalization element_normal[2] = 0; // not used in 2d but needed for normalization
......
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment