[Bug-gsl] gsl_sf_coupling_6j_INCORRECT_e

[Eric Diffenderfer]

GSL version 1.3

I geuss this is a bug:
There is a code error in function
gsl_sf_coupling_6j_INCORRECT_e line 200
of coupling.c in the specfunc directory.

The mistake is on line 207:


if(   two_ja < 0 || two_jb < 0 || two_jc < 0
     || two_jd < 0 || two_je < 0 || two_je < 0
     ) {
    DOMAIN_ERROR(result);
  }

It checks for (two_je < 0) twice.  It should be:
        
if(   two_ja < 0 || two_jb < 0 || two_jc < 0
     || two_jd < 0 || two_je < 0 || two_jf < 0
     ) {
    DOMAIN_ERROR(result);
  }