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

[Dave May]

You can assemble R^t and then use MatPtAP which supports MPIAIJ


On Wed, 21 Jun 2017 at 15:00, 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.
>
> 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
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20170621/07337b3f/attachment.html>