[Trilinos-Users] Problem with EpetraExt::MatrixMatrix::Multiply

[bender calculon]

Hi,

I think that my problem is very simple but I don't find the way to resolve it.

When I use 

EpetraExt::MatrixMatrix::Multiply(D,false,I,false,B_jor);   

[for example D = diagonal matrix, I = sparse matrix ]

in B_jor i don't have an identity matrix but only the first element equal to one. So i think that i have some problem whit EpetraExt::MatrixMatrix::Multiply.


The matrix I is generate by 

Trilinos_Util::CrsMatrixGallery Matrix(CLP.Get("-problem", "laplace_2d"), Comm);
Epetra_CrsMatrix I = Matrix.GetMatrixRef(); 
Epetra_Map Map = Matrice.GetMapRef(); 

and the matrix D by 

Epetra_CrsMatrix D(I);
Epetra_Vector Diag(Map);
D.PutScalar(0);
I.ExtractDiagonalCopy (Diag);
Diag.Reciprocal(Diag);    
D.ReplaceDiagonalValues(Diag);

Thanks

Mau

_________________________________________________________________
Con Windows Live, puoi organizzare, modificare e condividere le tue foto.
http://www.microsoft.com/italy/windows/windowslive/products/photo-gallery-edit.aspx
-------------- next part --------------
An HTML attachment was scrubbed...
URL: https://software.sandia.gov/pipermail/trilinos-users/attachments/20090901/2f6084e4/attachment.html