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: orientedSurface: multiple parts of surfaces not visited

Browse Source
This commit is contained in:
mattijs
2011-06-22 14:31:21 +01:00
parent e5f921e765
commit cf51528d49
1 changed files with 27 additions and 5 deletions
Download Patch File Download Diff File Expand all files Collapse all files

30
src/meshTools/triSurface/orientedSurface/orientedSurface.C
View File

@ -102,7 +102,7 @@ Foam::labelList Foam::orientedSurface::edgeToFace
{ {
if (flip[face1] == UNVISITED) if (flip[face1] == UNVISITED)
{ {
FatalErrorIn("orientedSurface::edgeToFace") << "Problem" FatalErrorIn("orientedSurface::edgeToFace(..)") << "Problem"
<< abort(FatalError); << abort(FatalError);
} }
else else
@ -283,8 +283,10 @@ bool Foam::orientedSurface::flipSurface
} }
} }
// Recalculate normals // Recalculate normals
if (hasFlipped)
{
s.clearOut(); s.clearOut();
}
return hasFlipped; return hasFlipped;
} }
@ -352,8 +354,28 @@ bool Foam::orientedSurface::orient
// FLIP: need to flip // FLIP: need to flip
labelList flipState(s.size(), UNVISITED); labelList flipState(s.size(), UNVISITED);
flipState[0] = NOFLIP; label faceI = 0;
walkSurface(s, 0, flipState); while (true)
{
label startFaceI = -1;
while (faceI < s.size())
{
if (flipState[faceI] == UNVISITED)
{
startFaceI = faceI;
break;
}
faceI++;
}
if (startFaceI == -1)
{
break;
}
flipState[startFaceI] = NOFLIP;
walkSurface(s, startFaceI, flipState);
}
anyFlipped = flipSurface(s, flipState); anyFlipped = flipSurface(s, flipState);
} }