zhyang-dev/openfoam
1
0
Fork 0
mirror of https://develop.openfoam.com/Development/openfoam.git synced 2025-11-28 03:28:01 +00:00
Code Issues Actions Packages Projects Releases Wiki Activity

ENH: reduce some isoSurfaceCell regularise failures (closes #870)

Browse Source 
- when the iso-surface (eg, a cut plane, or distanceSurface) passes
  exactly through the cell centre, zero-sized triangles can occur.
  These should be merged away.
This commit is contained in:
Mark Olesen
2018-06-14 14:31:07 +02:00
parent e7d81dfe04
commit 337aca5fd5
1 changed files with 11 additions and 9 deletions
Download Patch File Download Diff File Expand all files Collapse all files

20
src/sampling/surface/isoSurface/isoSurfaceCell.C
View File

@ -284,16 +284,18 @@ Foam::pointIndexHit Foam::isoSurfaceCell::collapseSurface
if (shared[0] != -1)
{
vector n0 = tri0.normal(localPoints);
n0 /= mag(n0);
vector n1 = tri1.normal(localPoints);
n1 /= mag(n1);
const vector n0 = tri0.normal(localPoints);
const vector n1 = tri1.normal(localPoints);
if ((n0 & n1) < 0)
{
// Too big an angle. Do not simplify.
}
else
// Merge any zero-sized triangles,
// or if they point in the same direction.
if
(
mag(n0) <= ROOTVSMALL
|| mag(n1) <= ROOTVSMALL
|| (n0 & n1) >= 0
)
{
info.setPoint
(