[Bug-gsl] Wigner 3j coupling coefficients return zero for high angular momenta

[Jakub Benda]

GSL Version: 1.15
OS: openSUSE 12.1, kernel 3.1.10-1.16-desktop, x86-64
Compiler: both GCC 4.6.2 and ICC 12.1

Wigner 3j coupling coefficients return zero for high angular momenta.

Sample code

---------------------------------------------------------------
#include <stdio.h>
#include <gsl/gsl_sf.h>
int main (*void*)
{

printf("%g\n", gsl_sf_coupling_3j(2*1,2*254,2*255,0,0,0));
printf("%g\n", gsl_sf_coupling_3j(2*1,2*255,2*256,0,0,0));
printf("%g\n", gsl_sf_coupling_3j(2*1,2*256,2*257,0,0,0));
printf("%g\n", gsl_sf_coupling_3j(2*1,2*257,2*258,0,0,0));
printf("%g\n", gsl_sf_coupling_3j(2*1,2*258,2*259,0,0,0));

return 0;

}
---------------------------------------------------------------

will output

-0.0313113
0.0312501
0
0
0

whereas the true values are (as computed by Maxima 5.26)

-0.0313112747451549
0.03125005960481531
-0.03118920210428185
0.03112869877462236
-0.03106854619382862

Thank you
Jakub