[petsc-users] addition of two matrices

[Umut Tabak]

Dear all,

How can I add two matrices? More openly, I have an expression

A^T * M* A  + B^T *K*B

I find M*A and K*B with MatMatMult then form the 1st and 2nd expression 
with MatMatMultTranspose and I would like to add the two eventually. I 
could not find a function in the documentation, maybe missed.

Any help is appreciated.

Umut