[petsc-users] Compute C*Ct using MatCreateTranspose for Ct

[Olivier Jamond]

Ahah, maybe I am using a too old version. I use the 3.12.3. I will try 
with a more recent one.

On 24/09/2020 12:47, Mark Adams wrote:
>
>
>     Is there a way to avoid the explicit transposition of the matrix?
>
>
> It does not look like we have A*B^T for mpiaij as the error message 
> says. I am not finding it in the code.
>
> Note, MatMatMult with a transpose shell matrix, I suspect that it does 
> an explicit transpose internally, or it could notice that you have 
> C^T*C and we might have that implemented in-place (I doubt it, but it 
> would be legal and fine to do).
>
>     Many thanks,
>     Olivier Jamond
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20200924/22d372ff/attachment.html>