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?

Mark

Picard iteration 0 max(resid) = 781020
  0 KSP Residual norm 2.499264167662e+07 
  1 KSP Residual norm 6.330915810095e+06 
  2 KSP Residual norm 1.856999435466e+06 
  3 KSP Residual norm 3.967325725053e+05 
  4 KSP Residual norm 9.037332775808e+04 
  5 KSP Residual norm 1.826960320400e+04 
Linear solve converged due to CONVERGED_RTOL iterations 5
Picard iteration 1 max(resid) = 727577 ------- Rate = 1.07345
  0 KSP Residual norm 6.885328126772e+06 
  1 KSP Residual norm 2.150319244570e+06 
  2 KSP Residual norm 2.297759417371e+05 
  3 KSP Residual norm 3.321819214451e+04 
  4 KSP Residual norm 4.155970071625e+03 
Linear solve converged due to CONVERGED_RTOL iterations 4
Picard iteration 2 max(resid) = 486788 ------- Rate = 1.49465
  0 KSP Residual norm 2.987417661953e+06 
  1 KSP Residual norm 5.046119985008e+05 
  2 KSP Residual norm 3.882757664605e+04 
  3 KSP Residual norm 4.106091678345e+03 
Linear solve converged due to CONVERGED_RTOL iterations 3
Picard iteration 3 max(resid) = 137348 ------- Rate = 3.54418
  0 KSP Residual norm 7.285359985607e+05 
  1 KSP Residual norm 9.132700545850e+04 
  2 KSP Residual norm 6.834608381922e+03 
Linear solve converged due to CONVERGED_RTOL iterations 2
Picard iteration 4 max(resid) = 25082.7 ------- Rate = 5.47582
  0 KSP Residual norm 1.207523847817e+05 
  1 KSP Residual norm 1.286174706329e+04 
Linear solve converged due to CONVERGED_RTOL iterations 1
Picard iteration 5 max(resid) = 5723.21 ------- Rate = 4.38264
  0 KSP Residual norm 2.211945822133e+04 
Linear solve converged due to CONVERGED_RTOL iterations 0


Reply via email to