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

Bugfixes in polyhedralRefinement.?

Browse source 
1. TypeName macro instead of ClassName
2. Need all points (original + added) when checking internal/boundary face
   orientation
This commit is contained in:
Vuko Vukcevic 2018-01-24 12:37:57 +01:00
parent 00b7c46868
commit 45d4642a2d
2 changed files with 31 additions and 11 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

40
src/dynamicMesh/dynamicMesh/polyMeshModifiers/polyhedralRefinement/polyhedralRefinement.C
View file

@ -2691,17 +2691,27 @@ void Foam::polyhedralRefinement::checkInternalOrientation
// Get list of polyAddPoint objects
const DynamicList<polyAddPoint>& polyAddedPoints(ref.addedPoints());
// Create a list of newly added points
pointField addedPoints(polyAddedPoints.size());
forAll (addedPoints, i)
// Create a list of all points
const label nOrigPoints = mesh_.nPoints();
pointField allPoints(nOrigPoints + polyAddedPoints.size());
// Set ordinary points first
const pointField& meshPoints = mesh_.points();
forAll (meshPoints, i)
{
addedPoints[i] = polyAddedPoints[i].newPoint();
allPoints[i] = meshPoints[i];
}
// Set newly added points next
forAll (polyAddedPoints, i)
{
allPoints[i + nOrigPoints] = polyAddedPoints[i].newPoint();
}
// Get compact points
const pointField compactPoints
(
IndirectList<point>(addedPoints, newFace)()
IndirectList<point>(allPoints, newFace)()
);
const vector n(compactFace.normal(compactPoints));
@ -2768,17 +2778,27 @@ void Foam::polyhedralRefinement::checkBoundaryOrientation
// Get list of polyAddPoint objects
const DynamicList<polyAddPoint>& polyAddedPoints(ref.addedPoints());
// Create a list of newly added points
pointField addedPoints(polyAddedPoints.size());
forAll (addedPoints, i)
// Create a list of all points
const label nOrigPoints = mesh_.nPoints();
pointField allPoints(nOrigPoints + polyAddedPoints.size());
// Set ordinary points first
const pointField& meshPoints = mesh_.points();
forAll (meshPoints, i)
{
addedPoints[i] = polyAddedPoints[i].newPoint();
allPoints[i] = meshPoints[i];
}
// Set newly added points next
forAll (polyAddedPoints, i)
{
allPoints[i + nOrigPoints] = polyAddedPoints[i].newPoint();
}
// Get compact points
const pointField compactPoints
(
IndirectList<point>(addedPoints, newFace)()
IndirectList<point>(allPoints, newFace)()
);
const vector n(compactFace.normal(compactPoints));

2
src/dynamicMesh/dynamicMesh/polyMeshModifiers/polyhedralRefinement/polyhedralRefinement.H
View file

@ -366,7 +366,7 @@ class polyhedralRefinement
public:
//- Runtime type information
ClassName("polyhedralRefinement");
TypeName("polyhedralRefinement");
// Constructors