diff --git a/libsrc/meshing/meshfunc2d.cpp b/libsrc/meshing/meshfunc2d.cpp
index 90ac9b68..88955a36 100644
--- a/libsrc/meshing/meshfunc2d.cpp
+++ b/libsrc/meshing/meshfunc2d.cpp
@@ -61,6 +61,7 @@ namespace netgen
 	      cerr << "Optimization code " << optstr[j-1] << " not defined" << endl;
 	    }  
 	}
+    mesh.Compress(); // better: compress in individual steps, if necessary
     if (secondorder)
       {
         mesh.GetGeometry()->GetRefinement().MakeSecondOrder(mesh);