From 12adfe13fa0a67bc2a0a2001de1814783e66bd1e Mon Sep 17 00:00:00 2001
From: Mark Olesen <Mark.Olesen@esi-group.com>
Date: Fri, 1 Mar 2019 14:50:54 +0100
Subject: [PATCH] ENH: update empty sampled surfaces storage (#1217)

- also store empty sampled surfaces, otherwise we miss geometry updates.
---
 .../sampledSurfaces/sampledSurfaces.C         | 49 ++++++++-----------
 src/surfMesh/polySurface/polySurface.H        |  4 --
 2 files changed, 21 insertions(+), 32 deletions(-)

diff --git a/src/sampling/sampledSurface/sampledSurfaces/sampledSurfaces.C b/src/sampling/sampledSurface/sampledSurfaces/sampledSurfaces.C
index a5682f042f..efe043f428 100644
--- a/src/sampling/sampledSurface/sampledSurfaces/sampledSurfaces.C
+++ b/src/sampling/sampledSurface/sampledSurfaces/sampledSurfaces.C
@@ -490,8 +490,7 @@ bool Foam::sampledSurfaces::read(const dictionary& dict)
 
 bool Foam::sampledSurfaces::performAction(unsigned request)
 {
-    // Update surfaces, writer associations etc.
-
+    // Update surfaces and store
     bool ok = false;
 
     forAll(*this, surfi)
@@ -508,12 +507,27 @@ bool Foam::sampledSurfaces::performAction(unsigned request)
             nFaces_[surfi] = returnReduce(s.faces().size(), sumOp<label>());
 
             ok = ok || nFaces_[surfi];
+
+
+            // Store surfaces (even empty ones) otherwise we miss geometry
+            // updates.
+            // Any associated fields will be removed if the size changes
+
+            if ((request & actions_[surfi]) & ACTION_STORE)
+            {
+                storeRegistrySurface(s);
+            }
+
+            if ((request & actions_[surfi]) & ACTION_SURF_MESH)
+            {
+                s.storeSurfMesh();
+            }
         }
     }
 
     if (!ok)
     {
-        // No surface with faces or an applicable action
+        // No surface with an applicable action or with faces to sample
         return true;
     }
 
@@ -522,36 +536,13 @@ bool Foam::sampledSurfaces::performAction(unsigned request)
     // Only seems to be needed for VTK legacy
     countFields();
 
+    // Update writers
 
     forAll(*this, surfi)
     {
         const sampledSurface& s = (*this)[surfi];
 
-        if (!(request & actions_[surfi]))
-        {
-            continue;
-        }
-
-        // TDB: do we store empty surfaces, skip them, or remove them
-        // from the registry?
-        // - For now, just skip touching them.
-
-        if (!nFaces_[surfi])
-        {
-            continue;
-        }
-
-        if ((request & actions_[surfi]) & ACTION_STORE)
-        {
-            storeRegistrySurface(s);
-        }
-
-        if ((request & actions_[surfi]) & ACTION_SURF_MESH)
-        {
-            s.storeSurfMesh();
-        }
-
-        if ((request & actions_[surfi]) & ACTION_WRITE)
+        if (((request & actions_[surfi]) & ACTION_WRITE) && nFaces_[surfi])
         {
             // Output writers
             surfaceWriter& outWriter = writers_[surfi];
@@ -588,6 +579,8 @@ bool Foam::sampledSurfaces::performAction(unsigned request)
         }
     }
 
+    // Sample fields
+
     const IOobjectList objects(obr_, obr_.time().timeName());
 
     performAction<volScalarField>(objects, request);
diff --git a/src/surfMesh/polySurface/polySurface.H b/src/surfMesh/polySurface/polySurface.H
index a7e51f946f..5956f5e187 100644
--- a/src/surfMesh/polySurface/polySurface.H
+++ b/src/surfMesh/polySurface/polySurface.H
@@ -140,10 +140,6 @@ protected:
         }
 
 
-        //- Update with new contents
-        void copyContents(const MeshedSurface<face>& surf, bool validate=false);
-
-
 public:
 
     // Public Typedefs