[Bug-gsl] Problem with gsl_sf_legendre_sphPlm

[James Casey]

Hello,

I am running GSL version 1.15.  I have originally ran into the problem
with pygsl and verified that the problem is actually in GSL by running
the C libraries.  When running:

gsl_sf_legendre_sphPlm(l,m,cos(t))

I have problems for large l, and small t.  Specifically, with l = 192,
m = 1, and t = 0.01,  the values of the spherical harmonics are
completely wrong.  The complex component is several orders of
magnitude too large and the real part is greater than one, which
should not happen for a normalized function.  I haven't tried to look
into the code myself, yet, but I wanted to know if anyone else had
reported this or if anyone else was experiencing similar problems.

Cheers,
James Casey

-- 
James Casey, Graduate Research Assistant
Center for Relativistic Astrophysics
School of Physics - Georgia Institute of Technology
837 State St. Atlanta GA, 30332-0430