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

BUG: tetDecomposer: illegal tet face. Fixes #2580

Browse Source
This commit is contained in:
mattijs
2022-09-12 16:34:51 +01:00
parent f9b85dbfbb
commit ff4c776b8b
1 changed files with 25 additions and 2 deletions
Download Patch File Download Diff File Expand all files Collapse all files

27
src/dynamicMesh/polyTopoChange/polyTopoChange/tetDecomposer.C
View File

@ -6,7 +6,7 @@
\\/ M anipulation |
-------------------------------------------------------------------------------
Copyright (C) 2012-2016 OpenFOAM Foundation
Copyright (C) 2015-2020 OpenCFD Ltd.
Copyright (C) 2015-2020,2022 OpenCFD Ltd.
-------------------------------------------------------------------------------
License
This file is part of OpenFOAM.
@ -158,19 +158,42 @@ const
const face& f = mesh_.faces()[facei];
const label fp0 = max(0, mesh_.tetBasePtIs()[facei]);
// Work out triangle index on this face. Assume 'fan' triangulation starting
// from fp0. E.g. 6 vertices on face -> 4 triangles. First and last triangle
// use consecutive vertices
//
// fp | triangle | vertices
// ------|----------|---------
// fp0 | 0 | fp0,fp0+1,fp0+2
// fp0+1 | 0 | ,,
// fp0+2 | 1 | fp0,fp0+2,fp0+3
// fp0+3 | 2 | fp0,fp0+3,fp0+4
// fp0+4 | 3 | fp0,fp0+4,fp0+3
// fp0+5 | 3 | ,,
// Work out triangle index on this face
label thisTrii;
if (fp == fp0)
{
thisTrii = 0;
}
else if (fp == f.fcIndex(fp0))
{
thisTrii = 0;
}
else if (fp == f.rcIndex(fp0))
{
thisTrii = f.size()-3;
}
else if (fp < fp0)
{
const label fpB = fp+f.size();
thisTrii = (fpB-fp0-1);
}
else
{
thisTrii = (fp-fp0-1) % (f.size()-2);
thisTrii = (fp-fp0-1);
}
return thisTrii;
}