<div class="gmail_quote">On Wed, Dec 28, 2011 at 03:52, Alexander Grayver <span dir="ltr"><<a href="mailto:agrayver@gfz-potsdam.de">agrayver@gfz-potsdam.de</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div id=":1xd">I wrote simple program to test different solvers. It reads matrix and rhs and solves system. That is it.<br>
I tested this program using PETSc's interface to MUMPS (LDLt), SuperLU_Dist (LU) and it works well giving right solution, however when I switch<br>
to PaStiX (LLt) the solution is zero. </div></blockquote><div><br></div><div>How are you specifying to use PaStiX?</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div id=":1xd">The output is attached.<br>
The matrix is SPD and has condition number of 16, there should be no problem.<br>
I use petsc-dev.<br>
How could I figure out where is the problem?</div></blockquote></div><br><div>This sounds like a small example matrix. Can you -ksp_view_binary and send us the output file ('binaryoutput').</div>