<div class="gmail_quote">On Tue, May 10, 2011 at 13:44, Danesh Daroui <span dir="ltr"><<a href="mailto:danesh.daroui@ltu.se">danesh.daroui@ltu.se</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
<div id=":74f">I also checked the SPAI web page and they have also mentioned that there<br>
is an interface for PETSc but nothing more.</div></blockquote><div><br></div><div>Configure with --download-spai. You might also try parasails from Hypre (another SPAI implementation), configure with --download-hypre, then run with -pc_type hypre -pc_hypre_type parasails.</div>
<div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;"><div id=":74f">
<br>
I have also another questions. In my code, even when the solution is<br>
converged withing few iterations, the result is far different from<br>
correct solution that I have. I am suspecting that there something wrong<br>
with matrix filling in my code. Is there anyway to use direct solver in<br>
PETSc so I will verify that my matrices are correct before trying<br>
iterative solvers?</div></blockquote></div><br><div>-pc_type lu</div>