This is the mail archive of the gsl-discuss@sources.redhat.com mailing list for the GSL project.


Index Nav: [Date Index] [Subject Index] [Author Index] [Thread Index]
Message Nav: [Date Prev] [Date Next] [Thread Prev] [Thread Next]
Other format: [Raw text]

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                                                |               
|                                                               |               
+---------------------------------------------------------------+


Index Nav: [Date Index] [Subject Index] [Author Index] [Thread Index]
Message Nav: [Date Prev] [Date Next] [Thread Prev] [Thread Next]