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 [[email protected]] 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>
