<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Tue, Jul 12, 2016 at 7:42 AM, Hoang Giang Bui <span dir="ltr"><<a href="mailto:hgbk2008@gmail.com" target="_blank">hgbk2008@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-style:solid;border-left-color:rgb(204,204,204);padding-left:1ex"><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></blockquote><div><br></div><div>1) In the log you sent, the linear solver converges due to the Relative Tolerance, 1.0e-9, not the breakdown tolerance 1e-30. Change rtol will affect the convergence.</div><div><br></div><div>2) What do you mean  -fieldsplit_wp_ksp_rtol 1.0e-8 does not work? ALWAYS send the view output.</div><div><br></div><div>3) I can't tell you anything about Newton convergence if you do not send the output, -snes_monitor -snes_view</div><div><br></div><div>4) If there is a difference between LU and an iterative solver with residual 1e-9, then your system is very ill-conditioned.</div><div><br></div><div>  Thanks,</div><div><br></div><div>     Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-style:solid;border-left-color:rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div>Sorry this problem is run with many time steps and is quite big so I cannot reproduce in a simple test case.</div><span class=""><font color="#888888"><div><br clear="all"><div><div data-smartmail="gmail_signature"><div dir="ltr">Giang</div><div dir="ltr"><br></div></div></div>
</div></font></span></div>
</blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature" data-smartmail="gmail_signature">What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>-- Norbert Wiener</div>
</div></div>