[Trilinos-Users] MatrixMatrix::Add and FillComplete

Williams, Alan B william at sandia.gov
Thu Oct 6 07:21:23 MDT 2005


Simone,

I can't recall any good reason for MatrixMatrix::Add to call
FillComplete on the B matrix, and in fact calling FillComplete without
specifying the domain and range maps is an error if the matrix is
rectangular.

I will remove the call to B.FillComplete().

Alan


> -----Original Message-----
> From: trilinos-users-bounces at software.sandia.gov 
> [mailto:trilinos-users-bounces at software.sandia.gov] On Behalf 
> Of Simone Deparis
> Sent: Wednesday, October 05, 2005 11:32 PM
> To: trilinos-users at software.sandia.gov
> Subject: [Trilinos-Users] MatrixMatrix::Add and FillComplete
> 
> 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
> 
> _______________________________________________
> Trilinos-Users mailing list
> Trilinos-Users at software.sandia.gov
> http://software.sandia.gov/mailman/listinfo/trilinos-users
> 



More information about the Trilinos-Users mailing list