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

[Patrick Alken]

Follow-up Comment #4, bug #39466 (project gsl):

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

(file #28555)
    _______________________________________________________

Additional Item Attachment:

File name: coupling_6_9_halfint_check_2.diff Size:2 KB


    _______________________________________________________

Reply to this item at:

  <http://savannah.gnu.org/bugs/?39466>

_______________________________________________
  Message sent via/by Savannah
  http://savannah.gnu.org/