Skip to content
Snippets Groups Projects
Commit 5f3d5870 authored by Karsten Rink's avatar Karsten Rink
Browse files

removed convexity check for zero-volume quads

parent 61c72181
No related branches found
No related tags found
No related merge requests found
...@@ -125,8 +125,10 @@ ElementErrorCode TemplateQuad<NNODES,CELLQUADTYPE>::validate() const ...@@ -125,8 +125,10 @@ ElementErrorCode TemplateQuad<NNODES,CELLQUADTYPE>::validate() const
ElementErrorCode error_code; ElementErrorCode error_code;
error_code[ElementErrorFlag::ZeroVolume] = this->hasZeroVolume(); error_code[ElementErrorFlag::ZeroVolume] = this->hasZeroVolume();
error_code[ElementErrorFlag::NonCoplanar] = (!GeoLib::pointsOnAPlane(*_nodes[0], *_nodes[1], *_nodes[2], *_nodes[3])); error_code[ElementErrorFlag::NonCoplanar] = (!GeoLib::pointsOnAPlane(*_nodes[0], *_nodes[1], *_nodes[2], *_nodes[3]));
error_code[ElementErrorFlag::NonConvex] = (!(GeoLib::dividedByPlane(*_nodes[0], *_nodes[2], *_nodes[1], *_nodes[3]) && // for collapsed quads (i.e. reduced to a line) this test might result "false" as all four points are actually located on a line.
GeoLib::dividedByPlane(*_nodes[1], *_nodes[3], *_nodes[0], *_nodes[2]))); if (!error_code[ElementErrorFlag::ZeroVolume])
error_code[ElementErrorFlag::NonConvex] = (!(GeoLib::dividedByPlane(*_nodes[0], *_nodes[2], *_nodes[1], *_nodes[3]) &&
GeoLib::dividedByPlane(*_nodes[1], *_nodes[3], *_nodes[0], *_nodes[2])));
return error_code; return error_code;
} }
......
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