<div dir="ltr"><div dir="ltr">On Thu, Jan 31, 2019 at 6:22 AM rickcha--- via petsc-users <<a href="mailto:petsc-users@mcs.anl.gov">petsc-users@mcs.anl.gov</a>> wrote:<br></div><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Hello,<br>
<br>
I have an issue with DMGetStatumSize() and DMGetStratumIS().<br>
<br>
My intention is to access all vertices (height = 0) that are part of a named label in my DM.<br></blockquote><div><br></div><div>Depth = 0 are the vertices, but this works with DMPlexGetDpethStratum(). You seems like</div><div>you are calling DMLabelGetStratumIS(). The value here is the label value, not the height or depth.</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">
The problem is that only vertices that are exclusively part of this named label seem to be included. The label  is assigned to a face, let’s call it “A” for simplicities sake. Adjacent faces have their own label (and Dirichlet BC’s prescribed to them if that matters) which I will refer to as (“B”,”C”,”D” and “E”). Vertices that lie on the edge should be part of more than one stratum. I.e. belonging to “A” AND “B”, or “A” AND “B” AND “C” if they lie on a corner where the three faces meet.<br>
However, those shared vertices seem not to be included when calling DMGetStatumSize() and DMGetStratumIS().<br>
<br>
For a 5x5x5 hex cube with faces labeled “A”, “B”, […], I would expect DMGetStratumSize(dm,”A”,0,&size) to yield size = 36. But I get 16 instead, missing the 20 “shared” vertices.<br>
<br>
Is there a way to change this behaviour? Or do I have to get the face stratum (height = 2) and then refer to the cone of the cone to get what I want?<br>
<br>
<br>
Thanks,<br>
Max</blockquote></div><br clear="all"><div><br></div>-- <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="http://www.cse.buffalo.edu/~knepley/" target="_blank">https://www.cse.buffalo.edu/~knepley/</a><br></div></div></div></div></div></div></div></div>