<div dir="ltr"><div dir="ltr">On Sat, Feb 15, 2025 at 11:40 AM Aldo Bonfiglioli <<a href="mailto:aldo.bonfiglioli@unibas.it">aldo.bonfiglioli@unibas.it</a>> wrote:</div><div class="gmail_quote gmail_quote_container"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Dear all,<br>
<br>
I am trying to identify the boundary vertices that belong to a given <br>
stratum of the Face sets. This is going to be used to prescribe <br>
Dirichlet-type bcs.<br></blockquote><div><br></div><div>The label "Face Sets" is designed to only contain faces.</div><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
I can select the boundary faces (edges in 2D) of a given stratum using <br>
PetscCall(DMLabelGetStratumIS(label, stratum, user%bndryfaces(i), ierr)) <br>
and that looks ok;<br></blockquote><div><br></div><div>Yes, you get the faces marked with some BC value.</div><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
I then try to identify the boundary vertices by looping over the points <br>
(edges in 2D/faces in 3D) of the face set of a given stratum and <br>
retrieve the vertices that make up each individual edge/face.<br></blockquote><div><br></div><div>You can do that. You can also do something like</div><div><br></div><div>  DMLabelDuplicate(faceSets, &newLabel);  </div><div>  DMPlexLabelComplete(dm, newLabel);</div><div><br></div><div>which will put in all the points in the transitive closure (such as vertices).</div><div>Then you can just loop over the points in the label, and check for vertices</div><div>using DMPlexGetPointDepth().</div><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
For reasons I fail to understand, the above procedure fails to identify <br>
certain vertices (those circled in the enclosed pdf where different <br>
colours mark different ranks) in a parallel environment.<br></blockquote><div><br></div><div>I do this all the time, so it should not happen. If the above fails, can you</div><div>send a small reproducer?</div><div><br></div><div>  Thanks,</div><div><br></div><div>     Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
Questions:<br>
<br>
1. is there an available function that does what I am trying to do? I <br>
know that the boundary points can be found in the "marker" label, but I <br>
need to discriminate among Face Sets of different strata.<br>
<br>
2. what might be wrong in the aforementioned approach?<br>
<br>
Thanks,<br>
<br>
Aldo<br>
<br>
-- <br>
Dr. Aldo Bonfiglioli<br>
Associate professor of Fluid Machines<br>
Scuola di Ingegneria<br>
Universita' della Basilicata<br>
V.le dell'Ateneo lucano, 10 85100 Potenza ITALY<br>
tel:+39.0971.205203 fax:+39.0971.205215<br>
web: <a href="https://urldefense.us/v3/__http://docenti.unibas.it/site/home/docente.html?m=002423__;!!G_uCfscf7eWS!YLL7mzBlK9H_sEdtXob1AVklf8cVLhT3NTDdExIZsdI3xMfNHBoLXr92BmYzCXOuJqtm2L6OU4DJKV_81AkjvQf-Kra-Ym8P0tU$" rel="noreferrer" target="_blank">https://urldefense.us/v3/__http://docenti.unibas.it/site/home/docente.html?m=002423__;!!G_uCfscf7eWS!YLL7mzBlK9H_sEdtXob1AVklf8cVLhT3NTDdExIZsdI3xMfNHBoLXr92BmYzCXOuJqtm2L6OU4DJKV_81AkjvQf-Kra-Ym8P0tU$</a> <br>
</blockquote></div><div><br clear="all"></div><div><br></div><span class="gmail_signature_prefix">-- </span><br><div dir="ltr" class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div>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><br></div><div><a href="https://urldefense.us/v3/__http://www.cse.buffalo.edu/*knepley/__;fg!!G_uCfscf7eWS!cwtigEBvRpcpoLeV3aRNqP6fHwoAHN_2SW4Rjbh_ZqKelJ54Nvgncprg0QimeBjNfY5Ox-8qG6AzP4FxXi_H$" target="_blank">https://www.cse.buffalo.edu/~knepley/</a><br></div></div></div></div></div></div></div></div>