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

[Mark Adams]

>
>
> 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/18fd8fff/attachment.html>