songtianlun/openfoam
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

ENH: polyMeshInitMesh.C: check on validity

Browse Source
This commit is contained in:
mattijs 2011-07-20 16:01:35 +01:00
parent 74e71fe1d7
commit 866c75a9a2
1 changed files with 22 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

22
src/OpenFOAM/meshes/polyMesh/polyMeshInitMesh.C
View File

@ -60,12 +60,26 @@ void Foam::polyMesh::initMesh()
forAll(owner_, facei) forAll(owner_, facei)
{ {
if (owner_[facei] < 0)
{
FatalErrorIn("polyMesh::initMesh()")
<< "Illegal cell label " << owner_[facei]
<< " in neighbour addressing for face " << facei
<< exit(FatalError);
}
nCells = max(nCells, owner_[facei]); nCells = max(nCells, owner_[facei]);
} }
// The neighbour array may or may not be the same length as the owner // The neighbour array may or may not be the same length as the owner
forAll(neighbour_, facei) forAll(neighbour_, facei)
{ {
if (neighbour_[facei] < 0)
{
FatalErrorIn("polyMesh::initMesh()")
<< "Illegal cell label " << neighbour_[facei]
<< " in neighbour addressing for face " << facei
<< exit(FatalError);
}
nCells = max(nCells, neighbour_[facei]); nCells = max(nCells, neighbour_[facei]);
} }
@ -113,6 +127,14 @@ void Foam::polyMesh::initMesh(cellList& c)
forAll(cellfaces, faceI) forAll(cellfaces, faceI)
{ {
if (cellfaces[faceI] < 0)
{
FatalErrorIn("polyMesh::initMesh(cellList&)")
<< "Illegal face label " << cellfaces[faceI]
<< " in cell " << cellI
<< exit(FatalError);
}
if (!markedFaces[cellfaces[faceI]]) if (!markedFaces[cellfaces[faceI]])
{ {
// First visit: owner // First visit: owner