<div dir="ltr">I plan to implement MatTransposeMatMult_MPIDense_MPIDense via <div><br></div><div>1. add MatTransposeMatMult_elemental_elemental()</div><div>2. C_dense = P_dense^T * B_dense</div><div>    via MatConvert_dense_elemental() and MatConvert_elemental_dense()</div><div><br></div><div>Let me know if you have better suggestions.</div><div><br></div><div>Hong<br><div class="gmail_extra"><br><div class="gmail_quote">On Thu, Oct 15, 2015 at 1:49 PM, Jed Brown <span dir="ltr"><<a href="mailto:jed@jedbrown.org" target="_blank">jed@jedbrown.org</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex"><span class="">Matthew Knepley <<a href="mailto:knepley@gmail.com">knepley@gmail.com</a>> writes:<br>
> I do not understand "unrelated". My understanding was that TAIJ could be<br>
> used (with T = I)<br>
> to get the action of A on a set of vectors, which I think would be A P. Why<br>
> specifically would<br>
> you not use it?<br>
<br>
</span>MAIJ is enough for that; it would involve packing P into a Vec, but<br>
Bikash might as well use MatMatMult to perform the same operation in a<br>
more direct way (MatMatMult_MPIAIJ_MPIDense exists).  The problem is<br>
that MatTransposeMatMult_MPIDense_MPIDense needs to be implemented.<br>
</blockquote></div><br></div></div></div>