[petsc-dev] MatTransposeMatMult

Jed Brown jed at jedbrown.org
Thu Sep 21 21:06:56 CDT 2017

Hong <hzhang at 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!

More information about the petsc-dev mailing list