[Trilinos-Users] MatrixMatrix::Add and FillComplete

Simone Deparis deparis at MIT.EDU
Wed Oct 5 23:31:48 MDT 2005


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

More information about the Trilinos-Users mailing list