<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Mon, Nov 10, 2014 at 10:28 PM, Jed Brown <span dir="ltr"><<a href="mailto:jed@jedbrown.org" target="_blank">jed@jedbrown.org</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><span class="">Matthew Knepley <<a href="mailto:knepley@gmail.com">knepley@gmail.com</a>> writes:<br>
>> at which point the nullspace is not attached, no?<br>
>><br>
><br>
> Crap, that is right. I think we should propagate the nullspace from pmat[1].<br>
<br>
</span>Uh, there is no reason pmat[1] can be expected to have the same null<br>
space.  pmat[1] is the zero matrix in some important cases, and the null<br>
space of the Schur complement (or its SIMPLE-type approximation) depends<br>
on boundary conditions in a different block (thus problematic for<br>
segregated assembly).  I've been unsuccessful so far in my search for an<br>
elegant solution, but passing the null space along in this way is<br>
encouraging overt lying and can be expected to break other algorithms as<br>
well as any consistency tests we might add in the future.<br>
</blockquote></div><br>I know it can break, however here is my thinking. The user is completely in control</div><div class="gmail_extra">here, and this is the only channel we have to pass the information. So I would make</div><div class="gmail_extra">them responsible for breakage. We have excellent diagnostics here, so we can tell</div><div class="gmail_extra">who has a null space right away. We have a lot of interface that allow a user to break</div><div class="gmail_extra">PETSc in order to allow them to customize it sufficiently.</div><div class="gmail_extra"><br></div><div class="gmail_extra">  Thanks,</div><div class="gmail_extra"><br></div><div class="gmail_extra">     Matt<br clear="all"><div><br></div>-- <br><div class="gmail_signature">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></div>