Re: [petsc-users] MUMPS with symmetric matrices

Gianluca Meneghello writes: > Dear all, > > I am trying to solve a linear system for a symmetric matrix with MUMPS. > > Is there a way to tell MUMPS that the matrix is indeed symmetric? > > The way I build the matrix is > > Mat A,AT,ATA >

Re: [petsc-users] MUMPS with symmetric matrices

Gianluca Meneghello writes: > That is correct... I will try with -pc_type cholesky and use > MatTransposeMatMult. > > Using cholesky I do not need to specify mumps as a solver, am I right? Of course you do. > A is a linearization of the Navier Stokes equation. Of the

Re: [petsc-users] MUMPS with symmetric matrices

That is correct... I will try with -pc_type cholesky and use MatTransposeMatMult. Using cholesky I do not need to specify mumps as a solver, am I right? A is a linearization of the Navier Stokes equation. Thanks! Gianluca On Wed, Nov 4, 2015 at 12:46 PM, Jed Brown wrote:

Re: [petsc-users] MUMPS with symmetric matrices

It is a discretization of the differential operator, of which I would need the inverse (or LU decomposition). My goal is frequency response (resolvant) analysis of the linearized Navier-Stokes operator. There was a reason I was not using MatTransposeMatMult, that is the matrix is complex and I

Re: [petsc-users] MUMPS with symmetric matrices

I have just read that there is no special algorithm for Hermitian matrices in MUMPS (sorry, I meant Hermitian, not symmetric... the matrix is complex). Sorry for this. In any case, if there is any suggestion it is more than welcome! Thanks for your help and your work, Gianluca On Wed, Nov 4,

[petsc-users] MUMPS with symmetric matrices

Dear all, I am trying to solve a linear system for a symmetric matrix with MUMPS. Is there a way to tell MUMPS that the matrix is indeed symmetric? The way I build the matrix is Mat A,AT,ATA MatHermitianTranspose(A,MAT_INITIAL_MATRIX,); MatMatMult(AT,A,MAT_INITIAL_MATRIX,7,);