<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div text="#000000" bgcolor="#FFFFFF"><br>
    Or do I have to use KSPBCGS or KSPGMRES, which is directly from
    PETSc? However, I ran KSPGMRES yesterday with the Poisson eqn and my
    ans didn't converge.<br></div></blockquote><div><br></div><div>As Matt said GMRES is not great for symmetric operators like Poisson and you can use CG for the KSP method. HYPRE and GAMG are both fine preconditioners. Try them both.</div><div><br></div><div>And I would suggest not worrying about scaling until you get a solver that converges. Get your solver working well on small problems and as long as you are not relying on fundamentally unscalable algorithms you should be able to get scalability later.</div></div></div></div>