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

Disable layers around corners

Browse source
This commit is contained in:
Hrvoje Jasak 2014-02-10 17:30:25 +00:00
parent 56e392ae55
commit d3c8528c07
1 changed files with 9 additions and 50 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

43
src/dynamicMesh/dynamicMesh/polyMeshModifiers/layerAdditionRemoval/setLayerPairing.C
View file

@ -144,48 +144,8 @@ bool Foam::layerAdditionRemoval::setLayerPairing() const
else
{
// Point mapped from some other face. Check the label
// Layers wrapped around cells are not allowed. HJ, 10/Feb/2014
if (ptc[clp] != lidFace[pointI])
{
// We may be dealing with a layer that is wrapped around
// a cell. Check second point neighbours for common point.
// This is most probably not a general solution for
// arbitray polyhedra!
// HR, 1/May/2011
const label curPoint = faces[mf[faceI]][pointI];
const labelListList& ppAddr = mesh.pointPoints();
const labelList& p1List = ppAddr[lidFace[pointI]];
const labelList& p2List = ppAddr[ptc[clp]];
bool found = false;
forAll (p1List, p1I)
{
const label p1 = p1List[p1I];
if (p1 != curPoint)
{
forAll (p2List, p2I)
{
label p2 = p2List[p2I];
if (p1 == p2)
{
ptc[clp] = p1;
found = true;
break;
}
}
}
if (found)
{
break;
}
}
if (!found)
{
nPointErrors++;
// Pout<< "Topological error in cell layer pairing. "
@ -198,7 +158,6 @@ bool Foam::layerAdditionRemoval::setLayerPairing() const
}
}
}
}
// Pout << "ptc: " << ptc << endl;
}