On Wed, Dec 14, 2011 at 12:45 PM, Mark F. Adams <span dir="ltr"><<a href="mailto:mark.adams@columbia.edu">mark.adams@columbia.edu</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">
I have a nonlinear problem and am finding the the convergence test for the interior PETSc linear solver is apparently testing against the norm of the first linear solve (see appended).  This is clearly not what anyone would want, is this the intended semantics?  Is there perhaps a method to reset rtol?<br>
</blockquote><div><br></div><div>-snes_view -ksp_view?</div><div><br></div><div>  Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Mark<br>
<br>
Picard iteration 0 max(resid) = 781020<br>
  0 KSP Residual norm 2.499264167662e+07<br>
  1 KSP Residual norm 6.330915810095e+06<br>
  2 KSP Residual norm 1.856999435466e+06<br>
  3 KSP Residual norm 3.967325725053e+05<br>
  4 KSP Residual norm 9.037332775808e+04<br>
  5 KSP Residual norm 1.826960320400e+04<br>
Linear solve converged due to CONVERGED_RTOL iterations 5<br>
Picard iteration 1 max(resid) = 727577 ------- Rate = 1.07345<br>
  0 KSP Residual norm 6.885328126772e+06<br>
  1 KSP Residual norm 2.150319244570e+06<br>
  2 KSP Residual norm 2.297759417371e+05<br>
  3 KSP Residual norm 3.321819214451e+04<br>
  4 KSP Residual norm 4.155970071625e+03<br>
Linear solve converged due to CONVERGED_RTOL iterations 4<br>
Picard iteration 2 max(resid) = 486788 ------- Rate = 1.49465<br>
  0 KSP Residual norm 2.987417661953e+06<br>
  1 KSP Residual norm 5.046119985008e+05<br>
  2 KSP Residual norm 3.882757664605e+04<br>
  3 KSP Residual norm 4.106091678345e+03<br>
Linear solve converged due to CONVERGED_RTOL iterations 3<br>
Picard iteration 3 max(resid) = 137348 ------- Rate = 3.54418<br>
  0 KSP Residual norm 7.285359985607e+05<br>
  1 KSP Residual norm 9.132700545850e+04<br>
  2 KSP Residual norm 6.834608381922e+03<br>
Linear solve converged due to CONVERGED_RTOL iterations 2<br>
Picard iteration 4 max(resid) = 25082.7 ------- Rate = 5.47582<br>
  0 KSP Residual norm 1.207523847817e+05<br>
  1 KSP Residual norm 1.286174706329e+04<br>
Linear solve converged due to CONVERGED_RTOL iterations 1<br>
Picard iteration 5 max(resid) = 5723.21 ------- Rate = 4.38264<br>
  0 KSP Residual norm 2.211945822133e+04<br>
Linear solve converged due to CONVERGED_RTOL iterations 0<br>
<br>
</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>