URL:
<http://savannah.gnu.org/bugs/?31362>
Summary: The Complete Elliptic Integrals
(gsl_sf_ellint_Ecomp and _Kcomp) Loop Forever with NaN Argument
Project: GNU Scientific Library
Submitted by: wmfarr
Submitted on: Mon 18 Oct 2010 01:34:26 PM GMT
Category: Runtime error
Severity: 3 - Normal
Operating System: Mac OS X 10.6.4
Status: None
Assigned to: None
Open/Closed: Open
Release: 1.14
Discussion Lock: Any
_______________________________________________________
Details:
Feeding a NaN to the complete elliptic integrals causes these functions to
loop forever. Compile and execute the following program:
-------------------
#include<gsl/gsl_sf.h>
int main() {
double nan = 0.0/0.0;
printf("Elliptic integral of nan = %g\n", gsl_sf_ellint_Ecomp(nan,
GSL_PREC_DOUBLE));
printf("Elliptic integral of nan = %g\n", gsl_sf_ellint_Kcomp(nan,
GSL_PREC_DOUBLE));
return 0;
}
--------------------
Nothing will ever be printed. I haven't checked any of the other elliptic
integral special functions. Probably these functions should return NaN with
NaN arguments, but they could also call gsl_error(...) I suppose. In any
case, the should not loop forever.
_______________________________________________________
Reply to this item at:
<http://savannah.gnu.org/bugs/?31362>
_______________________________________________
Message sent via/by Savannah
http://savannah.gnu.org/
_______________________________________________
Bug-gsl mailing list
[email protected]
http://lists.gnu.org/mailman/listinfo/bug-gsl
