<!DOCTYPE html>
<html>
<head>
<meta http-equiv="content-type" content="text/html; charset=UTF-8">
</head>
<body>
<p>Dear all,</p>
<p>I am having troubles in using DMPlexExtractSubMesh to extract the
different strata of the Face Sets of a given mesh.</p>
<p>When run on the enclosed tetrahedral mesh of the unit cube
generated with gmsh</p>
<p>
<blockquote type="cite">
<span style="font-family:monospace"><span
style="color:#000000;background-color:#ffffff;">Face Sets: 6
strata with value/size (1 (246), 2 (246), 3 (246), 4 (246),
5 (242), 6 (242))</span><br>
<span style="color:#000000;background-color:#ffffff;">
</span><br>
</span>
</blockquote>
I would expect 246 "points" on stratum 3, but when I DMview the
subdm (and plot it) the surface mesh looks incomplete</p>
<p>
<blockquote type="cite">
<span style="font-family:monospace"><span
style="color:#000000;background-color:#ffffff;">DM Object:
patch_03 1 MPI process</span><span
style="color:#000000;background-color:#ffffff;">
</span><br>
<span style="color:#000000;background-color:#ffffff;"> type:
plex</span><span
style="color:#000000;background-color:#ffffff;">
</span><br>
<span style="color:#000000;background-color:#ffffff;">patch_03
in 2 dimensions:</span><span
style="color:#000000;background-color:#ffffff;">
</span><br>
<span style="color:#000000;background-color:#ffffff;"> Cells
are at height 1</span><span
style="color:#000000;background-color:#ffffff;">
</span><br>
<span style="color:#000000;background-color:#ffffff;"> Number
of 0-cells per rank: 122</span><span
style="color:#000000;background-color:#ffffff;">
</span><br>
<span style="color:#000000;background-color:#ffffff;"> Number
of 1-cells per rank: 325</span><span
style="color:#000000;background-color:#ffffff;">
</span><br>
<span style="color:#000000;background-color:#ffffff;">Number
of 2-cells per rank: 204</span><span
style="color:#000000;background-color:#ffffff;">
</span><br>
<span style="color:#000000;background-color:#ffffff;">Number
of 3-cells per rank: 204 [204]</span><span
style="color:#000000;background-color:#ffffff;">
</span><br>
<span style="color:#000000;background-color:#ffffff;">Labels:</span><span
style="color:#000000;background-color:#ffffff;">
</span><br>
<span style="color:#000000;background-color:#ffffff;">celltype:
4 strata with value/size (0 (122), 1 (325), 3 (204), 12
(204))</span><span
style="color:#000000;background-color:#ffffff;">
</span><br>
<span style="color:#000000;background-color:#ffffff;">depth:
4 strata with value/size (0 (122), 1 (325), 2 (204), 3
(204))</span><span
style="color:#000000;background-color:#ffffff;">
</span><br>
<span style="color:#000000;background-color:#ffffff;">Cell
Sets: 1 strata with value/size (1 (204))</span><span
style="color:#000000;background-color:#ffffff;">
</span><br>
<span style="color:#000000;background-color:#ffffff;">Face
Sets: 1 strata with value/size (3 (204))</span><span
style="color:#000000;background-color:#ffffff;">
</span><br>
<span style="color:#000000;background-color:#ffffff;">Edge
Sets: 2 strata with value/size (1 (8), 5 (8))</span><br>
<span style="color:#000000;background-color:#ffffff;">
</span><br>
</span>
</blockquote>
see also patch_03.pdf</p>
<p>What am I doing wrong?</p>
<p>A simple reproducer (compiles with petsc-3.24.0) and the gmsh
mesh are enclosed.</p>
<p>Thanks,</p>
<p>Aldo</p>
<pre class="moz-signature" cols="72">--
Dr. Aldo Bonfiglioli
Associate professor of Fluid Mechanics
Dipartimento di Ingegneria
Universita' della Basilicata
V.le dell'Ateneo Lucano, 10 85100 Potenza ITALY
tel:+39.0971.205203 fax:+39.0971.205215
web: <a class="moz-txt-link-freetext" href="https://urldefense.us/v3/__http://docenti.unibas.it/site/home/docente.html?m=002423__;!!G_uCfscf7eWS!aMKmGG4aim9XcbNSnDyHUkDyhUkQHGZ-u-xX2C-sycYUMmtTij6AwqsQbZPXJSvPp9KUfgwRJK2Ok6Me2BLgO0en1w4QF2fHo7s$">http://docenti.unibas.it/site/home/docente.html?m=002423</a></pre>
</body>
</html>