<div dir="ltr">Sure! I send you both, one with the facets and one without.<br><div><br></div><div>Thanks,</div><div>Nicolas</div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">Le mar. 18 mai 2021 à 17:49, Matthew Knepley <<a href="mailto:knepley@gmail.com">knepley@gmail.com</a>> a écrit :<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div dir="ltr">On Tue, May 18, 2021 at 8:18 AM Karin&NiKo <<a href="mailto:niko.karin@gmail.com" target="_blank">niko.karin@gmail.com</a>> wrote:<br></div><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><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></blockquote><div><br></div><div>Can you also send the *.msh file? I do not have Gmsh on this machine.</div><div><br></div><div>  Thanks,</div><div><br></div><div>     Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div>Thanks for your help,</div><div>Nicolas</div><div><br></div><div><br></div></div>
</blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>-- Norbert Wiener</div><div><br></div><div><a href="http://www.cse.buffalo.edu/~knepley/" target="_blank">https://www.cse.buffalo.edu/~knepley/</a><br></div></div></div></div></div></div></div></div>
</blockquote></div>