[Getdp] Q: Lanczos error
Christophe Geuzaine
c.geuzaine at ulg.ac.be
Thu Jul 1 17:15:15 CEST 2004
Thanks a lot for the patch!
Best,
Christophe
ohyeahq wrote:
> Christophe,
>
> Thank you for your reply. I really like to use getdp for my problem (it's
> TOO nice to be left unused!), so I tried to make lanczos work. It seems
> ok,
> although I haven't done extensive checking against Numerical Recipes
> version. Here's my patch to Numeric/Lanczos.c and config.in. What I
> did are:
>
> (1) ran3() is replaced by GSL random number generator.
>
> (2) GSL doesn't have hqr() analogue, so I use a LAPACK routine for
> Hessenberg
> eigensolver. Checking for LAPACK availability is added to config.in.
> Relying on LAPACK is obviously minus in portability and performance (c to
> fortran interfacing), but LAPACK routines are likely to be useful in future
> extension of getdp (I hope).
>
--
Christophe Geuzaine
Applied and Computational Mathematics, Caltech
geuzaine at acm.caltech.edu - http://geuz.org