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

[Patrick Alken]

Hi, I've applied your patch to the bzr, thanks a lot for your work on 
this (you and Alexey). I'll close the two bug
reports but open a new one with Alexey's comment that there exists a 
better algorithm to compute these

Patrick

On 07/12/2013 03:22 AM, Håkan Johansson wrote:
> 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/
> >