[petsc-users] question about MatMatMultTranspose

[Gao Bin]

Hi Jed,

Thank you for your quick reply. But as pointed out at http://www.mcs.anl.gov/petsc/petsc-dev/docs/manualpages/Mat/MatMatTransposeMult.html:

This routine is currently only implemented for pairs of SeqAIJ matrices. C will be of type MATSEQAIJ.

Therefore I can not use it for dense matrix, am I right? If so, will MatMatTransposeMult be extended for other types of matrix later on? Thank you very much.

Cheers

Gao
________________________________
From: petsc-users-bounces at mcs.anl.gov [petsc-users-bounces at mcs.anl.gov] on behalf of Jed Brown [jedbrown at mcs.anl.gov]
Sent: Thursday, April 05, 2012 2:09 PM
To: PETSc users list
Subject: Re: [petsc-users] question about MatMatMultTranspose

On Thu, Apr 5, 2012 at 04:43, Gao Bin <bin.gao at uit.no<mailto:bin.gao at uit.no>> wrote:
I would like to use PETSc perform matrix-matrix multiplication C=A*B^T. I notice that MatMatMultTranspose could calculate C=A^T*B, but which is not I want.

Does anyone know how to perform C=A*B^T using PETSc (I know I could first calculate B^T and call MatMatMult with A and B^T, but it might not be efficient, is it?)? Thank you in advance.

Use petsc-dev. That routine is renamed (http://www.mcs.anl.gov/petsc/petsc-dev/docs/manualpages/Mat/MatTransposeMatMult.html) and the one you want is

http://www.mcs.anl.gov/petsc/petsc-dev/docs/manualpages/Mat/MatMatTransposeMult.html
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20120405/e2c9a14d/attachment.htm>