<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Thu, Nov 14, 2013 at 7:38 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">On Thu, Nov 14, 2013 at 5:35 PM, Matthew Knepley <<a href="mailto:knepley@gmail.com">knepley@gmail.com</a>> wrote:<br>

> On Thu, Nov 14, 2013 at 7:32 PM, Geoffrey Irving <<a href="mailto:irving@naml.us">irving@naml.us</a>> wrote:<br>
>><br>
>> If this line in ex12 is part of Dirichlet boundary conditions, why is<br>
>> it done unconditionally even for Neumann?  Is this a mistake that<br>
>> fails to test correct Neumann conditions?<br>
>><br>
>>     ierr = DMPlexProjectFunctionLocal(dm, user.fe, user.exactFuncs,<br>
>> INSERT_BC_VALUES, userJ.u);CHKERRQ(ierr);<br>
><br>
> No, its just idempotent in that case.<br>
<br>
That's what I meant by fail to test: if there was a bug related to<br>
null spaces, it wouldn't show up in this example.<br></blockquote><div><br></div><div>I do not understand. If you have Neumann conditions, nothing happens, so</div><div>it would check that case.</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">
Thanks,<br>
Geoffrey<br>
</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>