Hong <hzhang@mcs.anl.gov> writes:

> It was implemented a decade ago for small to medium matrices.
> Your matrix structure reveals the worst part of outer product.
>
> After discussing it with Barry, we decided to replace the algorithm (not
> API) using At=A^T and C=At*B  for sequential implementation.
> In parallel, I'll do local MatTranspose (used for PtAP).
> I'll give you this new implementation in few days.

Perfect, thanks!