Karl Rupp <[email protected]> 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?
pgpMyRWPEqdmM.pgp
Description: PGP signature
