From a7a50678309d0b047c5bb150b7ebdc8f4cb07fd9 Mon Sep 17 00:00:00 2001
From: Matthias Hochsteger <mhochsteger@cerbsim.com>
Date: Thu, 14 Apr 2022 17:40:22 +0200
Subject: [PATCH] Call FindOpenElements before each optimization step

---
 libsrc/meshing/meshfunc.cpp | 5 +++--
 1 file changed, 3 insertions(+), 2 deletions(-)

diff --git a/libsrc/meshing/meshfunc.cpp b/libsrc/meshing/meshfunc.cpp
index 553a1422..01e2a832 100644
--- a/libsrc/meshing/meshfunc.cpp
+++ b/libsrc/meshing/meshfunc.cpp
@@ -451,6 +451,7 @@ namespace netgen
             const char * optstr = "mcmstmcmstmcmstmcm";
             for (size_t j = 1; j <= strlen(optstr); j++)
             {
+               mesh.FindOpenElements();
                mesh.CalcSurfacesOfNode();
                mesh.FreeOpenElementsEnvironment(2);
                mesh.CalcSurfacesOfNode();
@@ -466,12 +467,12 @@ namespace netgen
 
             }
 
-            mesh.FindOpenElements();	      
+            mesh.FindOpenElements(domain);
             PrintMessage (3, "Call remove problem");
             // mesh.Save("before_remove.vol");
             RemoveProblem (mesh, domain);
             // mesh.Save("after_remove.vol");
-            mesh.FindOpenElements();
+            mesh.FindOpenElements(domain);
          }
          else
            {