This is the mail archive of the
gsl-discuss@sources.redhat.com
mailing list for the GSL project.
GSL_Wigner_6J error
- From: Peter Roche <P dot J dot P dot Roche at damtp dot cam dot ac dot uk>
- To: gsl-discuss at sources dot redhat dot com
- Date: Wed, 13 Nov 2002 19:36:40 +0000 (GMT)
- Subject: GSL_Wigner_6J error
Dear All,
I hope this isn't a repeat of an already dealt with problem. I have been
trying to use the GSL Wigner 6j coefficients routines
(gsl_sf_coupling_6j(2*ja, 2*jb, 2*jc, 2*jd,2*je,2*jf)) and have found,
I think, that some of the returned values from the GSL
routine are not correct, even for small input values.
I have included a table showing, in the first column, the input
parameter values, the angular mometa (which are multiplied by 2 in my
actual program), the second column shows the 'correct' values, obtained
from another program/published tables, and the third column shows the
return values from the GSL routine.
Input Correct values GSL return
(011100) 0.5773502691896257 0.0000000000000000
(101101) 0.3333333333333335 0.0000000000000000
(121101) 0.3333333333333334 0.0000000000000000
(101101) 0.3333333333333335 0.0000000000000000
(022200) 0.4472135954999579 0.0000000000000000
(132201) 0.2581988897471612 0.0000000000000000
(242202) 0.2000000000000000 0.0000000000000000
(213302) 0.1690308509457034 0.0000000000000000
(233302) 0.1690308509457033 0.0000000000000000
(253302) 0.1690308509457033 0.0000000000000000
Has anyone encountered the same problems, or can explain the differences.
Cheers,
Peter Roche
*********************************************************
Peter Roche
Department of Applied Mathematics and Theoretical Physics,
University of Cambridge, Silver Street, Cambridge, CB3 9EW
email: p.j.p.roche@damtp.cam.ac.uk
*********************************************************