bugGNU Scientific Library - Bugs: bug #24812, gsl_sf_hyperg_2F1(11, -1 ; 11/2;...

 
 
Show feedback again

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

bug #24812: gsl_sf_hyperg_2F1(11, -1 ; 11/2; 0.125) fails

Submitted by:  Brian Gough <bjg>
Submitted on:  Tue 11 Nov 2008 02:43:07 PM UTC  
 
Category: Runtime errorSeverity: 3 - Normal
Operating System: Status: Confirmed
Assigned to: NoneOpen/Closed: Open
Release: 1.11

Tue 11 Nov 2008 02:45:01 PM UTC, comment #1:

line 732:
if(GSL_MAX_DBL(fabs(a),1.0)fabs(bp)fabs(x) < 2.0*fabs(c)) {
line 740:
if(fabs(bpbpxx) < 0.001fabs(bp) && fabs(a) < 10.0) {

I think these should have a=>ap for consistency

Brian Gough <bjg>
Project Administrator
Tue 11 Nov 2008 02:43:07 PM UTC, original submission:

From: "Didier Pinchon" <pinchon.didier@orange.fr>
To: <bug-gsl@gnu.org>
Subject: [Bug-gsl] Bug in gsl_sf_hyperg_2F1 ??
Date: Sat, 8 Nov 2008 03:00:41 +0100

Hello,

I have tried to compute 2F1(11, -1 ; 11/2; 0.125) and I got an error message (below is my sample program).
However 2F1(-1, 11 ; 11/2; 0.125) provides the right result 0.75

Am I wrong somewhere ?
I did not find any indication in archives.

Thank you for your help.
All the best,
Didier

/* Compilation and execution:
$ gcc -o bug_gsl -I/usr/local/include/gsl bug_gsl.c -L/usr/local/lib -lgsl -lm

$ ./bug_gsl
gsl: hyperg_2F1.c:750: ERROR: error
Default GSL error handler invoked.
Abandon

*/

Brian Gough <bjg>
Project Administrator

 

Attached Files
file #16829:  hyperg2F1.c added by bjg (364B - application/octet-stream)

 

Depends on the following items: None found

Items that depend on this one: None found

 

Carbon-Copy List
  • -unavailable- added by bjg (Submitted the item)
  •  

    Please enter 421 here (basic anti-spam test):

     

     

    Follows 1 latest change.

    Date Changed By Updated Field Previous Value => Replaced By
    Tue 11 Nov 2008 02:43:07 PM UTCbjgAttached File-=>Added hyperg2F1.c, #16829
    Show feedback again

    Back to the top


    Powered by Savane 3.1-cleanup