<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex"><div class=""><div class="h5"><br>
<br>
</div></div>We could make what you have (finding rows with just diagonals) the<br>
default.  For, e.g., vector-valued problems with unaligned partial<br>
Dirichlet conditions, the user could specify something else.<br></blockquote><div><br></div><div>Sure.  But unaligned BCs are fine.  This problem just comes from users doing something stupid from our perspective (and not innately stupid) and putting a "big" fake number on the diagonal for the BC.</div><div><br></div><div>Note, unaligned BCs require the null space be adjusted accordingly so my syntactic sugar (that Jed hates) of taking coordinates and doing this for users is not quite correct.  Bit probably does not matter at BCs in practice anyway.</div><div><br></div></div></div></div>