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 <[email protected]> wrote:
> Matthew Knepley <[email protected]> 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.
>
