Hong <hzh...@mcs.anl.gov> writes: > Jed : > >> It is only implemented for SeqAIJ, but it should arguably have a default >> implementation that performs an explicit transpose and then calls >> MatMatMatMult or MatPtAP. >> > > We can implement > P= R^T (we have this for mpi matrices, expensive though) > C = R*A*P
Exactly.
signature.asc
Description: PGP signature