bugGNU Scientific Library - Bugs: bug #21838, gsl_cdf_fdist_P() nan return value

 
 

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

bug #21838: gsl_cdf_fdist_P() nan return value

Submitted by:  -Deleted Account- <bjg>
Submitted on:  Tue 18 Dec 2007 06:54:06 PM UTC  
 
Category: Runtime errorSeverity: 3 - Normal
Operating System: Status: Fixed
Assigned to: -Deleted Account- <bjg>Open/Closed: Closed
Release: 1.8

Wed 20 Feb 2008 05:57:52 PM UTC, comment #4:

added an asymptotic case
ae1c59b7e1c37f635f2352abd2ab23349981dfe8

-Deleted Account- <bjg>In charge of this item.
Sun 17 Feb 2008 01:49:05 PM UTC, comment #3:

no that doesn't fix it!

-Deleted Account- <bjg>In charge of this item.
Sun 17 Feb 2008 12:37:27 PM UTC, comment #2:

fixed by commit cb8cb7ea464cf5b6e6dfb30006db4b25846089cb

-Deleted Account- <bjg>In charge of this item.
Sun 10 Feb 2008 10:28:06 AM UTC, comment #1:

added some test cases for this

-Deleted Account- <bjg>In charge of this item.
Tue 18 Dec 2007 06:54:06 PM UTC, original submission:

Added test case
TEST (gsl_cdf_fdist_P, (3.479082213465832574, 1, 4040712), 0.93785072763723411967, TEST_TOL6);

From: "Daniel Neilson" <ddneilson@gmail.com>
To: -unavailable-
Subject: [Bug-gsl] gsl_cdf_fdist_P() nan return value
Date: Thu, 20 Sep 2007 17:05:23 -0600

Hi there,
I just spent the entire day tracking down the source of a nan in some code
of mine, and tracked it back to gsl_cdf_fdist_P(x, nu1, nu2). When called
with:

double x = 0x1.bd5290e3772c6p+1 (printf("%a", x) format to get the full
binary representation as hex)
int nu1 = 1
int nu2 = 4040712

this function returns nan. Furthermore, for Y \in {100,10000} calling:

gsl_cdf_fdist_P( (double)((unsigned int)x * Y + 0.5) / Y, nu1, nu2)

also returns nan.

However, for Y \in {1000, 100000} it returns a proper non-nan value.

This is all with the version of GSL in the Fedora 7 packages gsl-1.8-3.fc7&
gsl-devel-1.8-3.fc7.

-Deleted Account- <bjg>In charge of this item.

 

No files currently attached

 

Depends on the following items: None found

Items that depend on this one: None found

 

CC list is empty

 

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

 

 

Follow 9 latest changes.

Date Changed By Updated Field Previous Value => Replaced By
Wed 20 Feb 2008 05:57:52 PM UTCbjgStatusConfirmed=>Fixed
  Open/ClosedOpen=>Closed
Sun 17 Feb 2008 01:49:05 PM UTCbjgStatusFixed=>Confirmed
  Open/ClosedClosed=>Open
Sun 17 Feb 2008 12:37:27 PM UTCbjgStatusConfirmed=>Fixed
  Open/ClosedOpen=>Closed
Sun 10 Feb 2008 10:28:06 AM UTCbjgStatusNone=>Confirmed
  Assigned toNone=>bjg
  Summarygsl_cdf_fdist_P() nan return value =>gsl_cdf_fdist_P() nan return value

Back to the top


Powered by Savane 3.1-cleanup1