zhyang-dev/OpenFOAM-12
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

cuttingPlane: Corrected handling of intersection between a plane and a cellZone in parallel

Browse Source 
Patch contributed by Federico Piscaglia
Resolves bug-report https://bugs.openfoam.org/view.php?id=3283
This commit is contained in:
Henry Weller
2019-05-26 12:01:41 +01:00
parent 0a3b7a7932
commit 47d2e03ae3
1 changed files with 18 additions and 2 deletions
Download Patch File Download Diff File Expand all files Collapse all files

20
src/sampling/cuttingPlane/cuttingPlane.C
View File

@ -2,7 +2,7 @@
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Copyright (C) 2011-2018 OpenFOAM Foundation
\\ / A nd | Copyright (C) 2011-2019 OpenFOAM Foundation
\\/ M anipulation |
-------------------------------------------------------------------------------
License
@ -58,6 +58,12 @@ void Foam::cuttingPlane::calcCutCells
cutCells_.setSize(listSize);
label cutcelli(0);
// fp: check if the list (cellZone) is not empty.
const bool isZoneEmpty
(
(returnReduce(cellIdLabels.size(), sumOp<label>()) > 0) ? false : true
);
// Find the cut cells by detecting any cell that uses points with
// opposing dotProducts.
for (label listI = 0; listI < listSize; ++listI)
@ -68,7 +74,17 @@ void Foam::cuttingPlane::calcCutCells
{
celli = cellIdLabels[listI];
}
else
{
// fp: in parallel computation, if the cellZone exists globally
// but not locally, the postprocessing must be still be limited to
// the crossing plane.
if (!isZoneEmpty)
{
cutCells_.setSize(0);
return;
}
}
const labelList& cEdges = cellEdges[celli];
label nCutEdges = 0;