[petsc-users] About MatTransposeMatMult

[Joon Hee Choi]

Thank you. I think MatMatMult runs well for multiplication between seqaij and seqdense.

Joon

----- Original Message -----
From: "Jed Brown" <jedbrown at mcs.anl.gov>
To: "Karl Rupp" <rupp at mcs.anl.gov>, "Joon Hee Choi" <choi240 at purdue.edu>, petsc-users at mcs.anl.gov
Cc: "Hong Zhang" <hzhang at mcs.anl.gov>
Sent: Thursday, September 26, 2013 11:01:42 AM
Subject: Re: [petsc-users] About MatTransposeMatMult

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