<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Sun, Jul 19, 2015 at 9:45 AM, Khai Pham <span dir="ltr"><<a href="mailto:khaipham@utexas.edu" target="_blank">khaipham@utexas.edu</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><div><div><div>Dear Petsc-Users,<br><br></div>The SNES only converges for a small problem. I did the jacobian check. The result is sensitive to the
choice of differencing parameter.<br></div></div></div></div></blockquote><div><br></div><div>This looks like you probably have the correct Jacobian, but the problem is very sensitive. Then convergence might be</div><div>hard as you scale up because you have a bad initial guess. Its impossible to say anything else without knowing more</div><div>about your problem.</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"><div dir="ltr"><div><div><div>With the option "-snes_type test":<br><br>"Norm of matrix ratio 4.11569e-09 difference 7.64426e-07 (user-defined state)<br>Norm of matrix ratio 2.05846e-08 difference 3.89655e-06 (constant state -1.0)<br>Norm of matrix ratio 5.77446e-09 difference 1.09307e-06 (constant state 1.0)"<br><br></div>with the option "-snes_type test -mat_fd_type ds" :<br><br>"Norm of matrix ratio 1.05213e-07 difference 1.95418e-05 (user-defined state)<br>Norm of matrix ratio 2.51463e-07 difference 4.76004e-05 (constant state -1.0)<br>Norm of matrix ratio 2.81081e-08 difference 5.32071e-06 (constant state 1.0)".<br><br></div>Could you give me any hint where should I check for the error?<span class="HOEnZb"><font color="#888888"><br><br></font></span></div><span class="HOEnZb"><font color="#888888">Khai<br><div><div><div><div><br> <br></div></div></div></div></font></span></div>
</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>