[petsc-dev] MatMultHermitianTranspose

Jose E. Roman jroman at dsic.upv.es
Wed Jul 18 07:59:04 CDT 2012


Now the SVD in SLEPc correctly applies complex conjugation in complex scalar builds. We have two options:
- Implicit: uses MatMultHermitianTranspose
- Explicit: uses MatTranspose plus MatConjugate (the default)

The problem is that MatMultHermitianTranspose is defined only for SEQBAIJ (not for parallel BAIJ nor for AIJ sequential or parallel). So the implicit option is quite limited now for complex.

Conversely, the explicit option works in all cases except for BAIJ matrices since MatConjugate is not defined for BAIJ.

How difficult is to implement these operations?

Jose




More information about the petsc-dev mailing list