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: medialAxisMover: duplicate list

Browse Source
This commit is contained in:
mattijs
2013-12-17 11:59:31 +00:00
parent 39a8a088c0
commit 313841b989
1 changed files with 7 additions and 9 deletions
Download Patch File Download Diff File Expand all files Collapse all files

16
src/mesh/autoMesh/autoHexMesh/externalDisplacementMeshMover/medialAxisMeshMover.C
View File

@ -1317,16 +1317,14 @@ void Foam::medialAxisMeshMover::findIsolatedRegions
forAll(f, fp) forAll(f, fp)
{ {
label patchPointI = f[fp];
if (extrudeStatus[f[fp]] != autoLayerDriver::NOEXTRUDE) if (extrudeStatus[f[fp]] != autoLayerDriver::NOEXTRUDE)
{ {
if (islandPoint[faceI] == -1) if (islandPoint[faceI] == -1)
{ {
// First point to extrude // First point to extrude
islandPoint[faceI] = patchPointI; islandPoint[faceI] = f[fp];
} }
else else if (islandPoint[faceI] != -2)
{ {
// Second or more point to extrude // Second or more point to extrude
islandPoint[faceI] = -2; islandPoint[faceI] = -2;
@ -1336,12 +1334,11 @@ void Foam::medialAxisMeshMover::findIsolatedRegions
} }
// islandPoint: // islandPoint:
// -1 : no point extruded // -1 : no point extruded on face
// -2 : >= 2 points extruded // -2 : >= 2 points extruded on face
// >=0: label of point extruded // >=0: label of point extruded
// Check all surrounding faces that I am the islandPoint // Check all surrounding faces that I am the islandPoint
boolList keptPoints(pp.nPoints(), false);
forAll(pointFaces, patchPointI) forAll(pointFaces, patchPointI)
{ {
if (extrudeStatus[patchPointI] != autoLayerDriver::NOEXTRUDE) if (extrudeStatus[patchPointI] != autoLayerDriver::NOEXTRUDE)
@ -1422,6 +1419,7 @@ void Foam::medialAxisMeshMover::findIsolatedRegions
} }
} }
if (returnReduce(nChanged, sumOp<label>()) == 0) if (returnReduce(nChanged, sumOp<label>()) == 0)
{ {
break; break;
@ -1773,8 +1771,8 @@ void Foam::medialAxisMeshMover::calculateDisplacement
mesh().globalData().nTotalPoints() mesh().globalData().nTotalPoints()
); );
//- Use strick extrusionIsland detection //- Use strict extrusionIsland detection
const Switch detectExtrusionIsland = coeffDict.lookupOrDefault<label> const Switch detectExtrusionIsland = coeffDict.lookupOrDefault<Switch>
( (
"detectExtrusionIsland", "detectExtrusionIsland",
true true