Re: [Bug-gsl] seg fault from gsl_sf_legendre_sphPlm_e

[Brian Gough]

Kevin M. Huffenberger writes:
 > In the following bit of code, the first call works, but the second call 
 > dies.
 > 
 > #include <stdlib.h>
 > #include <stdio.h>
 > #include <math.h>
 > #include <gsl/gsl_sf_legendre.h>
 > 
 > int main() {
 >   gsl_sf_result y;
 >   gsl_sf_legendre_sphPlm_e(140,135,1,&y);
 >   gsl_sf_legendre_sphPlm_e(140,135,0.99998689456491752,&y);
 >   return(0);
 > }
 > 

Thanks for the bug report.  Can you tell me what the correct numerical
value for the second case should be.

regards,

-- 
Brian Gough

Network Theory Ltd            
15 Royal Park                 
Bristol BS8 3AL               
United Kingdom                

Tel: +44 (0)117 3179309 
Fax: +44 (0)117 9048108              
Web: http://www.network-theory.co.uk/