Re: [petsc-users] Exploiting symmetry with direct solvers

2015-04-09 Thread pascal . tremblay
té D3- Conservation <90j De :David Knezevic A : Hong , Cc :PETSc users list Date : 09/04/2015 04:17 Objet : Re: [petsc-users] Exploiting symmetry with direct solvers Envoyé par : On Wed, Apr 8, 2015 at 9:28 PM, Hong wrote: David , I'm curious about how to do a sy

Re: [petsc-users] Exploiting symmetry with direct solvers

2015-04-08 Thread David Knezevic
On Wed, Apr 8, 2015 at 11:06 PM, Hong wrote: > David: > >> Mumps supports Cholesky factorization for symmetric, and symmetric+spd >>> matrices. You may consult mumps user manual. >>> >> >> >> OK, thanks. But I was also wondering about the PETSc interface to MUMPS. >> >> I consulted the MUMP

Re: [petsc-users] Exploiting symmetry with direct solvers

2015-04-08 Thread Hong
David: > >>> Mumps supports Cholesky factorization for symmetric, and symmetric+spd >> matrices. You may consult mumps user manual. >> > > > OK, thanks. But I was also wondering about the PETSc interface to MUMPS. > > I consulted the MUMPS manual, but it just says that you need to specify > "SYM=

Re: [petsc-users] Exploiting symmetry with direct solvers

2015-04-08 Thread David Knezevic
On Wed, Apr 8, 2015 at 9:28 PM, Hong wrote: > David , > >> 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 understandin

Re: [petsc-users] Exploiting symmetry with direct solvers

2015-04-08 Thread Hong
David , > 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

[petsc-users] Exploiting symmetry with direct solvers

2015-04-08 Thread 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: MatSetOptio