Dear all,
I have run into a problem today with the gmm++ library. I am not very
sure where the problem could be. I coded a symmetric lanczos solver for
my problem and used gmm++ and MUMPS (as the linear solver). However, the
accuracy of the basic linear algebra operations matrix-vector product
and gmm::add, scale and so seem to be questionable to me.
Well it could be a bit cryptic to understand for you, maybe, but let me
rephrase with an example:
I am computing some coefficients to form a tridiagonal matrix and these
coefficients that I compute in gmm++ start to deviate from the ones I
computed in MATLAB. They seem to be the same at the starting iterations
then conducting these matrix-vector operations(for orthogonalizations),
they start to diverge, I am not sure what the problem could be.
I went over my code several times, however, the end results tell me that
there is indeed a problem in the orthogonalization phase, which is
basically a combination of matrix-vector operations.
Any comments are appreciated on this issue.
Best regards,
Umut
_______________________________________________
Getfem-users mailing list
[email protected]
https://mail.gna.org/listinfo/getfem-users
