<div dir="ltr"><br><div class="gmail_extra"><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex">
<div class=""><br>
<br>
</div>The most common cause of this is inaccurate eigenvalue bounds.  You can<br>
try using more iterations for estimation.<br>
<div class=""><br></div></blockquote><div><br></div><div>I have also found that CG converges much faster than GMRES (default) so if you are SPD, I would always use:  '-gamg_est_ksp_type cg' .  Also, '-pc_gamg_verbose 2' will print out the eigenvalues that GAMG computes for smoothing.  With -ksp_view you will see the eigenvalues that Cheby is using (GAMG sets these so they should be the same if you use GAMG with pc_gamg_nsmooths>0). The largest eigenvalue should be around 2-3.  If it is too low (<2) then the instability that Jed is referring to would be the problem.  </div>
<div><br></div><div>You can also use more iterations, and see if it has converged, with -gamg_est_ksp_max_it N.  I think the default is 5.</div></div></div></div>