Follow-up Comment #1, bug #39056 (group gsl): Wolpert's test case succeeds with the patch given in https://savannah.gnu.org/bugs/?21835.
With this patch Monajemi's test case still fails in GSL-2.8. Note that this
case may have been motivated by a and b both rational and differ by an
integer, as well as c in the neighborhood of a + b and x between 0 and 1. Such
cases typically fail.
For example, while the following two elegant results hold:
1- hyperg_2F1 (a, 1/2, a, 1/2) = sqrt(2)
2- hyperg_2F1 (a, -1/2, a, 1) = 0
they fail in GSL-2.8:
1- TEST_SF(s, gsl_sf_hyperg_2F1_e, (3, -1/2, 3, 1/2, &r), M_SQRT1_2,
TEST_TOL0, GSL_SUCCESS);
2- TEST_SF(s, gsl_sf_hyperg_2F1_e, (5, -1/2, 5, 1, &r), 0.0, TEST_TOL0,
GSL_SUCCESS);
_______________________________________________________
Reply to this item at:
<https://savannah.gnu.org/bugs/?39056>
_______________________________________________
Message sent via Savannah
https://savannah.gnu.org/
signature.asc
Description: PGP signature
