Hi,<br>
<br>
I am solving a version of the Grad-Shafranov equation, F(x)=0, which but for some extra spatial dependences is similar in form to the 2D Bratu equation in snes/examples/tutorials/ex5.c.&nbsp; I started with the ex5.c code, introducing just enough changes to model the new system.&nbsp; The analytic Jacobian function appears to be correct, with a Norm of matrix ratio &lt; 
1.e9 (found using -snes_type test).<br><br>The problem I am having is that in most cases the SNES solver does not converge to a solution x of F(x)=0.&nbsp; Rather, what happens is that the fnorm (obtained in a monitor function) converges to some large non-zero value, and F(x) seems to get &quot;stuck&quot;, 
i.e. it converges to a large non-zero result.&nbsp; Even though it is clearly not a solution, the&nbsp; -snes_converged_reason is reported as &quot;Nonlinear solve converged due to CONVERGED_TR_DELTA&quot;.&nbsp; The intermediate KSP steps have -ksp_converged_reason reported as &quot;Linear solve converged due to CONVERGED_STEP_LENGTH&quot;.&nbsp; I have been typically running with parameters -da_grid_x 100 -da_grid_y 101 -snes_converged_reason -ksp_converged_reason -snes_type tr -ksp_type gmres -snes_max_it 100.
<br><br>Does this sound like a familiar scenario with a familiar solution?&nbsp; Or can anyone point me to some documentation that describes the SNES tr and ls parameters in more detail than the manual.pdf?<br>Or can anyone recommend the best SNES and KSP parameters for the Bratu example?
<br><br>Any help or advice would be greatly appreciated.<br><br>Thanks,<br>Sean Dettrick<br>