<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN">
<html>
<head>
<meta content="text/html; charset=UTF-8" http-equiv="Content-Type">
</head>
<body bgcolor="#ffffff" text="#000000">
On 28.12.2011 13:55, Jed Brown wrote:
<blockquote
cite="mid:CAM9tzSkSOjS4oCz5tnZ2qh-RT+MrHPQY_ccSF=Ho9GjY6k22gg@mail.gmail.com"
type="cite">
<div class="gmail_quote">On Wed, Dec 28, 2011 at 03:52, Alexander
Grayver <span dir="ltr"><<a moz-do-not-send="true"
href="mailto:agrayver@gfz-potsdam.de">agrayver@gfz-potsdam.de</a>></span>
wrote:<br>
<blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt
0.8ex; border-left: 1px solid rgb(204, 204, 204);
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>
</blockquote>
<br>
Hi Jed,<br>
If you mean input parameters they are as following:<br>
-ksp_view -ksp_monitor -mat_type mpisbaij -ksp_type preonly -pc_type
cholesky -pc_factor_mat_solver_package pastix -mat_pastix_verbose 2
<br>
<br>
<blockquote
cite="mid:CAM9tzSkSOjS4oCz5tnZ2qh-RT+MrHPQY_ccSF=Ho9GjY6k22gg@mail.gmail.com"
type="cite">
<div class="gmail_quote">
<div> </div>
<blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt
0.8ex; border-left: 1px solid rgb(204, 204, 204);
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>
</blockquote>
<br>
Maybe it was too preliminary to say that it is positive definite
since CG diverges with reason DIVERGED_INDEFINITE_MAT, however GMRES
with default parameters converges successfully.<br>
Matrix is ~200000x200000, not so small, but zipped occupies 4.5 Mb.
Shall I send it?<br>
<br>
Regards,<br>
Alexander<br>
</body>
</html>