[Top][All Lists]
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
Re: [Bug-gsl] Bug in gsl_sf_coupling_3j_e
From: |
Brian Gough |
Subject: |
Re: [Bug-gsl] Bug in gsl_sf_coupling_3j_e |
Date: |
Mon, 15 Aug 2011 16:39:29 +0100 |
User-agent: |
Wanderlust/2.15.6 (Almost Unreal) Emacs/23.3 Mule/6.0 (HANACHIRUSATO) |
At Fri, 12 Aug 2011 18:10:17 +0200,
Martin Landriau wrote:
> I found a problem with the function gsl_sf_coupling_3j_e located in the file
> specfunc/coupling.c : the tests should pick up that
>
> (l1 l2 l3) = 0
> (0 0 0 )
>
> if l1+l2+l3 is not an even integer and it doesn't. It instead attempts to
> evaluate it and sometimes returns a non-zero value, albeit one with an error
> larger than the value itself.
Thanks for the email, could you send a small example program showing
some cases that return a non-zero value so I can reproduce the problem.
--
Brian Gough