<div class="gmail_quote">On Sun, Jun 3, 2012 at 8:37 AM, Niceno Bojan <span dir="ltr"><<a href="mailto:bojan.niceno@psi.ch" target="_blank">bojan.niceno@psi.ch</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div id=":5lw">I don't know.  Is there a way to find out?<br></div></blockquote><div><br></div><div>-ksp_monitor or -ksp_monitor_true_residual</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div id=":5lw">
<div class="im"><br>
<br>
> In that case, check that the system being solved is identical. You can save<br>
> it with -ksp_view_binary. Load it with<br>
> src/ksp/ksp/examples/tutorials/ex10.c as an independent test. Also be sure<br>
> to use the same initial guess.<br>
<br>
</div>Thanks.  I was checking norms of matrix "A" and "x" and "b" vectors (2nd order) and they ARE the same before and after restart.  Does the same norm ensure the matrices and and vectors are the same?</div>
</blockquote></div><br><div>They could be ordered differently. A different ordering would change the ILU algorithm. You can also run with -pc_type jacobi to avoid this source of variation (norms are still order-dependent up to rounding error).</div>