GNU Scientific Library - Bugs: bug #24146, bug in gsl_sf_ellint_Kcomp_e for k...
You are not allowed to post comments on this tracker with your current authentication level.
bug #24146: bug in gsl_sf_ellint_Kcomp_e for k near 1
Submitter: | -Deleted Account- <bjg> | ||
Submitted: | Tue 26 Aug 2008 04:36:55 PM UTC | ||
Category: | Accuracy problem | Severity: | 3 - Normal |
Operating System: | Status: | Fixed | |
Assigned to: | bjg | Open/Closed: | Closed |
Release: | 1.11 |
Tue 26 Aug 2008 04:39:58 PM UTC, comment #1: |
-Deleted Account- <bjg>![]() |
Tue 26 Aug 2008 04:36:55 PM UTC, original submission:
From: Thies Heidecke <theidecke@astrophysik.uni-kiel.de>
i'm using version 1.11 of the library.
gsl_sf_ellint_Kcomp_e(x, GSL_PREC_DOUBLE);
gives wrong results if x approaches the singularity at 1.0 as can be
#include <stdio.h>
int main (int argc, char const* argv[])
for (i=1;i<=16;++i) {
return 0;
the value should approach infinity but stays nearly constant for
looking at ellint.c the error is a simple typo in the part of the code
/* [Carlson, Numer. Math. 33 (1979) 1, (4.5)] */
line 489 reads:
but it should be
after recompiling the function now gives reasonable values.
Regards,
|
-Deleted Account- <bjg>![]() |
No files currently attached
Depends on the following items: None found
Items that depend on this one: None found
fixed by 163bd7269d3c74af6b0a81ac82d5396a2afd3b3e