[petsc-users] inplace MatPtAP

Patrick Lacasse patrick.m.lacasse at gmail.com
Tue Apr 28 13:31:58 CDT 2015

I have to compute a product Pt*A*P.

My projection matrix P won't change the structure of A and I don't need to
keep the original A.
Is it possible to do it inplace (modify A instead of creating a new C)?

Do you think it is feasible to expend MatPtAP to mpiaij and/or sbaij

There is more implementation of MatMatMult.
Is it possible to do Pt*A and A*P directly into A (without creating a copy)?

Patrick Lacasse
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20150428/ea3607f2/attachment.html>

More information about the petsc-users mailing list