Please, include -snes_view so that we can see exactly what solver options are being used. Dmitry.
On Mon, Dec 9, 2013 at 12:48 PM, John Peterson <[email protected]> wrote: > > On Mon, Dec 9, 2013 at 11:31 AM, Lorenzo Zanon > > <[email protected]>wrote: > > > >> I get faster convergence with e.g. -snes_linesearch_type basic -ksp_rtol > >> 1e-4: > >> > >> NL step 0, |residual|_2 = 3.464102e-05 > >> NL step 1, |residual|_2 = 2.417540e-04 > >> NL step 2, |residual|_2 = 6.174706e-08 > >> NL step 3, |residual|_2 = 3.577768e-12 > >> NL step 4, |residual|_2 = 7.687278e-17 > >> StVen system solved at nonlinear iteration 4 , final nonlinear residual > >> norm: 7.687278e-17 > > Your ksp_rtol might be hurting your convergence a bit near the root... > > Have you tried -snes_ksp_ew? There are several EW-specific options > within you can play with as well: > > > http://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/SNES/SNESSetFromOptions.html > > -- > John > > > ------------------------------------------------------------------------------ > Sponsored by Intel(R) XDK > Develop, test and display web and hybrid apps with a single code base. > Download it for free now! > > http://pubads.g.doubleclick.net/gampad/clk?id=111408631&iu=/4140/ostg.clktrk > _______________________________________________ > Libmesh-users mailing list > [email protected] > https://lists.sourceforge.net/lists/listinfo/libmesh-users > -- Dmitry Karpeev Mathematics and Computer Science Argonne National Laboratory Argonne, Illinois, USA and Computation Institute University of Chicago 5735 S. Ellis Avenue Chicago, IL 60637 ----------------------- Phone: 630-252-1229 Fax: 630-252-5986 ------------------------------------------------------------------------------ Sponsored by Intel(R) XDK Develop, test and display web and hybrid apps with a single code base. Download it for free now! http://pubads.g.doubleclick.net/gampad/clk?id=111408631&iu=/4140/ostg.clktrk _______________________________________________ Libmesh-users mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-users
