I add support for MatRARt_MPIAIJ https://bitbucket.org/petsc/petsc/commits/a3c14138ce0daf4ee55c7c10f1b4631a8ed2f13e
It is in branch hzhang/mpirart. Let me know your comments. Hong On Wed, Jun 21, 2017 at 12:53 PM, Barry Smith <[email protected]> wrote: > > > On Jun 21, 2017, at 8:00 AM, Franck Houssen <[email protected]> > wrote: > > > > How to compute RARt with A and R as distributed (MPI) matrices ? > > > > This works with sequential matrices. > > 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. > > 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 > > > > > Is this kind of matrix product possible with distributed matrices in > PETSc ? Or is this a known limitation ? > > Do I go the wrong way to do that (= should use another method) ? If yes, > what is the correct one ? > > > > Franck > > > > PS: running debian/testing + gcc-6.3 + bitbucket petsc. > > > > >> mpirun -n 2 matRARt.exe seq > > Mat Object: 1 MPI processes > > type: seqaij > > row 0: (0, 1.) (1, 0.) > > row 1: (0, 0.) (1, 1.) > > > > >> mpirun -n 2 matRARt.exe mpi > > [0]PETSC ERROR: --------------------- Error Message > -------------------------------------------------------------- > > [0]PETSC ERROR: No support for this operation for this object type > > [0]PETSC ERROR: Matrix of type <mpiaij> does not support RARt > > > > <matRARt.cpp> > >
