<!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 16:27, Jed Brown wrote:
<blockquote
cite="mid:CAM9tzS=0MPZs5wZotkAZpB7Q7YgkHxideNvdbipgxTrr9JGUCQ@mail.gmail.com"
type="cite">
<div class="gmail_quote">On Wed, Dec 28, 2011 at 09:10, 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 bgcolor="#ffffff" text="#000000">So if I understood it
right, I put diff you gave in file pastix_patch in petsc-dev
dir and called:<br>
<br>
lib/petsc-dev> patch -Np1 < pastix_patch<br>
<br>
And built petsc again. Now <br>
<br>
lib/petsc-dev> hg diff
src/mat/impls/aij/mpi/pastix/pastix.c </div>
</blockquote>
<div><br>
</div>
<div>Yes, the patch applied correctly.</div>
<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 bgcolor="#ffffff" text="#000000">
<div>
<div class="h5"><br>
</div>
</div>
So I guess it's patched correctly?<br>
Because I don't see much difference in output. Please find
it attached. <br>
Did I make something wrong?</div>
</blockquote>
</div>
<br>
<div>Hmm, it is hard-coded to say that the matrix is symmetric. We
are setting API_FACT_LU (you could check that it has not been
changed by attaching a debugger in MatFactorNumeric_PaStiX and
printing lu->iparm[IPARM_FACTORIZATION]). Maybe PaStiX is
ignoring this and always doing LLt instead. </div>
</blockquote>
<br>
No. I checked, the API_FACT_LU is set correctly. <br>
<br>
<blockquote
cite="mid:CAM9tzS=0MPZs5wZotkAZpB7Q7YgkHxideNvdbipgxTrr9JGUCQ@mail.gmail.com"
type="cite">
<div>We probably need to set API_PARM_LDLT (and possibly
API_PARM_LU, but they claim that only does static pivoting,
which might still be a problem). You could try adding this to
src/mat/impls/aij/mpi/pastix/pastix.c or ask the PaStiX
developers what they suggest for symmetric indefinite systems.</div>
</blockquote>
<br>
I tried API_PARM_LDLT, but it has no effect. <br>
<br>
Thanks for help, I will consider writing PaStiX team. <br>
<br>
Regards,<br>
Alexander<br>
</body>
</html>