<div dir="ltr">I add support for MatRARt_MPIAIJ<div><a href="https://bitbucket.org/petsc/petsc/commits/a3c14138ce0daf4ee55c7c10f1b4631a8ed2f13e">https://bitbucket.org/petsc/petsc/commits/a3c14138ce0daf4ee55c7c10f1b4631a8ed2f13e</a><br></div><div><br></div><div>It is in branch hzhang/mpirart.</div><div>Let me know your comments. </div><div><br></div><div>Hong</div></div><div class="gmail_extra"><br><div class="gmail_quote">On Wed, Jun 21, 2017 at 12:53 PM, Barry Smith <span dir="ltr"><<a href="mailto:bsmith@mcs.anl.gov" target="_blank">bsmith@mcs.anl.gov</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><span class=""><br>
> On Jun 21, 2017, at 8:00 AM, Franck Houssen <<a href="mailto:franck.houssen@inria.fr">franck.houssen@inria.fr</a>> wrote:<br>
><br>
> How to compute RARt with A and R as distributed (MPI) matrices ?<br>
><br>
> This works with sequential matrices.<br>
> 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>
<br>
</span> Yes, when we say AIJ we mean both SeqAIJ and MPIAIJ. The manual page here is wrong, probably because it got copied from MatPtAP page<br>
<span class=""><br>
><br>
> Is this kind of matrix product possible with distributed matrices in PETSc ? Or is this a known limitation ?<br>
> Do I go the wrong way to do that (= should use another method) ? If yes, what is the correct one ?<br>
><br>
> Franck<br>
><br>
> PS: running debian/testing + gcc-6.3 + bitbucket petsc.<br>
><br>
> >> 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>
> >> mpirun -n 2 matRARt.exe mpi<br>
> [0]PETSC ERROR: --------------------- Error Message ------------------------------<wbr>------------------------------<wbr>--<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>
</span>> <matRARt.cpp><br>
<br>
</blockquote></div><br></div>