On Mon, Aug 29, 2011 at 2:59 PM, Likun Tan <span dir="ltr"><<a href="mailto:likunt@andrew.cmu.edu">likunt@andrew.cmu.edu</a>></span> wrote:<br><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
Thanks,<br>
<br>
I read ex16.c, it seems that the two linear equations are solved<br>
sequentially?<br></blockquote><div><br></div><div>Yes, this is a test for reuse of preconditioner information. As I said, no one knows</div><div>how to do this for Krylov methods (and everyone has tried).</div><div><br></div>
<div> Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
On Mon, August 29, 2011 10:52 am, Satish Balay wrote:<br>
> On Mon, 29 Aug 2011, Matthew Knepley wrote:<br>
><br>
><br>
>> On Mon, Aug 29, 2011 at 2:40 PM, Likun Tan <<a href="mailto:likunt@andrew.cmu.edu">likunt@andrew.cmu.edu</a>><br>
>> wrote:<br>
>><br>
>><br>
>>> Dear all,<br>
>>><br>
>>><br>
>>> I have linear systems Ax=b with different right-hand-side vectors.<br>
>>> These<br>
>>> vectors are independent. Could i generate b simultaneously and solve<br>
>>> Ax=b<br>
>>> simultaneously?<br>
>>><br>
>><br>
>> Unfortunately, no one has figured out how to exploit this for Krylov<br>
>> methods.<br>
>><br>
><br>
> Also - check ksp/ksp/examples/tutorials/ex16.c<br>
><br>
><br>
> Satish<br>
><br>
><br>
><br>
<br>
<br>
<br>
<br>
</blockquote></div><br><br clear="all"><div><br></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>