domainDecomposition: Correct stitcher construction on read-update

src/parallel/parallel/domainDecomposition.C

@@ -109,6 +109,8 @@ void Foam::domainDecomposition::readComplete()
             false
         )
     );
+
+    completeMesh_->postConstruct(false, fvMesh::stitchType::none);
 }
 
 
@@ -132,6 +134,8 @@ void Foam::domainDecomposition::readProcs()
                 false
             )
         );
+
+        procMeshes_[proci].postConstruct(false, fvMesh::stitchType::none);
     }
 }
 
@@ -516,12 +520,7 @@ bool Foam::domainDecomposition::readDecompose()
 
 void Foam::domainDecomposition::postReadDecompose()
 {
-    completeMesh_->postConstruct(false, fvMesh::stitchType::nonGeometric);
+    completeMesh_->stitcher().connect(false, false, true);
-
-    forAll(procMeshes_, proci)
-    {
-        procMeshes_[proci].postConstruct(false, fvMesh::stitchType::none);
-    }
 }
 
 
@@ -663,15 +662,9 @@ bool Foam::domainDecomposition::readReconstruct()
 
 void Foam::domainDecomposition::postReadReconstruct()
 {
-    completeMesh_->postConstruct(false, fvMesh::stitchType::none);
-
     forAll(procMeshes_, proci)
     {
-        procMeshes_[proci].postConstruct
+        procMeshes_[proci].stitcher().connect(false, false, true);
-        (
-            false,
-            fvMesh::stitchType::nonGeometric
-        );
     }
 }
 
@@ -776,14 +769,6 @@ void Foam::domainDecomposition::postReadUpdateDecompose
     const fvMesh::readUpdateState stat
 )
 {
-    if (stat >= fvMesh::TOPO_CHANGE)
-    {
-        forAll(procMeshes_, proci)
-        {
-            procMeshes_[proci].postConstruct(false, fvMesh::stitchType::none);
-        }
-    }
-
     if (completeMesh_->stitcher().stitches() && stat != fvMesh::UNCHANGED)
     {
         procFaceAddressingBf_.clear();
@@ -883,11 +868,6 @@ void Foam::domainDecomposition::postReadUpdateReconstruct
     const fvMesh::readUpdateState stat
 )
 {
-    if (stat >= fvMesh::TOPO_CHANGE)
-    {
-        completeMesh_->postConstruct(false, fvMesh::stitchType::none);
-    }
-
     if (completeMesh_->stitcher().stitches() && stat != fvMesh::UNCHANGED)
     {
         procFaceAddressingBf_.clear();

src/parallel/parallel/domainDecompositionDecompose.C

@@ -930,6 +930,9 @@ void Foam::domainDecomposition::decompose()
         // Add patches to the mesh
         procMesh.addFvPatches(procPatches);
 
+        // Allocate a stitcher, but do not stitch
+        procMesh.postConstruct(false, fvMesh::stitchType::none);
+
         // Create point zones
         {
             const meshPointZones& pz = completeMesh().pointZones();

src/parallel/parallel/domainDecompositionReconstruct.C

@@ -2,7 +2,7 @@
   =========                 |
   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
    \\    /   O peration     | Website:  https://openfoam.org
-    \\  /    A nd           | Copyright (C) 2011-2023 OpenFOAM Foundation
+    \\  /    A nd           | Copyright (C) 2011-2024 OpenFOAM Foundation
      \\/     M anipulation  |
 -------------------------------------------------------------------------------
 License
@@ -516,6 +516,9 @@ void Foam::domainDecomposition::reconstruct()
     completeMesh_->setInstance(procMeshes()[0].facesInstance());
     completeMesh_->setPointsInstance(procMeshes()[0].pointsInstance());
 
+    // Allocate a stitcher, but do not stitch
+    completeMesh_->postConstruct(false, fvMesh::stitchType::none);
+
     Info<< decrIndent;
 }