Hi, 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. I started with the ex5.c code, introducing just enough changes to model the new system. The analytic Jacobian function appears to be correct, with a Norm of matrix ratio < 1.e9(found using -snes_type test).
The problem I am having is that in most cases the SNES solver does not converge to a solution x of F(x)=0. 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 "stuck", i.e. it converges to a large non-zero result. Even though it is clearly not a solution, the -snes_converged_reason is reported as "Nonlinear solve converged due to CONVERGED_TR_DELTA". The intermediate KSP steps have -ksp_converged_reason reported as "Linear solve converged due to CONVERGED_STEP_LENGTH". 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. Does this sound like a familiar scenario with a familiar solution? Or can anyone point me to some documentation that describes the SNES tr and ls parameters in more detail than the manual.pdf? Or can anyone recommend the best SNES and KSP parameters for the Bratu example? Any help or advice would be greatly appreciated. Thanks, Sean Dettrick -------------- next part -------------- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20080107/7ec37516/attachment.htm>
