[petsc-users] Exploiting symmetry with direct solvers

[David Knezevic]

I'm curious about how to do a symmetric LDL^T factorization (instead of LU)
with MUMPS and SuperLU. Based on this example:

http://www.mcs.anl.gov/petsc/petsc-3.4/src/ksp/ksp/examples/tutorials/ex52.c.html

my understanding is as follows:

- With MUMPS I gather that we need to specify:

MatSetOption(A,MAT_SPD,PETSC_TRUE);
PCSetType(pc,PCCHOLESKY);

I guess "-pc_type cholesky" on the command line is equivalent to the
PCSetType call, right? Is specifying MAT_SPD required in order for MUMPS to
do an LDL^T factorization?

- I gather that SuperLU doesn't provide a symmetric factorization.

Is this correct?

Thanks!
David
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20150408/63b8e975/attachment.html>