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
