<div dir="ltr">Thanks, I think I understand. <div>So at the moment, the user has to manage the link between the labels of the different strata and the names. </div><div>I am mainly willing to use DMPlex for AMR and names of groups are essential in our code - reading what Mat says </div><div>gives me the impression that this is doable. </div><div><br></div><div>Thanks again,</div><div>Nicolas</div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">Le mer. 19 mai 2021 à 16:02, 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 Wed, May 19, 2021 at 9:57 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">Thank you very much Lawrence.<div>Are the names lost or are they saved somewhere ?</div><div>If I do : "./ex2 -filename /tmp/test/Cube_with_facets.msh4 -dm_view vtk:/tmp/foo.vtk" , I only get the tets of the initial mesh.</div></div></blockquote><div><br></div><div>Yes, the names are lost. Jed, Lisandro and Stefano have also asked me about this. I have not done since everything is currently done by</div><div>number in Plex. It would involve either</div><div><br></div><div>a) making a separate table in the DM which translates names to numbers</div><div><br></div><div>or</div><div><br></div><div>b) putting such a translation table in the DMLabel itself</div><div><br></div><div>I am putting it off because doing some short is not hard, but getting it to work with everything else would be somewhat of a pain.</div><div>It would need to distribute properly, propagate to refined/coarsened/subset meshes, etc. which labels currently do. It would not</div><div>be properly passed on when you just set an int on a point. I think it will happen eventually, but I don't think I have time in the next</div><div>few weeks.</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, </div><div>Nicolas</div><div><br></div></div><div class="gmail_quote"><div dir="ltr" class="gmail_attr">Le mer. 19 mai 2021 à 15:49, Lawrence Mitchell <<a href="mailto:wence@gmx.li" target="_blank">wence@gmx.li</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"><br>
<br>
> On 19 May 2021, at 14:46, Karin&NiKo <<a href="mailto:niko.karin@gmail.com" target="_blank">niko.karin@gmail.com</a>> wrote:<br>
> <br>
> Dear Matthew,<br>
> <br>
> You are (again) right. This Gmsh test file is "dirty" : some triangles do not belong to tests. Sorry for that.<br>
> I have tried with another geo file (which is clean in that sense) and PETSc reads with no error.<br>
> <br>
> I take this opportunity to ask my initial question : are the labels of the Physical Entities saved somewhere ? <br>
> A gmsh test file with several Physical entities is attached to this email. <br>
<br>
Yes, plex represents Physical entities as "labels". In particular, for facets, they are loaded as the "Face Sets" label, plex also loads markers on cells in the "Cell Sets" label.<br>
<br>
e.g. here's the DMView output of your mesh.<br>
<br>
DM Object: DM_0x84000000_0 1 MPI processes<br>
type: plex<br>
DM_0x84000000_0 in 3 dimensions:<br>
0-cells: 64<br>
1-cells: 279<br>
2-cells: 378<br>
3-cells: 162<br>
Labels:<br>
celltype: 4 strata with value/size (0 (64), 6 (162), 3 (378), 1 (279))<br>
depth: 4 strata with value/size (0 (64), 1 (279), 2 (378), 3 (162))<br>
Cell Sets: 1 strata with value/size (7 (162))<br>
Face Sets: 6 strata with value/size (3 (18), 6 (18), 2 (18), 5 (18), 1 (18), 4 (18))<br>
<br>
<br>
Lawrence</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></div>
</blockquote></div>