<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Mon, Jan 27, 2014 at 4:27 PM, Danyang Su <span dir="ltr"><<a href="mailto:danyang.su@gmail.com" target="_blank">danyang.su@gmail.com</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">Hi All,<br>
<br>
I implement KSP solver in a Newton iteration. The convergence parameters are: rtol 1.0E-8, abstol 1.0E-50, dtol 1.0E5, maxits 1000.<br></blockquote><div><br></div><div>These tolerances use the preconditioned residual norm. You can sue the true residual norm for some solver types:</div>
<div><br></div><div>  <a href="http://www.mcs.anl.gov/petsc/petsc-dev/docs/manualpages/KSP/KSPSetNormType.html">http://www.mcs.anl.gov/petsc/petsc-dev/docs/manualpages/KSP/KSPSetNormType.html</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">
Does the following convergence monitoring make sense? I thought the solver should stop before iteration 36, but it stops at iteration 70.<br>
<br>
  0 KSP preconditioned resid norm 4.208619583781e-002 true resid norm 2.227983625699e-005 ||r(i)||/||b|| 1.000000000000e+000<br>
  1 KSP preconditioned resid norm 6.403275610900e-003 true resid norm 2.795175726580e-007 ||r(i)||/||b|| 1.254576422528e-002<br>
  2 KSP preconditioned resid norm 7.416300263102e-004 true resid norm 3.382118379502e-010 ||r(i)||/||b|| 1.518017610404e-005<br>
  3 KSP preconditioned resid norm 6.144409100696e-005 true resid norm 1.260466100377e-009 ||r(i)||/||b|| 5.657429820568e-005<br>
  4 KSP preconditioned resid norm 4.990158094666e-005 true resid norm 3.632496726225e-008 ||r(i)||/||b|| 1.630396509349e-003<br>
  5 KSP preconditioned resid norm 2.743769015994e-005 true resid norm 1.101666758939e-009 ||r(i)||/||b|| 4.944680680017e-005<br>
  6 KSP preconditioned resid norm 4.269805798759e-006 true resid norm 7.915494521121e-011 ||r(i)||/||b|| 3.552761532813e-006<br>
  7 KSP preconditioned resid norm 1.632413041808e-006 true resid norm 1.163794439302e-009 ||r(i)||/||b|| 5.223532282187e-005<br>
  8 KSP preconditioned resid norm 8.010590587547e-007 true resid norm 1.246839309256e-011 ||r(i)||/||b|| 5.596267831030e-007<br>
  9 KSP preconditioned resid norm 4.794695759395e-007 true resid norm 3.155943984006e-012 ||r(i)||/||b|| 1.416502324166e-007<br>
 10 KSP preconditioned resid norm 4.500033799628e-007 true resid norm 4.554787575528e-011 ||r(i)||/||b|| 2.044354151884e-006<br>
...<br>
 36 KSP preconditioned resid norm 2.752424650444e-009 true resid norm 2.990835953592e-016 ||r(i)||/||b|| 1.342395841286e-011<br>
...<br>
 51 KSP preconditioned resid norm 5.761764412500e-010 true resid norm 2.557370618510e-016 ||r(i)||/||b|| 1.147840849911e-011<br>
...<br>
 69 KSP preconditioned resid norm 4.324655127362e-010 true resid norm 1.982933485715e-014 ||r(i)||/||b|| 8.900125938281e-010<br>
 70 KSP preconditioned resid norm 3.923752626311e-010 true resid norm 1.065002491928e-015 ||r(i)||/||b|| 4.780118128535e-011<br>
<br>
Thanks,<br>
<br>
Danyang<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
</div></div>