Skip to content

fix/tet-geom

Michael Turner requested to merge meshing/nektar:fix/tetgeom into master

There is currently a tolerance used in the tetgeom.cpp code which is set to 1e-12.

When the tetrahedron is large (edge lengths of 50+) this small tolerance causes the orientation to fail.

Elements will report that they are invalid when they are not and if in debug any nektar code will assert out.

This MR puts in a small patch which normalises the vector so the test is always done on a order 1 problem and isnt scale dependant.

This has no impact on any of the existing tests.

Merge request reports