<div dir="ltr"><div><div><div>Thanks,<br><br></div>I got this message  <i>Linear solve converged due to CONVERGED_RTOL iterations 2</i>. <br><br></div>What is the definition of RTOL in petsc?<br><br></div>I doubt the RTOL is small.<i><br>
<br><br></i></div><div class="gmail_extra"><br><br><div class="gmail_quote">On Wed, Oct 30, 2013 at 6:38 PM, Jed Brown <span dir="ltr"><<a href="mailto:jedbrown@mcs.anl.gov" target="_blank">jedbrown@mcs.anl.gov</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div class="im">Fande Kong <<a href="mailto:fd.kong@siat.ac.cn">fd.kong@siat.ac.cn</a>> writes:<br>
<br>
> Hi all,<br>
><br>
> How to force the code run with a given number of iterations?<br>
><br>
> I use the following script to run code, but the code stop after doing<br>
</div>> *two*iterations:<br>
><br>
> * mpirun -n 2 ./linearElasticity  -ksp_type gmres -pc_type asm<br>
<div class="im">> -pc_asm_overlap 1 -sub_pc_type ilu -sub_pc_factor_levels 2 -sub_ksp_type<br>
> preonly -sub_pc_factor_mat_ordering_type rcm -ksp_view - -ksp_rtol 1e-5<br>
</div>> -ksp_max_it 500 -ksp_monitor -ksp_monitor_singular_value*<br>
<div class="im">><br>
> The residual history is listed as follows:<br>
><br>
</div>> *  0 KSP Residual norm 9.732477596672e+31<br>
<div class="im">>   0 KSP Residual norm 9.732477596672e+31 % max 1.000000000000e+00 min<br>
> 1.000000000000e+00 max/min 1.000000000000e+00<br>
>   1 KSP Residual norm 8.532845325800e+31<br>
>   1 KSP Residual norm 8.532845325800e+31 % max 4.821161757157e+55 min<br>
> 4.821161757157e+55 max/min 1.000000000000e+00<br>
>   2 KSP Residual norm 3.648279603061e+23<br>
>   2 KSP Residual norm 3.648279603061e+23 % max 1.234581012401e+79 min<br>
> 2.318814661396e+55 max/min 5.324190125905e+23<br>
> *<br>
>  The residual is not small enough to stop.<br>
<br>
</div>-ksp_converged_reason will tell you why it exited early.<br>
</blockquote></div><br></div>