[petsc-dev] [petsc-users] How to compute RARt with A and R as distributed (MPI) matrices ?

[Jed Brown]

Hong <hzhang at mcs.anl.gov> writes:

> I add support for MatRARt_MPIAIJ
> https://bitbucket.org/petsc/petsc/commits/a3c14138ce0daf4ee55c7c10f1b4631a8ed2f13e

Is it better this way or as a fallback when !A->ops->rart?  MatPtAP
handles other combinations like MAIJ.

> It is in branch hzhang/mpirart.
> Let me know your comments.
>
> Hong
>
> On Wed, Jun 21, 2017 at 12:53 PM, Barry Smith <bsmith at mcs.anl.gov> wrote:
>
>>
>> > On Jun 21, 2017, at 8:00 AM, Franck Houssen <franck.houssen at inria.fr>
>> 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>
>>
>>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: signature.asc
Type: application/pgp-signature
Size: 832 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20170621/3098db0a/attachment.sig>