[petsc-users] About MatTransposeMatMult

Karl Rupp rupp at mcs.anl.gov
Thu Sep 26 04:55:24 CDT 2013


Hi,

> I am trying to calculate M = X^T * B. X is a seqaij and B is a seqdense matrix.

This is not yet allowed in the latest 3.4. release.

> However, I used MatTransposeMatMult() but I got an error. Can someone why I got the error and how I should modify my code? I am attaching my simple code, error message, and result. Thank you very much.

Do you use the master branch from our Bitbucket repository? It contains 
the commit where Hong added the respective functionality:
https://bitbucket.org/petsc/petsc/commits/8949adfd119cb1d44c9fb5cf3dc01f4b6a02252d

Best regards,
Karli


More information about the petsc-users mailing list