neckel at in.tum.de
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?
More information about the petsc-users