Load balancing support, polyMesh

src/foam/meshes/polyMesh/polyMesh.C

@@ -719,6 +719,9 @@ void Foam::polyMesh::resetPrimitives
     // Clear addressing. Keep geometric props for mapping.
     clearAddressing();
 
+    // Clear everything
+    clearOut();
+
     // Take over new primitive data.
     // Optimized to avoid overwriting data at all
     if (!pts().empty())
@@ -748,19 +751,16 @@ void Foam::polyMesh::resetPrimitives
 
 
     // Reset patch sizes and starts
+    Pout<< "Resetting patches: starts: " << patchStarts << endl;
     forAll (boundary_, patchI)
     {
-        boundary_[patchI] = polyPatch
+        boundary_[patchI].resetPatch
         (
-            boundary_[patchI].name(),
             patchSizes[patchI],
-            patchStarts[patchI],
-            patchI,
-            boundary_
+            patchStarts[patchI]
         );
     }
 
-
     // Flags the mesh files as being changed
     setInstance(time().timeName());
 
@@ -804,12 +804,10 @@ void Foam::polyMesh::resetPrimitives
         }
     }
 
-
     // Set the primitive mesh from the owner_, neighbour_.
     // Works out from patch end where the active faces stop.
     initMesh();
 
-
     if (validBoundary)
     {
         // Note that we assume that all the patches stay the same and are

src/foam/meshes/polyMesh/polyMesh.H

@@ -539,8 +539,8 @@ public:
             void removeZones();
 
             //- Reset mesh primitive data. Assumes all patch info correct
-            //  (so does e.g. parallel communication). If not use
-            //  validBoundary=false
+            //  (so does e.g. parallel communication). If not, use
+            //  validBoundary = false
             //  (still assumes patchStarts[0] = nInternalFaces and last
             //  patch ends at nActiveFaces) and change patches with addPatches.
             void resetPrimitives

src/foam/meshes/polyMesh/polyPatches/polyPatch/polyPatch.C

@@ -579,6 +579,36 @@ void Foam::polyPatch::syncOrder() const
 {}
 
 
+void Foam::polyPatch::resetPatch
+(
+    const label newSize,
+    const label newStart
+)
+{
+    // Clear all data
+    clearAddressing();
+    
+    // Reset start and primitive patch
+    start_ = newStart;
+    
+    primitivePatch::operator=
+    (
+        primitivePatch
+        (
+            faceSubList
+            (
+                boundaryMesh_.mesh().allFaces(),
+                newSize,
+                newStart
+            ),
+            boundaryMesh_.mesh().allPoints()
+        )
+    );
+
+    Pout<< "Patch reset: faceCells: " << faceCells() << endl;
+}
+
+
 // * * * * * * * * * * * * * * * Member Operators  * * * * * * * * * * * * * //
 
 void Foam::polyPatch::operator=(const polyPatch& p)

src/foam/meshes/polyMesh/polyPatches/polyPatch/polyPatch.H

@@ -464,6 +464,13 @@ public:
         //  but is happening on other processors
         virtual void syncOrder() const;
 
+        //- Reset patch
+        virtual void resetPatch
+        (
+            const label newSize,
+            const label newStart
+        );
+
 
     // Member operators