On Wed, May 30, 2012 at 3:16 PM, John Fettig <span dir="ltr"><<a href="mailto:john.fettig@gmail.com" target="_blank">john.fettig@gmail.com</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">
On Wed, May 30, 2012 at 3:04 PM, Jed Brown <<a href="mailto:jedbrown@mcs.anl.gov">jedbrown@mcs.anl.gov</a>> wrote:<br>
> On Wed, May 30, 2012 at 1:58 PM, John Fettig <<a href="mailto:john.fettig@gmail.com">john.fettig@gmail.com</a>> wrote:<br>
>><br>
>> I have a question about the difference between SAME_PRECONDITIONER and<br>
>> DIFFERENT_NONZERO_PATTERN.  If you solve multiple RHS's (with<br>
>> different initial guesses) but the matrix does not change, should you<br>
>> expect to get identical results with SAME_PRECONDITIONER and<br>
>> DIFFERENT_NONZERO_PATTERN?<br>
><br>
><br>
> Why are you even calling KSPSetOperators? Just call KSPSolve() for each<br>
> right hand side.<br>
<br>
What I mean is if I run the code twice, once with SAME_PRECONDITIONER<br>
and once with DIFFERENT_NONZERO_PATTERN, the solution for the second<br>
solve in the code is different.  I am sure the matrix does not change.<br>
 There are a couple of calls to MatZeroRows but I have verified that<br>
the diagonal is set to the exact same value both time.  The difference<br>
in the solution is not insignificant (it is probably below the solver<br>
tolerance but still in the 3rd significant digit).</blockquote><div><br></div><div>Is the initial residual the same. It would be helpful to see the -ksp_monitor_true_residual</div><div>output</div><div><br></div><div>   Matt</div>
<div><br></div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><span class="HOEnZb"><font color="#888888"><br>
John<br>
</font></span></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>