diff --git a/libsrc/meshing/basegeom.cpp b/libsrc/meshing/basegeom.cpp
index 0e6900e2..d1dfb8e3 100644
--- a/libsrc/meshing/basegeom.cpp
+++ b/libsrc/meshing/basegeom.cpp
@@ -1271,9 +1271,6 @@ namespace netgen
         if (res != MESHING3_OK) return 1;
         if (multithread.terminate) return 0;
 
-        RemoveIllegalElements (*mesh);
-        if (multithread.terminate) return 0;
-
         MeshQuality3d (*mesh);
       }
 
diff --git a/libsrc/meshing/meshfunc.cpp b/libsrc/meshing/meshfunc.cpp
index 4ddc2844..451a44cb 100644
--- a/libsrc/meshing/meshfunc.cpp
+++ b/libsrc/meshing/meshfunc.cpp
@@ -515,6 +515,7 @@ namespace netgen
            throw NgException ("Stop meshing since surface mesh not consistent");
         }
      }
+     RemoveIllegalElements (mesh);
   }
 
   void MergeMeshes( Mesh & mesh, Array<MeshingData> & md )