<html dir="ltr">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<style type="text/css" id="owaParaStyle"></style>
</head>
<body fpstyle="1" ocsi="0">
<div style="direction: ltr;font-family: Tahoma;color: #000000;font-size: 10pt;">Dear Matt,
<div><br>
</div>
<div>Thanks for your reply and help. I replaced the code with DMViewFromOptions() as you suggested. Attached you find the code for your reference, and the output of -dm_view ::ascii_info_detail.</div>
<div><br>
</div>
<div>a) -dm_view of the original mesh</div>
<div><br>
</div>
<div><style type="text/css">
p.p1 {margin: 0.0px 0.0px 0.0px 0.0px; font: 11.0px Menlo}
span.s1 {font-variant-ligatures: no-common-ligatures}
</style>
<p class="p1"><span class="s1"><b>DM Object: 1 MPI processes</b></span></p>
<p class="p1"><span class="s1"><b> type: plex</b></span></p>
<p class="p1"><span class="s1"><b>DM_0x84000000_0 in 3 dimensions:</b></span></p>
<p class="p1"><span class="s1"><b> 0-cells: 16</b></span></p>
<p class="p1"><span class="s1"><b> 1-cells: 61</b></span></p>
<p class="p1"><span class="s1"><b> 2-cells: 80</b></span></p>
<p class="p1"><span class="s1"><b> 3-cells: 34</b></span></p>
<p class="p1"><span class="s1"><b>Labels:</b></span></p>
<p class="p1"><span class="s1"><b> Cell Sets: 2 strata of sizes (17, 17)</b></span></p>
<p class="p1"><span class="s1"><b> Face Sets: 3 strata of sizes (12, 8, 10)</b></span></p>
<p class="p1"><span class="s1"><b> depth: 4 strata of sizes (16, 61, 80, 34)</b></span></p>
<p class="p1"><span class="s1"><br>
</span></p>
<p class="p1"><span class="s1"><font face="Tahoma" size="2">b) -dm_view of the sub-mesh 45 obtained by</font></span><span style="font-size: small; font-variant-ligatures: no-common-ligatures; font-family: Tahoma;"> DMPlexFilter</span></p>
<style type="text/css">
p.p1 {margin: 0.0px 0.0px 0.0px 0.0px; font: 11.0px Menlo}
span.s1 {font-variant-ligatures: no-common-ligatures}
span.s2 {font-variant-ligatures: no-common-ligatures; color: #3d1d81}
span.s3 {font-variant-ligatures: no-common-ligatures; color: #272ad8}
span.s4 {font-variant-ligatures: no-common-ligatures; color: #78492a}
</style></div>
<div><br>
</div>
<div><style type="text/css">
p.p1 {margin: 0.0px 0.0px 0.0px 0.0px; font: 11.0px Menlo}
span.s1 {font-variant-ligatures: no-common-ligatures}
</style>
<p class="p1"><span class="s1"><b> sub-mesh 45 </b></span></p>
<p class="p1"><span class="s1"><b>DM Object: 1 MPI processes</b></span></p>
<p class="p1"><span class="s1"><b> type: plex</b></span></p>
<p class="p1"><span class="s1"><b>DM_0x84000000_1 in 3 dimensions:</b></span></p>
<p class="p1"><span class="s1"><b> 0-cells: 12</b></span></p>
<p class="p1"><span class="s1"><b> 1-cells: 38</b></span></p>
<p class="p1"><span class="s1"><b> 2-cells: 44</b></span></p>
<p class="p1"><span class="s1"><b> 3-cells: 17</b></span></p>
<p class="p1"><span class="s1"><b>Labels:</b></span></p>
<p class="p1"><span class="s1"><b> depth: 4 strata of sizes (12, 38, 44, 17)</b></span></p>
</div>
<div><br>
</div>
<div>c) <span class="s1" style="font-family: Menlo; font-size: 11px;"><font face="Tahoma" size="2">-dm_view of the sub-mesh 43 obtained by DMPlexFilter:</font></span></div>
<div><br>
</div>
<div><style type="text/css">
p.p1 {margin: 0.0px 0.0px 0.0px 0.0px; font: 11.0px Menlo}
span.s1 {font-variant-ligatures: no-common-ligatures}
</style>
<p class="p1"><span class="s1"><b> sub-mesh 43 </b></span></p>
<p class="p1"><span class="s1"><b>DM Object: 1 MPI processes</b></span></p>
<p class="p1"><span class="s1"><b> type: plex</b></span></p>
<p class="p1"><span class="s1"><b>DM_0x84000000_2 in 3 dimensions:</b></span></p>
<p class="p1"><span class="s1"><b> 0-cells: 12</b></span></p>
<p class="p1"><span class="s1"><b> 1-cells: 38</b></span></p>
<p class="p1"><span class="s1"><b> 2-cells: 44</b></span></p>
<p class="p1"><span class="s1"><b> 3-cells: 17</b></span></p>
<p class="p1"><span class="s1"><b>Labels:</b></span></p>
<p class="p1"><span class="s1"><b> depth: 4 strata of sizes (12, 38, 44, 17)</b></span></p>
</div>
<div><br>
</div>
<div>d) -dm_view of the "sub-mesh 40" (interface between sub-meshes 43 and 45) obtained by DMPlexCreateSubMesh:</div>
<div><br>
</div>
<div><style type="text/css">
p.p1 {margin: 0.0px 0.0px 0.0px 0.0px; font: 11.0px Menlo}
span.s1 {font-variant-ligatures: no-common-ligatures}
</style>
<p class="p1"><span class="s1"><b> interface 40 </b></span></p>
<p class="p1"><span class="s1"><b>DM Object: 1 MPI processes</b></span></p>
<p class="p1"><span class="s1"><b> type: plex</b></span></p>
<p class="p1"><span class="s1"><b>DM_0x84000000_3 in 2 dimensions:</b></span></p>
<p class="p1"><span class="s1"><b> 0-cells: 0</b></span></p>
<p class="p1"><span class="s1"><b> 1-cells: 0</b></span></p>
<p class="p1"><span class="s1"><b> 2-cells: 0 (-83895904)</b></span></p>
<p class="p1"><span class="s1"><b>Labels:</b></span></p>
<p class="p1"><span class="s1"><b> depth: 0 strata of sizes ()</b></span></p>
</div>
<div><br>
</div>
<div>Questions:</div>
<div><br>
</div>
<div>1) the output of (d) is empty. What am I missing?</div>
<div><br>
</div>
<div>2) I used the "Cell Sets" DMLabel along with value=45 or 43 to yield the output of (b) and (c), respectively, and they do not contain "Face Sets". Is it expected?</div>
<div><br>
</div>
<div>I will address the other question in a separate e-mail, because it is not related to the issues I described above. My apologies the confusion.</div>
<div><br>
</div>
<div>Regards,</div>
<div> Frederico.</div>
<div><br>
<div style="font-family: Times New Roman; color: #000000; font-size: 16px">
<hr tabindex="-1">
<div id="divRpF760684" style="direction: ltr;"><font face="Tahoma" size="2" color="#000000"><b>From:</b> Matthew Knepley [knepley@gmail.com]<br>
<b>Sent:</b> 26 February 2017 23:26<br>
<b>To:</b> Santos Teixeira Frederico<br>
<b>Cc:</b> petsc-users@mcs.anl.gov<br>
<b>Subject:</b> Re: [petsc-users] Multi-domain meshes with DMPLEX<br>
</font><br>
</div>
<div></div>
<div>
<div dir="ltr">
<div class="gmail_extra">
<div class="gmail_quote">On Sun, Feb 26, 2017 at 10:51 AM, Santos Teixeira Frederico
<span dir="ltr"><<a href="mailto:fsantost@student.ethz.ch" target="_blank">fsantost@student.ethz.ch</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex; border-left:1px #ccc solid; padding-left:1ex">
<div>
<div style="direction:ltr; font-family:Tahoma; color:#000000; font-size:10pt">Hi folks,
<div><br>
</div>
<div>Attached you find small.msh, a mesh with two prisms (physical domains 43 and 45) and the interface between them (physical domain 40).</div>
<div>test_dmplex.cpp loads the mesh and tries to filter each sub-mesh and the interface. I use PETSc master.</div>
<div><br>
</div>
<div><span style="font-size:13.3333px">When I load small.msh, all physical domains are correctly assigned to</span><span style="font-size:13.3333px"> "Cell Sets" or "Face Sets"</span><span style="font-size:10pt">. </span></div>
<div><br>
</div>
<div>1) Both sub-meshes (prisms with domains 43 and 45) are filtered with DMPlexFilter and saved correctly, but the interface isn't. What am I missing?</div>
</div>
</div>
</blockquote>
<div><br>
</div>
<div>What exactly does "the interface isn't" mean? Does it mean you are not getting what you expect from the submesh?</div>
<div><br>
</div>
<div>One way to more easily communicate about this is to replace your blocks of Viewer commands with one line</div>
<div>DMViewFromOptions() calls. That way you can show me just -dm_view or -dm_view ::ascii_info_detail, or output</div>
<div>VTK with -dm_view vtk:mesh.vtu, or HDF5 -dm_view hdf5:mesh.h5, etc.</div>
<div><br>
</div>
<div>We could start off with sending me just the plain ASCII output for the mesh and submeshes. I do recognize you sent</div>
<div>the code, but I have a conference deadline for Monday (I will be there until Friday), so I am unlikely to build it before them.</div>
<div> </div>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex; border-left:1px #ccc solid; padding-left:1ex">
<div>
<div style="direction:ltr; font-family:Tahoma; color:#000000; font-size:10pt">
<div>2) <span style="font-size:10pt">I can retrieve any specific strata of "Cell Sets" with DMPlexFilter, but this operation does not include any related strata of "Face Sets", as pointed by DMView. How can I do it? </span></div>
</div>
</div>
</blockquote>
<div><br>
</div>
<div>Do DMFilter is not propagating labels? Okay, that can go on the bug list. Is this what you mean?</div>
<div> </div>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex; border-left:1px #ccc solid; padding-left:1ex">
<div>
<div style="direction:ltr; font-family:Tahoma; color:#000000; font-size:10pt">
<div><span style="font-size:10pt">3) Physical domains 43 and 45 could be interpreted as fluid and solid domains in a F.S.I. simulation, for example, and I would like to define fluid and solid sub-problems only on the respective sub-mesh. </span><span style="font-size:10pt">Do
I really </span><span style="font-size:10pt">need</span><span style="font-size:10pt"> to have fully operational sub-DMPLEX's in order to achieve it? </span><span style="font-size:10pt">Is there an easier way?</span></div>
</div>
</div>
</blockquote>
<div><br>
</div>
<div>I am not sure I understand what is going wrong yet.</div>
<div><br>
</div>
<div> Thanks,</div>
<div><br>
</div>
<div> Matt</div>
<div> </div>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex; border-left:1px #ccc solid; padding-left:1ex">
<div>
<div style="direction:ltr; font-family:Tahoma; color:#000000; font-size:10pt">
<div><span style="font-size:10pt">My apologies if the explanations are confusing!</span></div>
<div><span style="font-size:10pt"><br>
</span></div>
<div>Best regards,</div>
<div> Frederico.</div>
</div>
</div>
</blockquote>
</div>
<br>
<br clear="all">
<div><br>
</div>
-- <br>
<div class="gmail_signature">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>
</div>
</div>
</div>
</div>
</div>
</body>
</html>