Commit d9c9d51c authored by Tom Fischer's avatar Tom Fischer

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

parent a7805df3
......@@ -9,7 +9,6 @@
#include "Utils.h"
#include "MathLib/Vector3.h"
#include "MeshLib/Elements/FaceRule.h"
namespace ProcessLib
......@@ -23,8 +22,9 @@ void computeNormalVector(MeshLib::Element const& e, unsigned const global_dim,
if (global_dim == 2)
{
assert(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[2] = 0; // not used in 2d but needed for normalization
......
Markdown is supported
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