<div dir="ltr">On Mon, Jan 7, 2013 at 8:12 PM, Barry Smith <span dir="ltr"><<a href="mailto:bsmith@mcs.anl.gov" target="_blank">bsmith@mcs.anl.gov</a>></span> wrote:<br><div class="gmail_extra"><div class="gmail_quote">
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div class="im"><br>
On Jan 7, 2013, at 8:05 PM, Jed Brown <<a href="mailto:jedbrown@mcs.anl.gov">jedbrown@mcs.anl.gov</a>> wrote:<br>
<br>
> It's sort of awkward to deal with PCREDUNDANT for coarse grids, and I have yet to see any evidence that it's better than gathering to rank 0, solving sequentially, and scattering the result. What about implementing LU for MPIAIJ to do the latter algorithm? Then we could trivially make it smart enough to skip the gather when Mark uses it on the coarse grid of GAMG (where the entire system is already on rank 0). In additon to being simpler for the user to interact with, this would fix a lot of fighting between the PCMG and its children.<br>

<br>
</div>   I thought PCREDUNDANT had this ability? What about adding it?</blockquote><div><br></div><div style>That wouldn't be redundant at all, so if we put it there, we'd have to at least rename it.</div><div> </div>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"> Have an option that cause PCREDUNDANT to only collect a single copy of the matrix from rank 0 to r-1 and in your case r is 1.  I hate the idea of having another separate place with similar stuff to PCREDUNDANT.  You could then implement LU for MPIAIJ using this variant of PCREDUNDANT, that is a call to PCREDUNDANT stuff with the right options.<br>

<span class="HOEnZb"><font color="#888888"><br>
   Barry<br>
<br>
<br>
<br>
</font></span></blockquote></div><br></div></div>