<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Thu, May 29, 2014 at 11:31 AM, Danny Lathouwers - TNW <span dir="ltr"><<a href="mailto:D.Lathouwers@tudelft.nl" target="_blank">D.Lathouwers@tudelft.nl</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex">





<div lang="NL" link="blue" vlink="purple">
<div>
<p class="MsoNormal">Dear users,<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal"><span lang="EN-GB">I have a problem where I step time and repeatedly solve a system with differing rhs.<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-GB">At some time step petsc solver returns that initial solution is good enough (converged reason = 2 with 0 iterations done).<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-GB">I do not expect this behaviour. I use rtol = 0.001, atol=0 and dtol =large number.<u></u><u></u></span></p>
<pre><span lang="EN-GB" style="font-size:11pt;font-family:Calibri,sans-serif"><u></u> <u></u></span></pre>
<pre><span lang="EN-GB" style="font-size:11pt;font-family:Calibri,sans-serif">The manual seems to suggest the criterion is: </span><span lang="EN-GB" style="font-size:11pt;font-family:Calibri,sans-serif">rnorm < MAX (rtol * rnorm_0, abstol) (probably based on preconditioned residual).<u></u><u></u></span></pre>

<pre><span lang="EN-GB" style="font-size:11pt;font-family:Calibri,sans-serif">How could this lead to zero iterations being done? Or is the criterion based on rnorm/bnorm instead (which I found in some reference on the internet concerning petsc and would explain the observed behaviour).</span></pre>
</div></div></blockquote><div>Its ||b||, not ||r_0||. You can change it to get the other behavior, as detailed here</div><div><br></div><div>  <a href="http://www.mcs.anl.gov/petsc/petsc-dev/src/ksp/ksp/interface/iterativ.c.html#KSPConvergedDefault">http://www.mcs.anl.gov/petsc/petsc-dev/src/ksp/ksp/interface/iterativ.c.html#KSPConvergedDefault</a></div>
<div><br></div><div>     Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex"><div lang="NL" link="blue" vlink="purple">
<div>
<pre><span lang="EN-GB" style="font-size:11pt;font-family:Calibri,sans-serif">Thanks,<u></u><u></u></span></pre>
<pre><span lang="EN-GB" style="font-size:11pt;font-family:Calibri,sans-serif">Danny.<u></u><u></u></span></pre>
<p class="MsoNormal"><span lang="EN-GB" style="color:rgb(31,73,125)"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-GB"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-GB"><u></u> <u></u></span></p>
</div>
</div>

</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
</div></div>