[petsc-dev] MatMultHermitianTranspose

[Jed Brown]

MatMultHermitianTranspose could have a default implementation in terms of
the others. Then it would work in all cases, right?

On Wed, Jul 18, 2012 at 7:59 AM, Jose E. Roman <jroman at dsic.upv.es> wrote:

> 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
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20120718/fc41fff5/attachment.html>