<div class="gmail_quote">On Tue, Apr 26, 2011 at 08:05, Tabrez Ali <span dir="ltr"><<a href="mailto:stali@purdue.edu">stali@purdue.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
I am trying to solve a system with constraints (0 on some diagonals). It works fine with AIJ but gives the following error (see below) with SBAIJ Matrices during KSPSolve (sequential). With GMRES it just segfaults.</blockquote>
</div><br><div>You cannot do LU with SBAIJ, try -pc_type cholesky.</div>