Hello numerical analysis experts,
I am in need of Fortran 90 code that implements the Lanzcos algorithm
for finding a few (chosen by the user) eigenvalues of a seminegative
(one eigenvalue 0, the rest negative) symmetric real matrix. There is a
nice parallel Fortran 90 implementation, TRLAN
(http://www.nersc.gov/~kewu/trlan.html), but I want to modify the code
and adopt it to global HPF, so something simpler (smaller) is preferred,
even if numerical robustness and speed are jeopardized. It is absolutely
essential that the routine accepts the matrix-vector product routine as
a user argument, so I can use my own multiplication routines.
Anybody aware of such code?
Thanks a lot,
Aleksandar
--
_____________________________________________
Aleksandar Donev
[log in to unmask]
http://www.pa.msu.edu/~donev/
Department of Physics and Astronomy
Michigan State University
East Lansing, MI 48825
(517) 432-6770
_____________________________________________
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
|