From e750e721a345f88377453c2af73c328b49093476 Mon Sep 17 00:00:00 2001
From: Christophe Geuzaine <cgeuzaine@uliege.be>
Date: Wed, 19 Jan 2022 21:58:46 +0200
Subject: [PATCH] testing new Gmsh partitioner behavior (interfaces might not
 belong to any physical group, so SaveAll is required)

---
 examples/helmholtzflow/nacelle3D/main.cpp | 4 +++-
 tutorials/helmholtz/scattering3D/main.cpp | 2 ++
 2 files changed, 5 insertions(+), 1 deletion(-)

diff --git a/examples/helmholtzflow/nacelle3D/main.cpp b/examples/helmholtzflow/nacelle3D/main.cpp
index 007ee2d6..df872089 100644
--- a/examples/helmholtzflow/nacelle3D/main.cpp
+++ b/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
diff --git a/tutorials/helmholtz/scattering3D/main.cpp b/tutorials/helmholtz/scattering3D/main.cpp
index df9d618f..008a284b 100644
--- a/tutorials/helmholtz/scattering3D/main.cpp
+++ b/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");
-- 
GitLab