[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