URL: <https://savannah.gnu.org/bugs/?58063>
Summary: Dropped NaNs in multiple cases Project: GNU Scientific Library Submitted by: jacksonvanover Submitted on: Sat 28 Mar 2020 04:58:51 PM UTC Category: Accuracy problem Severity: 3 - Normal Operating System: Ubuntu 18.04 Status: None Assigned to: None Open/Closed: Open Release: 2.6 Discussion Lock: Any _______________________________________________________ Details: Using version 2.6 of GSL from ftp://ftp.gnu.org/gnu/gsl/gsl-2.6.tar.gz OS is Ubuntu 18.04.2 LTS Hardware is a Dell workstation with a 3.60 GHz Intel i7-4790 and 32 GB of RAM Compiler is gcc 7.4.0 with the -w option Descriptions: when the second argument to gsl_sf_bessel_zero_Jnu is a NaN, that NaN is dropped and the function returns a zero without complaining. The same is true of an inf input. When the first argument is zero and second argument is a NaN, gsl_sf_hermite returns 1. When the first argument is nonzero, the error handler is invoked as would be expected. When the first argument is zero and third argument is a NaN, gsl_sf_gegenpoly_n returns 1. When the first argument is nonzero, the error handler is invoked as would be expected. #include <gsl/gsl_sf.h> #include <stdio.h> #include <gsl/gsl_math.h> int main (void){ double out; out = gsl_sf_bessel_zero_Jnu(1, GSL_NAN); printf("%f\n", out); printf("%.*e\n", out); out = gsl_sf_hermite(0, GSL_NAN); printf("%f\n", out); printf("%.21e\n", out); out = gsl_sf_gegenpoly_n(0, 2.5, GSL_NAN); printf("%f\n", out); printf("%.21e\n", out); return 0; } _______________________________________________________ Reply to this item at: <https://savannah.gnu.org/bugs/?58063> _______________________________________________ Message sent via Savannah https://savannah.gnu.org/