<div dir="ltr"><div dir="ltr">On Fri, May 5, 2023 at 10:55 AM Vilmer Dahlberg via petsc-users <<a href="mailto:petsc-users@mcs.anl.gov">petsc-users@mcs.anl.gov</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 class="msg-7327637652967978226">




<div dir="ltr">
<div id="m_-7327637652967978226divtagdefaultwrapper" style="font-size:12pt;color:rgb(0,0,0);font-family:Calibri,Helvetica,sans-serif" dir="ltr">
<p>Hi.</p>
<p><br>
</p>
<p>I'm trying to read a mesh of higher element order, in this example a mesh consisting of 10-node tetrahedral elements, from gmsh, into PETSC. But It looks like the mesh is not properly being loaded and converted into a DMPlex. gmsh tells me it has generated
 a mesh with 7087 nodes, but when I view my dm object it tells me it has 1081 0-cells. This is the printout I get</p></div></div></div></blockquote><div><br></div><div>Hi Vilmer,</div><div><br></div><div>Plex makes a distinction between topological entities, like vertices, edges and cells, and the function spaces used to represent fields, like velocity or coordinates. When formats use "nodes", they mix the two concepts together.</div><div><br></div><div>You see that if you add the number of vertices and edges, you get 7087, since for P2 there is a "node" on every edge. Is anything else wrong?</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 class="msg-7327637652967978226"><div dir="ltr"><div id="m_-7327637652967978226divtagdefaultwrapper" style="font-size:12pt;color:rgb(0,0,0);font-family:Calibri,Helvetica,sans-serif" dir="ltr">
<p>...</p>
<p></p>
<div><br>
</div>
<div>Info    : Done meshing order 2 (Wall 0.0169823s, CPU 0.016662s)<br>
</div>
<p></p>
<p></p>
<div>Info    : 7087 nodes 5838 elements</div>
<div><br>
</div>
<div>...</div>
<div><br>
</div>
<div>
<div>DM Object: DM_0x84000000_0 1 MPI process<br>
  type: plex<br>
DM_0x84000000_0 in 3 dimensions:<br>
  Number of 0-cells per rank: 1081<br>
  Number of 1-cells per rank: 6006<br>
  Number of 2-cells per rank: 9104<br>
  Number of 3-cells per rank: 4178<br>
Labels:<br>
  celltype: 4 strata with value/size (0 (1081), 6 (4178), 3 (9104), 1 (6006))<br>
  depth: 4 strata with value/size (0 (1081), 1 (6006), 2 (9104), 3 (4178))<br>
  Cell Sets: 1 strata with value/size (2 (4178))<br>
  Face Sets: 6 strata with value/size (12 (190), 21 (242), 20 (242), 11 (192), 22 (242), 10 (188))<br>
Field P2:<br>
  adjacency FEM<br>
...<br>
</div>
</div>
<p></p>
<p><br>
</p>
<p>To replicate the error try generating a mesh according to</p>
<p><br>
</p>
<p><a href="https://gmsh.info/doc/texinfo/gmsh.html#t5" id="m_-7327637652967978226LPlnk379974" target="_blank">https://gmsh.info/doc/texinfo/gmsh.html#t5</a></p>
<p><a href="https://gmsh.info/doc/texinfo/gmsh.html#t5" id="m_-7327637652967978226LPlnk379974" target="_blank"><br>
</a></p>
<p>setting the element order to 2, and then loading the mesh using</p>
<p><br>
</p>
<p>DMPlexCreateGmshFromFile</p>
<p><br>
</p>
<p>I don't have any issues when i set the element order to 1.<br>
</p>
<p><br>
</p>
<p>Thanks in advance,</p>
<p>Vilmer<br>
</p>
</div>
</div>

</div></blockquote></div><br clear="all"><div><br></div><span class="gmail_signature_prefix">-- </span><br><div dir="ltr" class="gmail_signature"><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>