<div dir="ltr">Dear PETSc team,<br><br>I have tried to load a test mesh available in Gmsh' s demos directory  (share/doc/gmsh/demos/simple_geo/filter.geo, attached to this email) as a DMPlex.<br>So I produced a msh4 file by doing : <br>        gmsh -3 filter.geo -o /tmp/test.msh4<br>Then I used src/dm/impls/plex/tutorials/ex2.c to load the mesh by doing :<br>        ./ex2 -filename /tmp/test.msh4<br><br>Unfortunately I get the error :<div>  <br>[0]PETSC ERROR: --------------------- Error Message --------------------------------------------------------------<br>[0]PETSC ERROR: No support for this operation for this object type<br>[0]PETSC ERROR: Could not determine Plex facet for Gmsh element 1268 (Plex cell 12681)<br></div><div><br></div><div>The error seems to come from the fact that the msh file contains tets *and* facets *only on the Physical entities* (aka parts of the mesh boundary where </div><div>the user will assign Dirichlet or Neuman conditions). </div><div>If I suppress these facets by commenting  the "Physical Surface" lines in the geo file and regenerating the mesh, everything is fine.</div><div><br></div><div>But the use of these "Physical" stuff is very common in lots of finite element codes in order to assign boundary conditions. </div><div>How should I do to keep  these boundary groups of 2D elements (with corresponding names) ?</div><div><br></div><div>Thanks for your help,</div><div>Nicolas</div><div><br></div><div><br></div></div>