Bugfix: adjusted handling of coupled interfaces on decomposition and reconstruction

src/foam/meshes/polyMesh/polyPatches/constraint/ggi/ggiPolyPatch.C

@@ -62,6 +62,15 @@ bool Foam::ggiPolyPatch::active() const
     polyPatchID shadow(shadowName_, boundaryMesh());
     faceZoneID zone(zoneName_, boundaryMesh().mesh().faceZones());
 
+    // For decomposition and reconstruction
+    // If not runing in parallel and the patch is not local, this is a serial
+    // operation on a piece of a parallel decomposition and is therefore
+    // inactive.  HJ, 5/Spe/2016
+    if (!Pstream::parRun() && !localParallel())
+    {
+        return false;
+    }
+
     return shadow.active() && zone.active();
 }
 

src/foam/meshes/polyMesh/polyPatches/constraint/mixingPlane/mixingPlanePolyPatch.C

@@ -58,6 +58,15 @@ bool Foam::mixingPlanePolyPatch::active() const
     polyPatchID shadow(shadowName_, boundaryMesh());
     faceZoneID zone(zoneName_, boundaryMesh().mesh().faceZones());
 
+    // For decomposition and reconstruction
+    // If not runing in parallel and the patch is not local, this is a serial
+    // operation on a piece of a parallel decomposition and is therefore
+    // inactive.  HJ, 5/Sep/2016
+    if (!Pstream::parRun() && !localParallel())
+    {
+        return false;
+    }
+
     return shadow.active() && zone.active();
 }
 

src/foam/meshes/polyMesh/polyPatches/constraint/overlapGgi/overlapGgiPolyPatch.C

@@ -53,6 +53,15 @@ bool Foam::overlapGgiPolyPatch::active() const
     polyPatchID shadow(shadowName_, boundaryMesh());
     faceZoneID zone(zoneName_, boundaryMesh().mesh().faceZones());
 
+    // For decomposition and reconstruction
+    // If not runing in parallel and the patch is not local, this is a serial
+    // operation on a piece of a parallel decomposition and is therefore
+    // inactive.  HJ, 5/Sep/2016
+    if (!Pstream::parRun() && !localParallel())
+    {
+        return false;
+    }
+
     return shadow.active() && zone.active();
 }
 

src/foam/meshes/polyMesh/polyPatches/constraint/regionCouple/regionCouplePolyPatch.C

@@ -70,10 +70,10 @@ bool Foam::regionCouplePolyPatch::active() const
     {
         // Shadow region present
         const polyMesh& sr = boundaryMesh().mesh().db().parent().
-        objectRegistry::lookupObject<polyMesh>
-        (
-            shadowRegionName_
-        );
+            objectRegistry::lookupObject<polyMesh>
+            (
+                shadowRegionName_
+            );
 
         polyPatchID shadowPatch(shadowPatchName_, sr.boundaryMesh());