[Trilinos-Users] MatrixMatrix::Add and FillComplete

[Simone Deparis]

Hi,

In the method
EpetraExt::MatrixMatrix::Add(A, transposeA, scalarA, B, scalarB )
the matrix B is the one to which A is added
(kind of B*=scalarB; B+=scalarA*A).
at input, B may be filled or not.

Is there a reason why, if B is not filled, then B is filled at the end?
I would have expected that the FillComplete is left to the user, also to 
allow the choice of to FillComplete(map1,map2) or to perform another add 
where the entries are not yet fixed.

Does somebody have an explanation, please?

thank you very much
Best
Simone