Fri 27 Aug 2010 08:42:16 AM UTC, comment #2:
From: Andrew Benson <abenson@caltech.edu>
To: -unavailable-
Subject: [Bug-gsl] Bug in gsl_sf_hyperg_2F1
Date: Tue, 24 Aug 2010 21:00:48 -0700
[1 <text/plain; us-ascii (7bit)>]
Under GSL 1.14 (compiled with either gcc version 4.4.4 or gcc version 4.1.2)
the gsl_sf_hyperg_2F1 function crashes with the following arguments:
gsl_sf_hyperg_2F1(-10.34,2.05,3.05,0.1725)
The expected answer (computed using Maple for example) is 0.3104735522
I've attached an example case:
$ g++ test.cpp -o test.x -lgsl -lgslcblas
$ test.x
gsl: hyperg_2F1.c:755: ERROR: error
Default GSL error handler invoked.
Abort
The crash seems to occur only when the first argument is less than or equal to
-10. For example,
gsl_sf_hyperg_2F1(-10.0,2.05,3.05,0.1725)
crashes, but
gsl_sf_hyperg_2F1(-9.99999999999,2.05,3.05,0.1725)
evaluates correctly to 0.321419346301974773.
|