From e7fbb7e55749c5ebaf526d287c38282c1957f49d Mon Sep 17 00:00:00 2001
From: mattijs <mattijs>
Date: Fri, 20 May 2022 17:35:04 +0100
Subject: [PATCH] BUG: snappyHexMesh: unset moving flag. Fixes #2477

We want to use movePoints to update all the geometry
but do not want to have a meshPhi
---
 .../meshRefinement/meshRefinementProblemCells.C     | 13 +++++++++++++
 .../snappyHexMeshDriver/snappyLayerDriver.C         |  6 ++++++
 .../snappyHexMeshDriver/snappyRefineDriver.C        |  6 ++++++
 3 files changed, 25 insertions(+)

diff --git a/src/mesh/snappyHexMesh/meshRefinement/meshRefinementProblemCells.C b/src/mesh/snappyHexMesh/meshRefinement/meshRefinementProblemCells.C
index fa541ba8a6..3581ac263c 100644
--- a/src/mesh/snappyHexMesh/meshRefinement/meshRefinementProblemCells.C
+++ b/src/mesh/snappyHexMesh/meshRefinement/meshRefinementProblemCells.C
@@ -606,6 +606,10 @@ void Foam::meshRefinement::markFacesOnProblemCells
             const_cast<Time&>(mesh_.time())++;
             pointField oldPoints(mesh_.points());
             mesh_.movePoints(newPoints);
+
+            // Unset any moving state
+            mesh_.moving(false);
+
             Pout<< "Writing newPoints mesh to time " << timeName()
                 << endl;
             write
@@ -615,6 +619,9 @@ void Foam::meshRefinement::markFacesOnProblemCells
                 mesh_.time().path()/"newPoints"
             );
             mesh_.movePoints(oldPoints);
+
+            // Unset any moving state
+            mesh_.moving(false);
         }
     }
 
@@ -1122,6 +1129,9 @@ void Foam::meshRefinement::markFacesOnProblemCellsGeometric
         );
 
         mesh_.movePoints(newPoints);
+
+        // Unset any moving state
+        mesh_.moving(false);
     }
 
 
@@ -1259,6 +1269,9 @@ void Foam::meshRefinement::markFacesOnProblemCellsGeometric
     // Restore points.
     mesh_.movePoints(oldPoints);
 
+    // Unset any moving state
+    mesh_.moving(false);
+
 
     Info<< "markFacesOnProblemCellsGeometric : marked "
         << returnReduce(nBaffleFaces, sumOp<label>())
diff --git a/src/mesh/snappyHexMesh/snappyHexMeshDriver/snappyLayerDriver.C b/src/mesh/snappyHexMesh/snappyHexMeshDriver/snappyLayerDriver.C
index e8758feb26..10cb92283f 100644
--- a/src/mesh/snappyHexMesh/snappyHexMeshDriver/snappyLayerDriver.C
+++ b/src/mesh/snappyHexMesh/snappyHexMeshDriver/snappyLayerDriver.C
@@ -4086,6 +4086,9 @@ void Foam::snappyLayerDriver::addLayers
 
             pp.movePoints(mesh.points());
 
+            // Unset any moving state
+            mesh.moving(false);
+
             // Update patchDisp (since not all might have been honoured)
             patchDisp = oldPatchPos - pp.localPoints();
         }
@@ -4390,6 +4393,9 @@ void Foam::snappyLayerDriver::addLayers
         pp.movePoints(mesh.points());
         medialAxisMoverPtr().movePoints(mesh.points());
 
+        // Unset any moving state
+        mesh.moving(false);
+
         // Grow out region of non-extrusion
         for (label i = 0; i < layerParams.nGrow(); i++)
         {
diff --git a/src/mesh/snappyHexMesh/snappyHexMeshDriver/snappyRefineDriver.C b/src/mesh/snappyHexMesh/snappyHexMeshDriver/snappyRefineDriver.C
index a33f6f8965..cfdf5597d9 100644
--- a/src/mesh/snappyHexMesh/snappyHexMeshDriver/snappyRefineDriver.C
+++ b/src/mesh/snappyHexMesh/snappyHexMeshDriver/snappyRefineDriver.C
@@ -2586,6 +2586,9 @@ Foam::label Foam::snappyRefineDriver::directionalSmooth
             );
             baseMesh.movePoints(baseMesh.points()+disp);
 
+            // Reset any moving state
+            baseMesh.moving(false);
+
             if (debug&meshRefinement::MESH)
             {
                 const_cast<Time&>(baseMesh.time())++;
@@ -2745,6 +2748,9 @@ Foam::label Foam::snappyRefineDriver::directionalSmooth
             );
             baseMesh.movePoints(baseMesh.points()+dispSmooth);
 
+            // Reset any moving state
+            baseMesh.moving(false);
+
             if (debug&meshRefinement::MESH)
             {
                 const_cast<Time&>(baseMesh.time())++;