<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>Dear developers,</p>

<p><br>

</p>

<p>I have a mesh with topological dimension 2 (with triangular faces) embedded in a 3D space stored in a hdf5 file. This file contains the coordinates (x,y,z) of each vertex and also, it contains a list of vertices that make each element. I was able to read

the mesh into a DMPlex using DMPlexCreateFromFile (without interpolation, i.e. interpolate=PETSC_FALSE ). When I do this, the elements are stored at stratum height 0 and the vertices are stored at stratum height 1 (as checked via DMPlexGetHeightStratum).

 However, when I interpolate the mesh (<span><span>interpolate</span>=PETSC_TRUE</span>), it seems that some entities are created in stratum height 2 and 3. From my understanding of DMPlex, this should not happen for 2D elements because there are edges but

 no 3D cells. <br>

</span></p>

<p><span><br>

</span></p>

<p><span>Does DMPlex assume that the mesh in 3D has a topological dimension 3 and is it possible to use 2D mesh embeded in 3D space?</span></p>

<p><span><br>

</span></p>

<p><span>Best regards<br>

</span></p>

<p><br>

</p>

<div id="Signature">

<div class="BodyFragment"><font size="2">

<div class="PlainText">François Fillion-Gourdeau, PhD<br>

INRS-EMT</div>

</font></div>

</div>

</div>

</body>

</html>