<div class="gmail_quote">On Sun, Jun 3, 2012 at 8:30 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=":2aw">Residual norm obtained from a call to "KSPGetResidualNorm(ksp, &rnorm);", gives:<br>
<br>
- First run: 1.98077e-02<br>
- Second run: 2.61329e-02<br>
<br>
The tolerances I request from KSP's gmres solver are:<br>
<br>
- Linear solver relative tolerance = 0.001;<br>
- Linear solver absolute tolerance = 1.0e-5;<br>
<br>
Does it tell you something?  It does not to me :-(<br></div></blockquote><div><br></div><div>Is the initial residual about 30?</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div id=":2aw"><div class="im">
> For some codes, the domain decomposition will not be identical on restart;<br>
> that changes the algorithm.<br>
<br>
</div>That is clear, but what I am comparing are two sequential runs.</div></blockquote></div><br><div>In that case, check that the system being solved is identical. You can save it with -ksp_view_binary. Load it with src/ksp/ksp/examples/tutorials/ex10.c as an independent test. Also be sure to use the same initial guess.</div>