Matrix addition

Tobias Neckel neckel at
Thu Dec 28 07:02:45 CST 2006


I have a small (perhaps even silly) question concerning an addition of two matrices. 
Sorry in advance since I am quite new in using petsc.

Currently I am successfully using MatPtAP for the computation of a matrix Q = M^T*A*M 
with matrices M and A in my C++ code.
Now, I need to enlarge the functionality to use two matrices M1 and M2 together with 
A in order to calculate Q = M1^T*A*M1 + M2^T*A*M2.

I have been looking in the manual pages for a routine for adding up two matrices 
without success. Is there any way to do this and I just didn't see it? Or is there 
another possibility in order to compute my desired matrix Q?

Thank you

More information about the petsc-users mailing list