<html><body><div style="font-family: times new roman, new york, times, serif; font-size: 12pt; color: #000000"><div>How to compute RARt with A and R as distributed (MPI) matrices ?</div><div><br></div><div>This works with sequential matrices.</div><div>The doc say "currently only implemented for pairs of AIJ matrices and classes which inherit from AIJ": I supposed that MPIAIJ was someway inheriting from AIJ, seems that it doesn't.<br></div><div><br></div><div>Is this kind of matrix product possible with distributed matrices in PETSc ? Or is this a known limitation ?</div><div>Do I go the wrong way to do that (= should use another method) ? If yes, what is the correct one ?<br></div><div><br></div><div>Franck<br></div><div><br></div><div>PS: running debian/testing + gcc-6.3 + bitbucket petsc.<br></div><div><br></div><div>>> mpirun -n 2 matRARt.exe seq<br>Mat Object: 1 MPI processes<br> type: seqaij<br>row 0: (0, 1.) (1, 0.) <br>row 1: (0, 0.) (1, 1.) <br><br></div><div>>> mpirun -n 2 matRARt.exe mpi<br>[0]PETSC ERROR: --------------------- Error Message --------------------------------------------------------------<br>[0]PETSC ERROR: No support for this operation for this object type<br>[0]PETSC ERROR: Matrix of type <mpiaij> does not support RARt<br><br></div></div></body></html>