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

Bug fix: topological changes

Browse source
This commit is contained in:
Hrvoje Jasak 2011-03-11 16:32:42 +00:00
parent bb20fc0848
commit fd95610120
1 changed files with 15 additions and 6 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

21
src/OpenFOAM/meshes/polyMesh/zones/faceZone/faceZone.C
View file

@ -192,26 +192,35 @@ void Foam::faceZone::calcCellLayers() const
if (!faceFlip[faceI])
{
// Face is oriented correctly, no flip needed
// Master is the neighbour, with normal pointing into it
// Bug fix, HJ, 7/Mar/2011
if (mesh.isInternalFace(mf[faceI]))
{
curMc = nei[mf[faceI]];
curSc = own[mf[faceI]];
}
if (mf[faceI] < mesh.nFaces())
else if (mf[faceI] < mesh.nFaces())
{
curMc = -1;
curSc = own[mf[faceI]];
}
}
else
{
if (mf[faceI] < mesh.nFaces())
{
curMc = own[mf[faceI]];
}
// Face flip
// Master is the owner, with normal pointing into it
// Bug fix, HJ, 7/Mar/2011
if (mesh.isInternalFace(mf[faceI]))
{
curMc = own[mf[faceI]];
curSc = nei[mf[faceI]];
}
else if (mf[faceI] < mesh.nFaces())
{
curMc = own[mf[faceI]];
curSc = -1;
}
}
mc[faceI] = curMc;