Re: [Bug-gsl] [bug #39466] errors in gsl_sf_coupling_6j_e and gsl_sf_coupling_9j_e

[Håkan Johansson]

Yes, this is a duplicate of #29606.  Sorry for not finding that.

Attached is a new patch (based on Alexey's suggestion) that now should 
apply cleanly to bzr trunk. The previous was against v1.15.  It also 
includes the tested case in #29606 in test_sf.c.

For completeness, while the case in #29606 is correct (6j to be 0 as 
half-integer with integer cannot couple to integer), I think there were 
a typo in the report: jc=21 should be jd=22, giving

abs(1-21.5)=20.5 <= 22 <= 22.5.

instead of

abs(1-21.5)=20.5 <= 21 <= 22.5.

Cheers,
Håkan


On Thu, 11 Jul 2013, Patrick Alken wrote:

> Follow-up Comment #3, bug #39466 (project gsl):
>
> This may be a repeat of bug #29606
>
>    _______________________________________________________
>
> Reply to this item at:
>
>  <http://savannah.gnu.org/bugs/?39466>
>
> _______________________________________________
>  Message sent via/by Savannah
>  http://savannah.gnu.org/

coupling_6_9_halfint_check_2.diff
Description: Text Data