This is the mail archive of the
gsl-discuss@sources.redhat.com
mailing list for the GSL project.
gsl_sf_coupling_3j
- From: Filip Floegel <floegel at itp dot uni-hannover dot de>
- To: gsl-discuss at sources dot redhat dot com
- Date: Wed, 2 Oct 2002 17:45:52 +0200
- Subject: gsl_sf_coupling_3j
hello gsl community,
gsl_sf_coupling_3j doesn't work properly for large
angular mementa like j1=j2=j3=40 m1=m2=m3=0 !!!
gsl_sf_coupling_3j_e(80,80,80,0,0,0,&result)
which correspond to j1=j2=j3=40 m1=m2=m3=0
gives me:
result.val = 0.0625
result.err = 0.172766
mathematica gives me:
N[ThreeJSymbol[{40,0},{40,0},{40,0}]] = 0.0149685
my guess is that gamma.c which is used for calculating fractional
numbers probably diverges for large integers.
cheers,
ff
--
+---------------------------------------------------------------+
| dipl.phys. Filip Floegel e-mail:floegel@itp.uni-hannover.de |
| ITP/Universitaet Hannover phone: 0511/762-3381 |
| Appelstr.2 fax: 0511/762-4945 |
| 30167 Hannover |
| |
+---------------------------------------------------------------+