Re: [petsc-users] Inertia of Hermitian Matrix

2018-03-08 Thread Adam Denchfield
You can verify there is a generalization of the inertia theorem thanks to Ikramov. Furthermore, recall the inertia theorem is about eigenvalues; for a Hermitian matrix, you can diagonalize it in a form that has real entries, and that preserves the eigenvalues. Of course, being able to apply the

Re: [petsc-users] Inertia of Hermitian Matrix

2018-02-23 Thread Jose E. Roman
Unfortunately MUMPS does not return inertia with complex matrices, it seems that it is not implemented. See the note "Usage with Complex Scalars" in section 3.4.5 of SLEPc users manual. You could use multi-communicators with as many partitions as MPI processes, so that each process performs a

[petsc-users] Inertia of Hermitian Matrix

2018-02-22 Thread Anthony Ruth
Hello, I am trying to diagonalize a hermitian matrix using the Eigen Problem Solver in SLEPc, I run into errors on calls to MatGetInertia() with complex hermitian matrices that I did not see with real matrices. The complex and real versions were done with separate PETSC_ARCH. I do not know if the