[petsc-users] question about special matrix-matrix-matrix product

[Kubis, Tillmann C]

Hello,
I need the diagonal of GR*Gamma*GR^dagger , where Gamma is a sparse and GR is a dense matrix. Is there a good/already coded way to do that?
I am up to write most of the product myself using the method MatGetRow on GR and Gamma and multiplying the elements directly. So far these are serial complex matrices, but all matrices will be distributed in the end.
Thanks,
Tillmann Kubis
____________________________________
Tillmann Kubis, PhD
Research Assistant Professor
Network for Computational Nanotechnology
207 S Martin Jischke Drive
Purdue University, DLR, room 441-1
West Lafayette, Indiana 47907-1971
phone: +1-765-496-7312
fax: +1-765-496-6026

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20130729/10e76637/attachment.html>