From c095eb4c36715b47a7c433be7705b8cbb52c71f6 Mon Sep 17 00:00:00 2001
From: graham <g.macpherson@opencfd.co.uk>
Date: Fri, 24 Jun 2011 14:17:04 +0100
Subject: [PATCH] ENH: Use reduce to ensure parallel filtering is synchronised.

---
 .../conformalVoronoiMeshCalcDualMesh.C             | 14 +++++++++++++-
 1 file changed, 13 insertions(+), 1 deletion(-)

diff --git a/src/mesh/conformalVoronoiMesh/conformalVoronoiMesh/conformalVoronoiMeshCalcDualMesh.C b/src/mesh/conformalVoronoiMesh/conformalVoronoiMesh/conformalVoronoiMeshCalcDualMesh.C
index 9a54c1361a..fb10902af7 100644
--- a/src/mesh/conformalVoronoiMesh/conformalVoronoiMesh/conformalVoronoiMeshCalcDualMesh.C
+++ b/src/mesh/conformalVoronoiMesh/conformalVoronoiMesh/conformalVoronoiMeshCalcDualMesh.C
@@ -82,6 +82,8 @@ void Foam::conformalVoronoiMesh::calcDualMesh
     {
         label nInitialBadQualityFaces = checkPolyMeshQuality(points).size();
 
+        reduce(nInitialBadQualityFaces, sumOp<label>());
+
         Info<< nl << "Initial check before face collapse, found "
             << nInitialBadQualityFaces << " bad quality faces"
             << endl;
@@ -145,10 +147,16 @@ void Foam::conformalVoronoiMesh::calcDualMesh
 
                 nBadQualityFaces = wrongFaces.size();
 
+                reduce(nBadQualityFaces, sumOp<label>());
+
                 Info<< nl << "Found " << nBadQualityFaces
                     << " bad quality faces" << endl;
 
-                if (lastWrongFaces == wrongFaces)
+                bool sameFacesAsLastTime(lastWrongFaces == wrongFaces);
+
+                reduce(sameFacesAsLastTime, andOp<bool>());
+
+                if (sameFacesAsLastTime)
                 {
                     Info<< nl << "Consecutive iterations found the same set "
                         << "of bad quality faces." << endl;
@@ -523,6 +531,8 @@ void Foam::conformalVoronoiMesh::smoothSurface
             dualPtIndexMap
         );
 
+        reduce(nCollapsedFaces, sumOp<label>());
+
         reindexDualVertices(dualPtIndexMap);
 
         mergeCloseDualVertices(pts, boundaryPts);
@@ -759,6 +769,8 @@ void Foam::conformalVoronoiMesh::collapseFaces
             deferredCollapseFaces
         );
 
+        reduce(nCollapsedFaces, sumOp<label>());
+
         reindexDualVertices(dualPtIndexMap);
 
         mergeCloseDualVertices(pts, boundaryPts);