[Top][All Lists]
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
[Bug-gsl] Wigner 3j coupling coefficients return zero for high angular m
From: |
Jakub Benda |
Subject: |
[Bug-gsl] Wigner 3j coupling coefficients return zero for high angular momenta |
Date: |
Mon, 23 Jul 2012 21:36:43 +0200 |
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
[Prev in Thread] |
Current Thread |
[Next in Thread] |
- [Bug-gsl] Wigner 3j coupling coefficients return zero for high angular momenta,
Jakub Benda <=