<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Wed, Nov 20, 2013 at 10:38 AM, 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"><div class="HOEnZb"><div class="h5">On Tue, Nov 19, 2013 at 10:32 PM, Jed Brown <<a href="mailto:jedbrown@mcs.anl.gov">jedbrown@mcs.anl.gov</a>> wrote:<br>

> Geoffrey Irving <<a href="mailto:irving@naml.us">irving@naml.us</a>> writes:<br>
>> I was missing the obvious INSERT_BC_VALUES parameter in the call to<br>
>> DMPlexProjectFunctionLocal.  It's a shame one has to loop over the<br>
>> interior just to set the boundary, but it does seem like the code<br>
>> produces the correct answer.<br>
><br>
> You can of course build an index to loop over any chosen feature fast.<br>
> But with no index, an iterator for a feature would still involve<br>
> visiting everything.<br>
<br>
</div></div>That was part of why I was confused: in snes ex12, we *do* build a<br>
boundary index set (SetupSection).  I haven't traced the code to see<br>
if the IS itself is discarded, but the DM at least has access to it in<br>
some form as the boundary section of its default section.</blockquote><div><br></div><div>The problem here is generality. The INSERT_BC_VALUES sets anything that is</div><div>constrained in the PetscSection. The constraints can come from many labels</div>
<div>and many values, so that no one of those ISes is an index. Adding support for an</div><div>index over constrained points in a PetscSection is easy (basically the same as the</div><div>DMLabel code).</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"><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>