<div dir="ltr">Dear Matthew,<div><br></div><div>You are (again) right. This Gmsh test file is "dirty" : some triangles do not belong to tests. Sorry for that.</div><div>I have tried with another geo file (which is clean in that sense) and PETSc reads with no error.</div><div><br></div><div>I take this opportunity to ask my initial question : are the labels of the Physical Entities saved somewhere ? <br></div><div>A gmsh test file with several Physical entities is attached to this email. </div><div><br></div><div>Thank you again,</div><div>Nicolas</div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">Le mer. 19 mai 2021 à 15:12, 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">Okay, what Plex is saying is that Gmsh gave it a facet (a set of vertices), but those are not the vertices<div>of any tetrahedral face in the mesh. It is hard to verify such a thing with 12K tets. I have definitely read in</div><div>GMsh files before with facets, so I don't think the code is completely wrong.</div><div><br></div><div>Would it be possible to make a smaller mesh (10-20 cells) that shows the same bug?</div><div><br></div><div> Thanks,</div><div><br></div><div> Matt</div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Tue, May 18, 2021 at 12:25 PM Karin&NiKo <<a href="mailto:niko.karin@gmail.com" target="_blank">niko.karin@gmail.com</a>> wrote:<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">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" target="_blank">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>
</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>
</blockquote></div>