[petsc-dev] PaStiX solution

[Jed Brown]

On Wed, Dec 28, 2011 at 03:52, Alexander Grayver <agrayver at gfz-potsdam.de>wrote:

> I wrote simple program to test different solvers. It reads matrix and rhs
> and solves system. That is it.
> 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
> to PaStiX (LLt) the solution is zero.
>

How are you specifying to use PaStiX?


> The output is attached.
> The matrix is SPD and has condition number of 16, there should be no
> problem.
> I use petsc-dev.
> How could I figure out where is the problem?
>

This sounds like a small example matrix. Can you -ksp_view_binary and send
us the output file ('binaryoutput').
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20111228/70945acf/attachment.html>