diff --git a/libsrc/meshing/meshfunc.cpp b/libsrc/meshing/meshfunc.cpp
index ff2af770..0368a5d4 100644
--- a/libsrc/meshing/meshfunc.cpp
+++ b/libsrc/meshing/meshfunc.cpp
@@ -275,7 +275,6 @@ namespace netgen
     auto & mesh = *md.mesh;
     auto domain = md.domain;
     MeshingParameters & mp = md.mp;
-    cout << "try to close open quads " << md.domain << endl;
 
     int oldne;
     if (multithread.terminate)
@@ -283,7 +282,6 @@ namespace netgen
     
     mesh.CalcSurfacesOfNode();
     mesh.FindOpenElements(domain);
-    cout << "\thave close open elements " << mesh.GetNOpenElements() << endl;
     
     if (!mesh.GetNOpenElements())
       return;
@@ -294,7 +292,6 @@ namespace netgen
        
        if (mesh.HasOpenQuads())
          {
-           cout << "\thave close open quads" << endl;
            string rulefile = ngdir;
            
            const char ** rulep = NULL;
@@ -356,7 +353,6 @@ namespace netgen
              << "mesh has " << mesh.GetNE() << " prism/pyramid elements" << endl;
            
            mesh.FindOpenElements(domain);
-           cout << "\thave close open elements " << mesh.GetNOpenElements() << endl;
          }
      }
    
@@ -364,8 +360,8 @@ namespace netgen
    if (mesh.HasOpenQuads())
    {
       if(debugparam.write_mesh_on_error) {
-        md.mesh->Save("open_quads_"+ToString(md.domain)+".vol.gz");
-        GetOpenElements(*md.mesh, md.domain)->Save("open" + ToString(md.domain)+".vol");
+        md.mesh->Save("open_quads_starting_mesh_"+ToString(md.domain)+".vol.gz");
+        GetOpenElements(*md.mesh, md.domain)->Save("open_quads_rest_" + ToString(md.domain)+".vol.gz");
       }
       PrintSysError ("mesh has still open quads");
       throw NgException ("Stop meshing since too many attempts");
diff --git a/libsrc/meshing/python_mesh.cpp b/libsrc/meshing/python_mesh.cpp
index 298a2951..1a50a250 100644
--- a/libsrc/meshing/python_mesh.cpp
+++ b/libsrc/meshing/python_mesh.cpp
@@ -1298,7 +1298,7 @@ DLL_HEADER void ExportNetgenMeshing(py::module &m)
                CreateMPfromKwargs(mp, kwargs);
              }
              MeshVolume (mp, self);
-             // OptimizeVolume (mp, self);
+             OptimizeVolume (mp, self);
            }, py::arg("mp")=nullptr,
           meshingparameter_description.c_str(),
           py::call_guard<py::gil_scoped_release>())
diff --git a/libsrc/visualization/vsmesh.cpp b/libsrc/visualization/vsmesh.cpp
index 8fb523af..7e439d0f 100644
--- a/libsrc/visualization/vsmesh.cpp
+++ b/libsrc/visualization/vsmesh.cpp
@@ -456,7 +456,7 @@ namespace netgen
                 const auto& seg = (*mesh)[si];
                 Point<3> c = Center((*mesh)[seg[0]], (*mesh)[seg[1]]);
 		glRasterPos3d (c[0], c[1], c[2]);
-		snprintf (buf, size(buf),  "%d", int(seg.edgenr));
+		snprintf (buf, size(buf),  "%d", int(si));
 		MyOpenGLText (buf);
               }
             }
@@ -507,10 +507,7 @@ namespace netgen
                                (*mesh)[sel[2]],
                                (*mesh)[sel[3]]);
                   glRasterPos3d (c[0], c[1], c[2]);
-                  // auto & fd = mesh->GetFaceDescriptor(sel.GetIndex());
-                  // string s = ToString(fd.DomainIn()) + "/" + ToString(fd.DomainOut());
-                  // snprintf (buf, size(buf),  "%s", s.c_str());
-                  snprintf (buf, size(buf),  "%d", int(sel.GetIndex()));
+                  snprintf (buf, size(buf),  "%d", int(sei));
                   MyOpenGLText (buf);
                 }
             }