[petsc-users] About MatTransposeMatMult

[Jed Brown]

Karl Rupp <rupp at mcs.anl.gov> writes:

> 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.

You can explicitly transpose X:

  MatTranspose(X,MAT_INITIAL_MATRIX,&Xt);

and then use MatMatMult.

>> 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

This is only for MPI.  Hong, can you expose MatTransposeMatMult_SeqAIJ_SeqDense?
-------------- next part --------------
A non-text attachment was scrubbed...
Name: not available
Type: application/pgp-signature
Size: 835 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20130926/3f3bb60d/attachment-0001.pgp>