[petsc-users] How do I know it is steady state?

Barry Smith bsmith at mcs.anl.gov
Mon Nov 2 19:29:23 CST 2015


> On Oct 30, 2015, at 12:23 PM, Zou (Non-US), Ling <ling.zou at inl.gov> wrote:
> 
> Hi All,
> 
> From physics point of view, I know my simulation converges if nothing changes any more.
> 
> I wonder how normally you do to detect if your simulation reaches steady state from numerical point of view.
> Is it a good practice to use SNES convergence as a criterion, i.e.,
> SNES converges and it takes 0 iteration(s)

   Depends on the time integrator and SNES tolerance you are using. If you use a -snes_rtol 1.e-5 it will always try to squeeze 5 MORE digits out of the residual so won't take 0 iterations even if there is only a small change in the solution. 
> 
> Thanks,
> 
> Ling



More information about the petsc-users mailing list