bugGNU Scientific Library - Bugs: bug #24704, error in gsl_cdf_chisq_Pinv for...

 
 

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

bug #24704: error in gsl_cdf_chisq_Pinv for very large degrees of freedom

Submitted by:  -Deleted Account- <bjg>
Submitted on:  Thu 30 Oct 2008 07:11:29 PM UTC  
 
Category: Runtime errorSeverity: 3 - Normal
Operating System: Status: Fixed
Assigned to: NoneOpen/Closed: Closed
Release: 1.11

Wed 03 Dec 2008 03:04:01 PM UTC, comment #1:

fixed by
d774e37a260b7187c578f3f304946b2dcf3e3ccf
e5e1db0a50b6ee0a0f6afec23528f3634360b2b6

possibly this will need to be revisited for larger orders of magnitude but this is enough to fix the problem reported with nu O(1e6).

-Deleted Account- <bjg>
Thu 30 Oct 2008 07:11:29 PM UTC, original submission:

Dear GSL,

I've encountered a potential bug in the routine "gsl_cdf_chisq_Pinv" defined in the file gsl-1.11/cdf/chisqinv.c.

When calling this routine with the degrees of freedom > 1263131, I get the following error message:

gsl: gamma_inc.c:99: ERROR: error
Default GSL error handler invoked.
Abort

For degrees of freedom <= 1263131, it seems to perform fine.

I'm running GSL version 1.11 and calling the routine as follows:
v = gsl_cdf_chisq_Pinv(0.05, nu);
where nu = the degrees of freedom.

Thanks,
John Halley Gotway
-unavailable-

-Deleted Account- <bjg>

 

Attached Files
file #16760:  chisq3.c added by bjg (508B - text/x-csrc)

 

Depends on the following items: None found

Items that depend on this one: None found

 

Carbon-Copy List
  • -unavailable- added by bjg
  •  

    Please enter the title of George Orwell's famous dystopian book (it's a date):

     

     

    Follow 4 latest changes.

    Date Changed By Updated Field Previous Value => Replaced By
    Wed 03 Dec 2008 03:04:01 PM UTCbjgStatusConfirmed=>Fixed
      Open/ClosedOpen=>Closed
    Thu 30 Oct 2008 07:11:29 PM UTCbjgAttached File-=>Added chisq3.c, #16760
      Carbon-Copy-=>Added -unavailable-

    Back to the top


    Powered by Savane 3.1-cleanup1