missing call to setAllVolumesPositive

Geo/GModel.cpp

@@ -540,13 +540,11 @@ int GModel::mesh(int dimension)
 
 bool GModel::setAllVolumesPositive()
 {
-  std::vector<GEntity*> entities;
-  getEntities(entities);
   bool ok = true;
-  for(std::set<GRegion*,GEntityLessThan>::iterator itReg=regions.begin();
-                                          itReg != regions.end(); ++itReg) {
-    int nbEl = (*itReg)->getNumMeshElements();
-    for (int iEl=0; iEl<nbEl; ++iEl) ok = ok && (*itReg)->getMeshElement(iEl)->setVolumePositive();
+  for(riter it = regions.begin(); it != regions.end(); ++it) {
+    int nbEl = (*it)->getNumMeshElements();
+    for (int iEl = 0; iEl < nbEl; ++iEl)
+      ok = (ok && (*it)->getMeshElement(iEl)->setVolumePositive());
   }
   return ok;
 }

Mesh/Generator.cpp

@@ -671,6 +671,8 @@ static void Mesh3D(GModel *m)
     }
   }
 
+  m->setAllVolumesPositive();
+
   double t2 = Cpu();
   CTX::instance()->meshTimer[2] = t2 - t1;
   Msg::StatusBar(true, "Done meshing 3D (%g s)", CTX::instance()->meshTimer[2]);