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>
