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

[PL/BC] Use Eigen::Vector3d instead of auto, set face normal.

Set higher dimensions components of the face normal to zero.
parent 43f67600
No related branches found
No related tags found
No related merge requests found
......@@ -106,7 +106,11 @@ private:
// At the moment (2016-09-28) the surface normal is not oriented
// according to the right hand rule
// for correct results it is necessary to multiply the normal with -1
return -MeshLib::FaceRule::getSurfaceNormal(&e).normalized();
Eigen::Vector3d surface_normal =
-MeshLib::FaceRule::getSurfaceNormal(&e).normalized();
auto const zeros_size = 3 - _data.process.getMesh().getDimension();
surface_normal.tail(zeros_size).setZero();
return surface_normal;
}
HCNonAdvectiveFreeComponentFlowBoundaryConditionData const& _data;
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment