[Bug-gsl] [bug #31362] The Complete Elliptic Integrals (gsl_sf_ellint_Ecomp and _Kcomp) Loop Forever with NaN Argument

[Will M. Farr]

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/