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

polyMesh::findCell: Call tetBasePtIs for CELL_TETS to ensure parallel consistency

Browse Source 
Resolves additional bug-report http://www.openfoam.org/mantisbt/view.php?id=1544
This commit is contained in:
Henry 2015-06-11 17:45:25 +01:00
parent a5b1708e63
commit d71b839320
1 changed files with 15 additions and 9 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

24
src/OpenFOAM/meshes/polyMesh/polyMesh.C
View File

@ -1483,6 +1483,21 @@ Foam::label Foam::polyMesh::findCell
return -1;
}
if
(
Pstream::parRun()
&& (decompMode == FACE_DIAG_TRIS || decompMode == CELL_TETS)
)
{
// Force construction of face-diagonal decomposition before testing
// for zero cells.
//
// If parallel running a local domain might have zero cells so never
// construct the face-diagonal decomposition which uses parallel
// transfers.
(void)tetBasePtIs();
}
if (decompMode == CELL_TETS)
{
// Advanced search method utilizing an octree
@ -1501,15 +1516,6 @@ Foam::label Foam::polyMesh::findCell
// Approximate search avoiding the construction of an octree
// and cell decomposition
if (Pstream::parRun() && decompMode == FACE_DIAG_TRIS)
{
// Force construction of face-diagonal decomposition before testing
// for zero cells. If parallel running a local domain might have
// zero cells so never construct the face-diagonal decomposition
// (which uses parallel transfers)
(void)tetBasePtIs();
}
// Find the nearest cell centre to this location
label celli = findNearestCell(p);