<div dir="ltr">I cannot surely say my Jacobian for this particular problem is correct, as I have not checked it. For a smaller problem, I've already checked its correctness using -snes_type test or -snes_compare_explicit (but linear solve and nonlinear solve with FD one need a few more iterations than with my Jacobian). To make it sure, now I started -snes_type test for the problem and will update you once it finished. By the way, I'm waiting the calculation for more than three hours now. Is it usual for a large problem (>1e6 dof) or is there something wrong? <br>
<div><br></div><br><div><br></div></div><div class="gmail_extra"><br><br><div class="gmail_quote">On Mon, Apr 28, 2014 at 6:34 AM, 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>
I have run your code. I changed to use -snes_type newtonls and also -snes_mf_operator there is something wrong with your Jacobian:<br>
<br>
Without -snes_mf_operator<br>
0 SNES Function norm 1.821611413735e+03<br>
0 KSP Residual norm 1821.61<br>
<div class=""> Linear solve converged due to CONVERGED_ITS iterations 1<br>
</div> Linear solve converged due to CONVERGED_RTOL iterations 1<br>
1 KSP Residual norm 0.000167024<br>
<div class=""> Linear solve converged due to CONVERGED_ITS iterations 1<br>
</div> Linear solve converged due to CONVERGED_RTOL iterations 1<br>
2 KSP Residual norm 7.66595e-06<br>
<div class=""> Linear solve converged due to CONVERGED_ITS iterations 1<br>
</div> Linear solve converged due to CONVERGED_RTOL iterations 1<br>
3 KSP Residual norm 4.4581e-07<br>
<div class=""> Linear solve converged due to CONVERGED_ITS iterations 1<br>
</div> Linear solve converged due to CONVERGED_RTOL iterations 1<br>
4 KSP Residual norm 3.77537e-08<br>
<div class=""> Linear solve converged due to CONVERGED_ITS iterations 1<br>
</div> Linear solve converged due to CONVERGED_RTOL iterations 1<br>
5 KSP Residual norm 2.20453e-09<br>
<div class=""> Linear solve converged due to CONVERGED_ITS iterations 1<br>
</div> Linear solve converged due to CONVERGED_RTOL iterations 1<br>
6 KSP Residual norm 1.711e-10<br>
Linear solve converged due to CONVERGED_RTOL iterations 6<br>
<br>
with -snes_mf_operator<br>
<br>
0 SNES Function norm 1.821611413735e+03<br>
0 KSP Residual norm 1821.61<br>
<div class=""> Linear solve converged due to CONVERGED_ITS iterations 1<br>
</div> Linear solve converged due to CONVERGED_RTOL iterations 1<br>
1 KSP Residual norm 1796.39<br>
<div class=""> Linear solve converged due to CONVERGED_ITS iterations 1<br>
</div> Linear solve converged due to CONVERGED_RTOL iterations 1<br>
2 KSP Residual norm 1786.2<br>
<div class=""> Linear solve converged due to CONVERGED_ITS iterations 1<br>
</div> Linear solve converged due to CONVERGED_RTOL iterations 1<br>
3 KSP Residual norm 1741.11<br>
<div class=""> Linear solve converged due to CONVERGED_ITS iterations 1<br>
</div> Linear solve converged due to CONVERGED_RTOL iterations 1<br>
4 KSP Residual norm 1733.92<br>
<div class=""> Linear solve converged due to CONVERGED_ITS iterations 1<br>
</div> Linear solve converged due to CONVERGED_RTOL iterations 1<br>
5 KSP Residual norm 1726.57<br>
<div class=""> Linear solve converged due to CONVERGED_ITS iterations 1<br>
</div> Linear solve converged due to CONVERGED_RTOL iterations 1<br>
6 KSP Residual norm 1725.35<br>
<div class=""> Linear solve converged due to CONVERGED_ITS iterations 1<br>
</div> Linear solve converged due to CONVERGED_RTOL iterations 1<br>
7 KSP Residual norm 1723.89<br>
<div class=""> Linear solve converged due to CONVERGED_ITS iterations 1<br>
</div> Linear solve converged due to CONVERGED_RTOL iterations 1<br>
8 KSP Residual norm 1715.41<br>
<div class=""> Linear solve converged due to CONVERGED_ITS iterations 1<br>
</div> Linear solve converged due to CONVERGED_RTOL iterations 1<br>
9 KSP Residual norm 1713.72<br>
<div class=""> Linear solve converged due to CONVERGED_ITS iterations 1<br>
</div> Linear solve converged due to CONVERGED_RTOL iterations 1<br>
10 KSP Residual norm 1702.84<br>
<div class=""> Linear solve converged due to CONVERGED_ITS iterations 1<br>
<br>
</div>…<br>
<br>
This means your Jacobian is wrong. Your first order of business is to fix your Jacobian. I noticed in previous emails your discussion with Jed about switching to MatGetLocalSubMatrix() and using -snes_type test YOU NEED TO DO THIS. You will get no where with an incorrect Jacobian. You need to fix your Jacobian before you do anything else! No amount of other options or methods will help you with a wrong Jacobian! Once you have a correct Jacobian if you still have convergence problems let us know and we can make further suggestions.<br>
<span class="HOEnZb"><font color="#888888"><br>
Barry<br>
</font></span><div class="im HOEnZb"><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>
</div><div class="HOEnZb"><div class="h5">> Hi,<br>
><br>
> In my simulation, nonlinear solve with the trust regtion method got stagnent after linear solve (see output below). Is it possible that the method goes to inifite loop? Is there any parameter to avoid this situation?<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>