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

[Jed Brown]

Hong <hzhang at mcs.anl.gov> 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

Exactly.
-------------- 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-users/attachments/20170621/09052c06/attachment.pgp>