<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Mon, Nov 2, 2015 at 7:29 PM, Barry Smith <span dir="ltr"><<a href="mailto:bsmith@mcs.anl.gov" target="_blank">bsmith@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"><br>
> On Oct 30, 2015, at 12:23 PM, Zou (Non-US), Ling <<a href="mailto:ling.zou@inl.gov">ling.zou@inl.gov</a>> wrote:<br>
><br>
> Hi All,<br>
><br>
> From physics point of view, I know my simulation converges if nothing changes any more.<br>
><br>
> I wonder how normally you do to detect if your simulation reaches steady state from numerical point of view.<br>
> Is it a good practice to use SNES convergence as a criterion, i.e.,<br>
> SNES converges and it takes 0 iteration(s)<br>
<br>
   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.<br></blockquote><div><br></div><div>There are two different situations here:</div><div><br></div><div>  1) Solving for a mathematical steady state. You remove the time derivative and solve the algebraic system with SNES. Then</div><div>      the SNES tolerance is a good measure.</div><div><br></div><div>  2) Use timestepping to advance until nothing looks like it is changing. This is a "physical" steady state.</div><div><br></div><div>You can use 1) with a timestepping preconditioner TSPSEUDO, which is what I would recommend if you</div><div>want a true steady state.</div><div><br></div><div>  Thanks,</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">
><br>
> Thanks,<br>
><br>
> Ling<br>
<br>
</blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature">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></div>