<div class="gmail_quote">On Thu, Dec 29, 2011 at 19:30, Dave Nystrom <span dir="ltr"><<a href="mailto:Dave.Nystrom@tachyonlogic.com">Dave.Nystrom@tachyonlogic.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div id=":6jr">Generally I use CG and LU from petsc.  Cholesky runs slower than LU and I was<br>
told in a previous email that it was because of the extra data movement for<br>
Cholesky versus LU.  I have also tried minres instead of CG and tend to get<br>
reduced iteration counts.  But I think minres may have an extra dot product<br>
and so is a little more expensive than cg.  But minres is usually cheaper<br>
over all because of the reduction in iterations.</div></blockquote></div><br><div>If you are concerned about solve time for serial runs, I would either use Cholesky from MUMPS or CHOLMOD (install suitesparse, then -pc_type cholesky -pc_factor_mat_solver_type cholmod).</div>