> On Jun 21, 2017, at 8:00 AM, Franck Houssen 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 w
Hong writes:
> Jed :
>
>> It is only implemented for SeqAIJ, but it should arguably have a default
>> implementation that performs an explicit transpose and then calls
>> MatMatMatMult or MatPtAP.
>>
>
> We can implement
> P= R^T (we have this for mpi matrices, expensive though)
> C = R*A*P
Exac
> >
> >> De: "Dave May"
> >> À: "Franck Houssen" , "PETSc users list"
> >> , "petsc-dev"
> >> Envoyé: Mercredi 21 Juin 2017 15:11:27
> >> Objet: Re: [petsc-users] How to compute RARt with A and
ages are the same with R replaced with Pt). Where was it possible to
> find this information ?
>
> Franck
>
> - Mail original -
>
>> De: "Dave May"
>> À: "Franck Houssen" , "PETSc users list"
>> , "petsc-de
t"
> , "petsc-dev"
> Envoyé: Mercredi 21 Juin 2017 15:11:27
> Objet: Re: [petsc-users] How to compute RARt with A and R as distributed
> (MPI) matrices ?
> You can assemble R^t and then use MatPtAP which supports MPIAIJ
> On Wed, 21 Jun 2017 at 15:00, Franck H
You can assemble R^t and then use MatPtAP which supports MPIAIJ
On Wed, 21 Jun 2017 at 15:00, Franck Houssen
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
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
