<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<style type="text/css" style="display:none;"><!-- P {margin-top:0;margin-bottom:0;} --></style>
</head>
<body dir="ltr">
<div id="divtagdefaultwrapper" style="font-size:12pt;color:#000000;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>
<p><br>
</p>
<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" class="OWAAutoLink" id="LPlnk379974">https://gmsh.info/doc/texinfo/gmsh.html#t5</a></p>
<p><a href="https://gmsh.info/doc/texinfo/gmsh.html#t5" class="OWAAutoLink" id="LPlnk379974"><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>
</body>
</html>