exaFOAM/foam-extend4.1-coherent-io
Archived
1
1
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity

Bugfix in protectedInitialRefinement strategy

Browse source 
Protect only if the current cell level is equal to the initial cell level,
otherwise allow to unrefine.
This commit is contained in:
Vuko Vukcevic 2018-07-26 07:35:16 +02:00
parent 56144efed0
commit 6974c3f99f
1 changed files with 23 additions and 2 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

25
src/dynamicMesh/topoChangerFvMesh/dynamicPolyRefinementFvMesh/refinementSelection/protectedInitialRefinement/protectedInitialRefinement.C
View file

@ -155,13 +155,19 @@ Foam::protectedInitialRefinement::unrefinementPointCandidates() const
// Loop through all cells and count number of protected points
label nProtectedPoints = 0;
// Get current cell level
const labelIOList& curCellLevel =
mesh().lookupObject<labelIOList>("cellLevel");
forAll (initialCellLevelIn, cellI)
{
const scalar& cl = initialCellLevelIn[cellI];
const scalar curcl = curCellLevel[cellI];
if (cl > SMALL)
if ((cl > SMALL) && equal(cl, curcl))
{
// Cell has been refined during meshing, mark all of its points
// Cell has been refined during meshing and the original level is
// equal to the current level, mark all of its points
const labelList& curCellPoints = meshCellPoints[cellI];
forAll (curCellPoints, i)
@ -178,6 +184,21 @@ Foam::protectedInitialRefinement::unrefinementPointCandidates() const
}
}
}
else if (cl > (curcl + SMALL))
{
FatalErrorIn
(
"Xfer<labelList>\n"
"protectedInitialRefinement::unrefinementPointCandidates()"
" const"
) << "Detected a cell with initial refinement level greater"
<< " than current refinement level." << nl
<< "cellI: " << cellI << nl
<< "current level: " << curcl << nl
<< "initial level: " << cl << nl
<< "This should not happen."
<< abort(FatalError);
}
}
// Create the list for unrefinement candidates