<div dir="ltr"><br><div class="gmail_extra"><br><br><div class="gmail_quote">On Thu, Sep 19, 2013 at 2:53 PM, Barry Smith <span dir="ltr"><<a href="mailto:bsmith@mcs.anl.gov" target="_blank">bsmith@mcs.anl.gov</a>></span> wrote:<br>

<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"><br>
   It is essentially 3 Laplacians so I think the default null space of 3 constant vectors is fine. The problem is without the block information presumably GAMG is using only a single constant vector over all variables?  So maybe we need to construct a 3 vector null space which just marks in the reduced vector from which of the 3 components each entry came from.<br>

</blockquote><div> </div><div> SNESVIGetActiveSetIS() currently doesn't try to detect or set the inactive set IS block structure, so MatGetSubMatrix() can't take advantage of that.  This (blocking) information is geometric in nature -- since it is mesh nodes that really come "into contact" -- so unless we involve a DM somehow, I don't know if it is possible to retain the right blocking structure in general.  Or add a user callback hook to compute the inactive set?</div>

<div><br></div><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">
<span class=""><font color="#888888"><br>
   Barry<br>
</font></span><div class=""><div class="h5"><br>
On Sep 19, 2013, at 3:48 PM, Jed Brown <<a href="mailto:jedbrown@mcs.anl.gov">jedbrown@mcs.anl.gov</a>> wrote:<br>
<br>
> Barry Smith <<a href="mailto:bsmith@mcs.anl.gov">bsmith@mcs.anl.gov</a>> writes:<br>
>>   That is what it is doing and apparently it doesn't result in a good<br>
>>   preconditioner; I don't know why off hand. One thing is it no<br>
>>   longer knows about the block structure.<br>
><br>
> How is the near null space being specified?<br>
<br>
</div></div></blockquote></div><br></div><br clear="all"><div><br></div><br></div>