<div>You can assemble R^t and then use MatPtAP which supports MPIAIJ </div><div><br></div><div><br><div class="gmail_quote"><div>On Wed, 21 Jun 2017 at 15:00, Franck Houssen <<a href="mailto:franck.houssen@inria.fr">franck.houssen@inria.fr</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><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></div></blockquote></div></div>