<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Wed, Aug 31, 2016 at 10:32 AM, Ivano Barletta <span dir="ltr"><<a href="mailto:ibarletta@inogs.it" target="_blank">ibarletta@inogs.it</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div><div><div><div><div><div><div><div><div><div><div><div><div><div><div><div>Dear Petsc Users<br><br></div>I'm using Petsc to solve an elliptic equation<br><br></div>The code can be run in parallel but I'm running<br></div>some tests in sequential by the moment<br><br></div>When I look at the output, what it looks odd to <br></div>me is that the relative tolerance that I set is <br></div>not fulfilled. <br></div>I've set -ksp_rtol 1e-8 in my runtime options<br></div>but the solver stops when the ratio<br></div>|| r || / || b || is still 9e-8, then almost<br></div>one order of magnitude greater of the rtol <br></div>that I set (as you can see in the txt in attachment).<br><br></div>My question is, isn't the solver supposed to <br></div>make other few iterations to reach the relative tolerance?<br></div></div></div></div></blockquote><div><br></div><div>Here is the code:</div><div><br></div><div>  <a href="https://bitbucket.org/petsc/petsc/src/a4f377b2def412a6138fd768d05e3a45fdbb68fb/src/ksp/ksp/interface/iterativ.c?at=master&fileviewer=file-view-default#iterativ.c-725">https://bitbucket.org/petsc/petsc/src/a4f377b2def412a6138fd768d05e3a45fdbb68fb/src/ksp/ksp/interface/iterativ.c?at=master&fileviewer=file-view-default#iterativ.c-725</a></div><div><br></div><div>If you have a zero initial guess, we do not use ||b||, but ||r0||. In this case,</div><div>the ratio</div><div><br></div><div>  ||r_75|| / ||r_0|| < 1.0e-8</div><div><br></div><div>when using preconditioned residuals.</div><div><br></div><div>  Thanks,</div><div><br></div><div>     Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div><div><div></div>Thanks in advance for replies and suggestions<br></div><div>Kind Regards <br></div><span class=""><font color="#888888">Ivano <br></font></span><br></div>P.S. my runtime options are these:<br>-ksp_monitor_true_residual -ksp_type cg -ksp_converged_reason -ksp_view -ksp_rtol 1e-8 <br><div><div><div><div><div><div><div><div><div><div><div><br></div></div></div></div></div></div></div></div></div></div></div></div>
</blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature" data-smartmail="gmail_signature">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></div>