[petsc-dev] MatMultHermitianTranspose

Jed Brown jedbrown at mcs.anl.gov
Wed Jul 18 13:07:36 CDT 2012

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

> Er... well, now that I see better: it is not always safe, it will error
> out if mat->ops->conjugate is not available, for instance in BAIJ.

There is MatMultHermitianTranspose_SeqBAIJ so that would only an issue for
parallel BAIJ, but you were just going to use VecConjugate, right? In
Matlab notation, A'*x = (x'*A)', so you only need VecConjugate.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20120718/d816ee0d/attachment.html>

More information about the petsc-dev mailing list