<div class="gmail_quote">On Wed, Dec 28, 2011 at 09:10, 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 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:0 0 0 .8ex;border-left:1px #ccc solid;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. 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>