<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Tue, Nov 19, 2013 at 5:27 PM, Geoffrey Irving <span dir="ltr"><<a href="mailto:irving@naml.us" target="_blank">irving@naml.us</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">This question isn't about changing ex12, I'm just making sure I<br>
understand what's going on:<br>
<br>
SetupSection in snes ex12 calls DMPlexLabelComplete on the "marker"<br>
label if Dirichlet conditions are chosen (not if Neumann is chosen).<br>
Specifically,<br>
<br>
    if (user->bcType == DIRICHLET) {ierr = DMPlexLabelComplete(dm,<br>
label);CHKERRQ(ierr);}<br>
<br>
It then uses the same label ("marker") to create an index set:<br>
<br>
    if (user->bcType == DIRICHLET) {ierr  = DMPlexGetStratumIS(dm,<br>
bdLabel, 1, &bcPoints[0]);CHKERRQ(ierr);}<br>
<br>
I believe 1 is the codimension.  If the index set is restricted to<br>
codimension 1, and DMPlexLabelComplete labels vertices based on edge<br>
labels, isn't the DMPlexLabelComplete call unnecessary?<br></blockquote><div><br></div><div>1) The value 1 is just a convention that I use for marking the boundary. In a real simulation,</div><div>     there are probably many makers and they are user specified. This is what happens in</div>
<div>     TS ex11.</div><div><br></div><div>2) I should probably complete the label in the Neumann case, but I just forgot. The completion</div><div>     is not necessary the way my defaults mesh generators work, but is sometimes necessary</div>
<div>     when reading stuff in.</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">
Hopefully either all that made sense and the answer is yes, or someone<br>
can point out the hole in my understanding.<br>
<br>
Thanks for the help with all of my recent emails!  I am getting closer<br>
to the actual solver part of ex12, which I understand from previous<br>
use, so hopefully I will have a working rewritten example in not too<br>
long.<br>
<span class="HOEnZb"><font color="#888888"><br>
Geoffrey<br>
</font></span></blockquote></div><br><br clear="all"><div><br></div>-- <br>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>