<div dir="ltr">Hello<div><br></div><div>I encountered different convergence behaviour of Newton Raphson when using different solver settings with PETSc</div><div><br></div><div>For the first solver configuration, I used direct solver</div><div><div>-ksp_type preonly</div><div>-pc_type lu</div><div>-pc_factor_mat_solver_package mumps</div><div>-mat_mumps_icntl_1 6</div><div>-mat_mumps_icntl_4 3</div><div>-mat_mumps_icntl_7 4</div><div>-mat_mumps_icntl_14 40</div><div>-mat_mumps_icntl_23 0</div></div><div><br></div><div>The simulation can run completely and the NR typically converged after 6/7 iterations. Of course, it's very slow. For the second solver configuration:</div><div><div>-ksp_type gmres</div><div>-ksp_max_it 300</div><div>-ksp_gmres_restart 300</div><div>-ksp_gmres_modifiedgramschmidt</div><div>-pc_view</div><div>-pc_fieldsplit_type multiplicative</div><div>-fieldsplit_u_pc_type hypre</div><div>-fieldsplit_u_pc_hypre_type boomeramg</div><div>-fieldsplit_u_pc_hypre_boomeramg_coarsen_type PMIS</div><div>-fieldsplit_u_pc_hypre_boomeramg_strong_threshold 0.6</div><div>-fieldsplit_u_pc_hypre_boomeramg_max_levels 25</div><div>-fieldsplit_wp_ksp_rtol 1.0e-8</div><div>-fieldsplit_wp_pc_type hypre</div><div>-fieldsplit_wp_pc_hypre_type boomeramg</div><div>-fieldsplit_wp_pc_hypre_boomeramg_coarsen_type PMIS</div><div>-fieldsplit_wp_pc_hypre_boomeramg_strong_threshold 0.6</div><div>-fieldsplit_wp_pc_hypre_boomeramg_max_levels 25</div></div><div><br></div><div>The solver runs much faster, but the NR does not converge in 30 iterations after some time steps. I thought setting the solver tolerance -ksp_rtol 1.0e-12 but it doesn't help much because GMRES already terminate with tolerance 1e-30 (see sample log file). Can we set the tolerance of the sub-ksp of the Fieldsplit? I tried -fieldsplit_wp_ksp_rtol 1.0e-8 but it doesn't work.</div><div><br></div><div>Sorry this problem is run with many time steps and is quite big so I cannot reproduce in a simple test case.</div><div><br clear="all"><div><div class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr">Giang</div><div dir="ltr"><br></div></div></div>
</div></div>