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

Generalisation of polyRef::getAnchorLevel

Browse source 
The function now returns the correct anchor level for arbitrary polygonal faces,
i.e. there is no assumption that the face comes from possibly spit hex cell.
This commit is contained in:
Vuko Vukcevic 2017-09-08 15:07:53 +02:00
parent 7372970256
commit 3ebae28896
2 changed files with 5 additions and 36 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

37
src/dynamicMesh/dynamicMesh/directTopoChange/directTopoChange/directActions/polyRef.C
View file

@ -799,39 +799,6 @@ Foam::label Foam::polyRef::getAnchorLevel(const label faceI) const
}
}
}
if (debug)
{
// Print out number of points per level
Info<< "nPointsPerLevel: " << nPointsPerLevel << endl;
}
// Find level with maximum number of points
const label maxLevelI = findMax(nPointsPerLevel);
// Check whether there are at least three points with the level <= of this
// level with maximum number of occurences
label nPointsBelow = 0;
for (label belowLevelI = 0; belowLevelI <= maxLevelI; ++belowLevelI)
{
nPointsBelow += nPointsPerLevel[belowLevelI];
}
if (nPointsBelow < 3)
{
// For some reason, the minimum requirement that at least three points
// have level below the anchor level is not met
WarningIn("label polyRef::getAnchorLevel(const label faceI) const")
<< "Found less than three points below an anchor level."
<< nl
<< "Face: " << faceI
<< ", anchor level: " << allLevels[maxLevelI]
<< ", n points below level: " << nPointsBelow
<< ", n total points of this face: " << nPoints
<< endl;
}
return allLevels[maxLevelI];
}
@ -3842,8 +3809,8 @@ Foam::labelListList Foam::polyRef::setRefinement
// Unsplit face that hasn't been handled. Add edge splits to
// face.
const face& f = meshFaces[faceI];
const labelList& fEdges = meshFaceEdges[faceI];
const face& f = mesh_.faces()[faceI];
const labelList& fEdges = mesh_.faceEdges()[faceI];
dynamicLabelList newFaceVerts(f.size());

4
src/dynamicMesh/dynamicMesh/directTopoChange/directTopoChange/directActions/polyRef.H
View file

@ -385,7 +385,9 @@ public:
}
// Refinement
//- Gets level such that the face has four points <= level.
//- Gets level such that the face has n points <= level, where n is
// the number of points (or edges) for a face.
label getAnchorLevel(const label faceI) const;
//- Helper: get points of a cell without using cellPoints addressing