testing new Gmsh partitioner behavior (interfaces might not belong to any...

testing new Gmsh partitioner behavior (interfaces might not belong to any physical group, so SaveAll is required)

examples/helmholtzflow/nacelle3D/main.cpp

@@ -101,8 +101,10 @@ int main(int argc, char **argv)
     if(!getMPIRank()) {
       // build geometry
       Nacelle(h_near, h_far, D, R, L_duct, 2.*h_far, 2.*h_near, nDom);
-      gmsh::option::setNumber("Mesh.ElementOrder", FEMorder); // FIXME
+      gmsh::option::setNumber("Mesh.ElementOrder", FEMorder);
       gmsh::option::setNumber("Mesh.SecondOrderLinear", 1);
+      // save all, as partitioning will create some entities without physical groups
+      gmsh::option::setNumber("Mesh.SaveAll", 1);
       gmsh::model::mesh::generate(3);
       gmsh::model::mesh::partition(nDom);
       // save partitioned mesh in single file for mono-process runs

tutorials/helmholtz/scattering3D/main.cpp

@@ -106,6 +106,8 @@ int main(int argc, char **argv)
       sphere(Rin, Rext, lc);
       // partition
       gmsh::model::mesh::partition(nDom);
+      // save all, as partitioning will create some entities without physical groups
+      gmsh::option::setNumber("Mesh.SaveAll", 1);
       // save partitioned mesh in single file for sequential runs
       gmsh::option::setNumber("Mesh.PartitionSplitMeshFiles", 0);
       gmsh::write("sphere.msh");