On Tue, May 26, 2009 at 12:29 PM, Jed Brown <span dir="ltr"><<a href="mailto:jed@59a2.org">jed@59a2.org</a>></span> wrote:<br><div class="gmail_quote"><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
<div class="im">Matthew Knepley wrote:<br>
<br>
> It not computing Ab. Its computing Ab, A^2b, A^3b, etc.<br>
<br>
</div>Of course, but that's all possible if A*b is computable.  Lisandro<br>
demonstrated that A*x was fine, which it better be because of how x was<br>
obtained.  Had he tried A*b, which is the first thing the KSP needs, the<br>
problem would have been obvious.</blockquote><div><br>This is the whole point. b and x cannot have different local sizes in a Krylov<br>method. It makes no sense.<br><br>  Matt<br> </div><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
<font color="#888888"><br>
Jed<br></font></blockquote></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<br>