[Trilinos-Users] MatrixMatrix & repeated indeces

[Simone Deparis]

Hi,

just FYI: I did very small tests using MatrixMatrix:: sum and multiply
using repeated indices.
My experience is:

let A be a CsrMatrix whose rowMap has repeated indeces
let B be a CsrMatrix whose rowMap has unique indeces.

if C has the same rowMap as B:
sum(A,..,C,..) does _not_ work
multiply(A,..,A,...,C) does not work because of the repeated indices

if C has the same rowMap as A:
sum(A,..,C,..) works
multiply(A,..,A,...,C) does not work because of the repeated indices

Simone


Michael A Heroux wrote:
> Simone,
> 
> Good to hear this approach works.  
> 
> In EpetraExt (the Epetra extensions package) there is a MatrixMatrix
> subpackage.  It supports matrix addition.  If the operation you want is not
> supported, then it is fairly easy to write a simple loop expression to
> accomplish the task.