<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" data-smartmail="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>