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