<div class="gmail_quote">On Wed, Feb 16, 2011 at 16:33, Juha Jäykkä <span dir="ltr">&lt;<a href="mailto:juhaj@iki.fi">juhaj@iki.fi</a>&gt;</span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
<div id=":4z4">Sorry, I looked at a wrong output. For the exact solution, it is the line<br>
search, which fails. KSP converges with CONVERGED_RTOL, but SNES quits with<br>
DIVERGED_LS_FAILURE.</div></blockquote></div><br><div>Your Jacobian may be incorrect, try running with -snes_mf_operator -pc_type lu -ksp_monitor. The linear solves should converge in 1 iteration. You can also try -mat_mffd_type ds if the residual is ill-conditioned.</div>
<div><br></div><div>If you can make the problem small, use -snes_type test to check the correctness of the Jacobian directly and -snes_test_display to show how the entries.</div>