GNU Scientific Library - Bugs: bug #28267, poor convergence region for...
You are not allowed to post comments on this tracker with your current authentication level.
bug #28267: poor convergence region for gsl_sf_hyperg_1F1
Submitter: | -Deleted Account- <bjg> | ||
Submitted: | Fri 11 Dec 2009 08:35:17 PM UTC | ||
Category: | Accuracy problem | Severity: | 3 - Normal |
Operating System: | Status: | Postponed | |
Assigned to: | None | Open/Closed: | Open |
Release: | 1.13 |
Thu 25 Feb 2010 03:42:25 PM UTC, comment #2: |
-Deleted Account- <bjg> |
Fri 11 Dec 2009 08:36:23 PM UTC, comment #1: It's inherently difficult to compute the value in this region either way as there is massive cancellation in both the series and the Kummer transformed series. |
-Deleted Account- <bjg> |
Fri 11 Dec 2009 08:35:17 PM UTC, original submission:
The magnitude of the error is greater than the value itself for
gsl_sf_hyperg1F1(-3.78e+01, 2, 1.035e+02) => -7.00055e+18 +/- 3.77654e+19
-----------------------------------------------
Hi, guys
I experienced bugs with gsl_sf_hyperg_1F1. The version is gsl_1.12, but
#include<stdlib.h>
int main(int argc, char **argv)
gsl_sf_result r;
gsl_sf_hyperg_1F1_e(1.0-v0,2,2.0*Ri/v0,&r);
}
The output of above code shows an extremely large error. by increasing
Thank you very much.
Best
Weibin
|
-Deleted Account- <bjg> |
Depends on the following items: None found
Items that depend on this one: None found
I cannot find any algorithm which handles this case.