GNU Scientific Library - Bugs: bug #23101, gsl_cdf_chisq_Pinv inverse failed...
You are not allowed to post comments on this tracker with your current authentication level.
bug #23101: gsl_cdf_chisq_Pinv inverse failed to converge
Submitter: | -Deleted Account- <bjg> | ||
Submitted: | Tue 29 Apr 2008 01:36:20 PM UTC | ||
Category: | Accuracy problem | Severity: | 3 - Normal |
Operating System: | Status: | Fixed | |
Assigned to: | bjg | Open/Closed: | Closed |
Release: | 1.11 |
Tue 29 Apr 2008 01:39:36 PM UTC, comment #1: |
-Deleted Account- <bjg>![]() |
Tue 29 Apr 2008 01:36:20 PM UTC, original submission:
While testing my use of gsl_cdf_chisq_Pinv I ran into the following
gsl: gammainv.c:105: ERROR: inverse failed to converge
I looked into this and discovered the issue is in the routine
if (dP == 0.0 || n++ > 32)
to
if (dP == 0.0 || n++ > 40)
and it no longer fails. In fact, the problem is most likely in line
double x0 = (xg < -sqrt (a)) ? a : sqrt (a) * xg + a;
For nu = 1.5 and around P = 1.93238145206123590e-01, xg is very close
Thanks,
|
-Deleted Account- <bjg>![]() |
Depends on the following items: None found
Items that depend on this one: None found
fixed by 6d5d363d069aa108b351dff86872f7e10137f496