<div dir="ltr">I mean "it's <span class="GINGER_SOFTWARE_mark" id="514aeb7c-1374-4293-9b72-aaf4d913d8c4">keep</span> <span class="GINGER_SOFTWARE_mark" id="315b78c7-5731-417e-a387-a446fe8c0fab">running but</span> not print anything". I guess the program is running in while<span class="GINGER_SOFTWARE_mark" id="015fe25b-a8df-45e4-b7b4-062280a2c541">(</span>1) loop after linear solve in SNESSolve_NEWTONTR<span class="GINGER_SOFTWARE_mark" id="dde5a4c4-0430-4f37-befc-5d91cf451416">(</span>) in <span class="GINGER_SOFTWARE_mark" id="d932ef5e-9585-40b7-8a03-44e169f20317">src</span>/<span class="GINGER_SOFTWARE_mark" id="bafca501-abe5-4aaa-8e78-d22d2423a838">snes</span>/<span class="GINGER_SOFTWARE_mark" id="2c51629b-bf08-4746-992f-9046be55e819">impls</span>/<span class="GINGER_SOFTWARE_mark" id="1df2f64e-cae9-4e6f-bb24-8cf5352e40fd">tr</span>/<span class="GINGER_SOFTWARE_mark" id="ff00d495-0a41-4f7d-ac44-65358d69ba8f">tr.c</span>. But I'm not sure if the breaking condition (if (rho > <span class="GINGER_SOFTWARE_mark" id="cd17af90-a650-4b2f-af99-79ffd0c6494c">neP</span>->sigma) break;<span class="GINGER_SOFTWARE_mark" id="1f716fea-70da-4e83-8367-29f93ee1e458"> )</span>  can be always satisfied at the end.<br>
<pre width="80"><strong><font color="#4169E1"><br></font></strong></pre></div><div class="gmail_extra"><br><br><div class="gmail_quote">On Fri, Apr 25, 2014 at 2:49 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"><div class=""><br>
On Apr 25, 2014, at 7:31 AM, Norihiro Watanabe <<a href="mailto:norihiro.w@gmail.com">norihiro.w@gmail.com</a>> wrote:<br>
<br>
> Hi,<br>
><br>
> In my simulation, nonlinear solve with the trust regtion method got stagnent after linear solve (see output below).<br>
<br>
</div>   What do you mean, get stagnant? Does the code just hang, that is keep running but not print anything.<br>
<div class=""><br>
> Is it possible that the method goes to inifite loop?<br>
<br>
</div>   This is not suppose to be possible.<br>
<div class=""><br>
> Is there any parameter to avoid this situation?<br>
<br>
</div>  You need to determine what it is “hanging” on. Try running with -start_in_debugger and when it “hangs” hit control C and type where to determine where it is.<br>
<span class="HOEnZb"><font color="#888888"><br>
   Barry<br>
</font></span><div class="HOEnZb"><div class="h5"><br>
><br>
>   0 SNES Function norm 1.828728087153e+03<br>
>     0 KSP Residual norm 91.2735<br>
>     Linear solve converged due to CONVERGED_ITS iterations 1<br>
>     Linear solve converged due to CONVERGED_RTOL iterations 3<br>
>     1 KSP Residual norm 3.42223<br>
>   Linear solve converged due to CONVERGED_STEP_LENGTH iterations 1<br>
><br>
><br>
> Thank you in advance,<br>
> Nori<br>
<br>
</div></div></blockquote></div><br><br clear="all"><br>-- <br>Norihiro Watanabe
</div>