[Bug-gsl] seg fault from gsl_sf_legendre_sphPlm_e

[Kevin M. Huffenberger]

Sorry I was missing some information.

GSL version 1.13 installed from RPM from 
http://www.network-theory.co.uk/gsl/freedownloads.html

Compiler:
Reading specs from /usr/lib/gcc-lib/i386-redhat-linux/3.2.2/specs
Configured with: ../configure --prefix=/usr --mandir=/usr/share/man 
--infodir=/usr/share/info --enable-shared --enable-threads=posix 
--disable-checking --with-system-zlib --enable-__cxa_atexit 
--host=i386-redhat-linux
Thread model: posix
gcc version 3.2.2 20030222 (Red Hat Linux 3.2.2-5)
Used to compile: gcc -o test test.c -lm -lgsl -lgslcblas

Linux kernal 2.4.20-19.9 on an i686 (Compaq laptop)

After a bit more investigation, the routine appears to fail when x is 
close to 1 and the m gets large.  This does not appear to depend on l.  
For example,   
gsl_sf_legendre_sphPlm_e(140,131,0.99999,&y);  
gsl_sf_legendre_sphPlm_e(10000,131,0.99999,&y);
are okay, while
gsl_sf_legendre_sphPlm_e(140,132,0.99999,&y);
gsl_sf_legendre_sphPlm_e(10000,132,0.99999,&y);
cause a segmentation fault.

Thanks,

---Kevin