bugGNU Scientific Library - Bugs: bug #39056, gsl_sf_hyperg_2F1_e fails for some...

 
 

You are not allowed to post comments on this tracker with your current authentication level.

bug #39056: gsl_sf_hyperg_2F1_e fails for some test cases

Submitter:  Patrick Alken <psa>
Submitted:  Thu 23 May 2013 08:26:48 PM UTC
   
 
Category:  Runtime error Severity:  4 - Important
Operating System:  all Status:  Confirmed
Assigned to:  None Open/Closed:  Open
Release: 

Thu 30 Jan 2025 09:46:20 AM UTC, comment #1: 

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);

Fermé le Lundi <fermelelundi>
Thu 23 May 2013 08:26:48 PM UTC, original submission:  

The following test cases are confirmed to fail for the gsl_sf_hyperg_2F1_e function. They have been commented out of specfunc/test.c to enable 'make check' to succeed.


  /* Test case from Hatef Monajemi <monajemi@stanford.edu> */

  TEST_SF(s, gsl_sf_hyperg_2F1_e, (3.5, -0.5, 5.0, 0.9, &r), 0.5923981284370653465208973272, TEST_TOL2, GSL_SUCCESS);

  /* Test case from Robert L Wolpert <Wolpert@stat.duke.edu> */

  TEST_SF(s, gsl_sf_hyperg_2F1_e, (-1.0, -10.0, 1.0, 0.5, &r), 6.0, TEST_TOL0, GSL_SUCCESS);

Patrick Alken <psa>
Group administrator

 

(Note: upload size limit is set to 16384 kB, after insertion of the required escape characters.)

No files currently attached

 

Depends on the following items: None found

Items that depend on this one: None found

 

Carbon-Copy List
  • -email is unavailable- added by fermelelundi (Posted a comment)
  • -email is unavailable- added by psa (Submitted the item)
  •  

    No changes have been made to this item

    Back to the top

    Powered by Savane 3.14-9aa3.
    Corresponding source code