<div dir="ltr">Also, -snes_stol does not seem to work for vinewtonrsls?<br></div><div class="gmail_extra"><br><div class="gmail_quote">On Wed, Nov 30, 2016 at 8:54 PM, Justin Chang <span dir="ltr"><<a href="mailto:jychang48@gmail.com" target="_blank">jychang48@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div>Generally speaking, yes. But this is VIRS where I am enforcing maximum principles. From what I have seen, the residual typically exhibit this behavior if there are very few violations in the first place. <br><br></div>How would I "terminate on stagnation"?<br></div><div class="HOEnZb"><div class="h5"><div class="gmail_extra"><br><div class="gmail_quote">On Wed, Nov 30, 2016 at 8:48 PM, Matthew Knepley <span dir="ltr"><<a href="mailto:knepley@gmail.com" target="_blank">knepley@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div class="gmail_extra"><div class="gmail_quote"><span>On Wed, Nov 30, 2016 at 8:38 PM, Justin Chang <span dir="ltr"><<a href="mailto:jychang48@gmail.com" target="_blank">jychang48@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr">By manually terminating I meant setting -snes_max_it to 5 if I know <span class="m_-305007415447110955m_-3912563911207077406gmail-m_3282469470788906451gmail-im"><span class="m_-305007415447110955m_-3912563911207077406gmail-m_3282469470788906451gmail-m_2891339697193986227m_-9010523953303203197gmail-st">DIVERGED_LINE_SEARCH occurs after 6 iterations. In a transient simulation I cannot do this</span></span></div></blockquote><div><br></div></span><div>Let me elaborate. I mean that using DIVERGED_LINE_SEARCH as an indication of convergence is dicey. It may be that in the problem</div><div>you looked at this was true, but I see no reason to believe its true in general. Why does your residual quit decreasing?</div><div><br></div><div>  Thanks,</div><div><br></div><div>     Matt</div><span><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div class="m_-305007415447110955m_-3912563911207077406gmail-HOEnZb"><div class="m_-305007415447110955m_-3912563911207077406gmail-h5"><div class="gmail_extra"><div class="gmail_quote">On Wed, Nov 30, 2016 at 8:35 PM, Matthew Knepley <span dir="ltr"><<a href="mailto:knepley@gmail.com" target="_blank">knepley@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div class="gmail_extra"><div class="gmail_quote"><span>On Wed, Nov 30, 2016 at 8:27 PM, Justin Chang <span dir="ltr"><<a href="mailto:jychang48@gmail.com" target="_blank">jychang48@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div><div>Hi all,<br><br></div>I am running some transient simulations using <span class="m_-305007415447110955m_-3912563911207077406gmail-m_3282469470788906451m_2891339697193986227m_-9010523953303203197gmail-st">SNESVINEWTONRSLS. At certain timesteps, I get a "DIVERGED_LINE_SEARCH" which essentially "resets" my solution to zero and messes everything up. I notice that this happens when the SNES Function norm no longer decreases, and if I were to manually terminate the solver right before the final iteration I get the answer I want. <br></span></div></div></blockquote><div><br></div></span><div>Yes, its possible, however isn't that a dangerous way to terminate? Couldn't you terminate on stagnation?</div><div><br></div><div>  Thanks,</div><div><br></div><div>     Matt</div><span><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div><span class="m_-305007415447110955m_-3912563911207077406gmail-m_3282469470788906451m_2891339697193986227m_-9010523953303203197gmail-st">Is there a way to "detect" this error and use the solution from the previous non-failing iteration? Setting a fixed maximum iteration doesn't seem reasonble because every time level will require different numbers of iterations to converge.<br><br></span></div><span class="m_-305007415447110955m_-3912563911207077406gmail-m_3282469470788906451m_2891339697193986227m_-9010523953303203197gmail-st">Thanks,<br>Justin</span><span class="m_-305007415447110955m_-3912563911207077406gmail-m_3282469470788906451m_2891339697193986227m_-9010523953303203197gmail-st"></span></div>
</blockquote></span></div><span class="m_-305007415447110955m_-3912563911207077406gmail-m_3282469470788906451HOEnZb"><font color="#888888"><br><br clear="all"><div><br></div>-- <br><div class="m_-305007415447110955m_-3912563911207077406gmail-m_3282469470788906451m_2891339697193986227gmail_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>
</font></span></div></div>
</blockquote></div><br></div>
</div></div></blockquote></span></div><span><br><br clear="all"><div><br></div>-- <br><div class="m_-305007415447110955m_-3912563911207077406gmail_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>
</span></div></div>
</blockquote></div><br></div>
</div></div></blockquote></div><br></div>