<div dir="ltr">What is the first "-pc_type none" in your example?<div>Is it overridden by "-pc_type lu"?</div><div><br><div><br></div></div></div><div class="gmail_extra"><br><br><div class="gmail_quote">
2014-04-29 0:21 GMT-04:00 Barry Smith <span dir="ltr"><<a href="mailto:bsmith@mcs.anl.gov" target="_blank">bsmith@mcs.anl.gov</a>></span>:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<br>
ex10 -f0 ~/Datafiles/Matrices/Foad -ksp_type gmres -pc_type none -ksp_gmres_restart 2000 -ksp_plot_eigenvalues -draw_pause -1  -ksp_monitor_true_residual -ksp_rtol 1.e-6 -ksp_max_it 1990 -geometry 0,0,600,600 -draw_save -pc_type lu -ksp_rtol 1.e-12 -pc_factor_mat_solver_package superlu<br>

  0 KSP preconditioned resid norm 1.752626139140e+05 true resid norm 1.860478790525e+07 ||r(i)||/||b|| 1.000000000000e+00<br>
  1 KSP preconditioned resid norm 3.174750685865e+04 true resid norm 1.816271753408e+07 ||r(i)||/||b|| 9.762388921916e-01<br>
  2 KSP preconditioned resid norm 5.903085310063e-07 true resid norm 1.783619436500e-05 ||r(i)||/||b|| 9.586884008478e-13<br>
<br>
superlu does lots of pivoting for stability. What does it mean that the first iteration barely affects the residual and then the second iteration wipes it out?<br>
<br>
Right precondition wipes out the residual on one iteration<br>
<br>
gmres_restart 2000 -ksp_plot_eigenvalues -draw_pause -1  -ksp_monitor_true_residual -ksp_rtol 1.e-6 -ksp_max_it 1990 -geometry 0,0,600,600 -draw_save -pc_type lu -ksp_rtol 1.e-12 -pc_factor_mat_solver_package superlu -ksp_pc_side right<br>

  0 KSP unpreconditioned resid norm 1.860478790525e+07 true resid norm 1.860478790525e+07 ||r(i)||/||b|| 1.000000000000e+00<br>
  1 KSP unpreconditioned resid norm 2.559955879548e-07 true resid norm 2.505789584416e-07 ||r(i)||/||b|| 1.346852002386e-14<br>
Number of iterations =   1<br>
Residual norm 2.50579e-07<br>
<br>
<br>
</blockquote></div><br></div>