Try (e.g., -pc_type gamg) and additive Schwarz. Also try using ASM as a smoother.<div class="gmail_extra"><br><br><div class="gmail_quote">On Mon, Nov 12, 2012 at 10:21 AM, Rob Kunz <span dir="ltr"><<a href="mailto:rfk102@arl.psu.edu" target="_blank">rfk102@arl.psu.edu</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">





<div lang="EN-US" link="blue" vlink="purple">
<div>
<p class="MsoNormal">Hi:<u></u><u></u></p>
<p class="MsoNormal">I have a pressure based CFD code using PETSC for nearly SPD pressure equation.<u></u><u></u></p>
<p class="MsoNormal">PETSC default solver + preconditioner (GMRES+BJ for n-proc) is very reliable for me.<u></u><u></u></p>
<p class="MsoNormal">However I have a very tough bearing problem with extremely high aspect ratio cells and the 1-processor run works fine (GMRES+ILU0), but for nproc>1 the code does not converge.<u></u><u></u></p>
<p class="MsoNormal">Of course I can get the code to return identical 1-proc and 12-proc convergence histories if I choose a point preconditioner like Jacobi, but that kills the convergence of both.<u></u><u></u></p>
<p class="MsoNormal">Any recommendations?<u></u><u></u></p>
<p class="MsoNormal">Thanks<span class="HOEnZb"><font color="#888888"><br>
Rob Kunz<u></u><u></u></font></span></p>
</div>
</div>

</blockquote></div><br></div>