[petsc-users] MatPtAP for involving MPIDENSE and MPIAIJ matrices

[Hong]

I plan to implement MatTransposeMatMult_MPIDense_MPIDense via

1. add MatTransposeMatMult_elemental_elemental()
2. C_dense = P_dense^T * B_dense
    via MatConvert_dense_elemental() and MatConvert_elemental_dense()

Let me know if you have better suggestions.

Hong

On Thu, Oct 15, 2015 at 1:49 PM, Jed Brown <jed at jedbrown.org> wrote:

> Matthew Knepley <knepley at gmail.com> writes:
> > I do not understand "unrelated". My understanding was that TAIJ could be
> > used (with T = I)
> > to get the action of A on a set of vectors, which I think would be A P.
> Why
> > specifically would
> > you not use it?
>
> MAIJ is enough for that; it would involve packing P into a Vec, but
> Bikash might as well use MatMatMult to perform the same operation in a
> more direct way (MatMatMult_MPIAIJ_MPIDense exists).  The problem is
> that MatTransposeMatMult_MPIDense_MPIDense needs to be implemented.
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20151015/65dd92a6/attachment.html>