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>
