diff options
Diffstat (limited to 'sysdeps/ieee754/dbl-64/s_erf.c')
-rw-r--r-- | sysdeps/ieee754/dbl-64/s_erf.c | 19 |
1 files changed, 7 insertions, 12 deletions
diff --git a/sysdeps/ieee754/dbl-64/s_erf.c b/sysdeps/ieee754/dbl-64/s_erf.c index b4975a8af8..5f820604f8 100644 --- a/sysdeps/ieee754/dbl-64/s_erf.c +++ b/sysdeps/ieee754/dbl-64/s_erf.c @@ -115,7 +115,10 @@ static char rcsid[] = "$NetBSD: s_erf.c,v 1.8 1995/05/10 20:47:05 jtc Exp $"; #include <errno.h> #include <float.h> #include <math.h> +#include <math-narrow-eval.h> #include <math_private.h> +#include <math-underflow.h> +#include <libm-alias-double.h> #include <fix-int-fp-convert-zero.h> static const double @@ -201,7 +204,7 @@ __erf (double x) ix = hx & 0x7fffffff; if (ix >= 0x7ff00000) /* erf(nan)=nan */ { - i = ((u_int32_t) hx >> 31) << 1; + i = ((uint32_t) hx >> 31) << 1; return (double) (1 - i) + one / x; /* erf(+-inf)=+-1 */ } @@ -294,11 +297,7 @@ __erf (double x) else return r / x - one; } -weak_alias (__erf, erf) -#ifdef NO_LONG_DOUBLE -strong_alias (__erf, __erfl) -weak_alias (__erf, erfl) -#endif +libm_alias_double (__erf, erf) double __erfc (double x) @@ -309,7 +308,7 @@ __erfc (double x) ix = hx & 0x7fffffff; if (ix >= 0x7ff00000) /* erfc(nan)=nan */ { /* erfc(+-inf)=0,2 */ - double ret = (double) (((u_int32_t) hx >> 31) << 1) + one / x; + double ret = (double) (((uint32_t) hx >> 31) << 1) + one / x; if (FIX_INT_FP_CONVERT_ZERO && ret == 0.0) return 0.0; return ret; @@ -421,8 +420,4 @@ __erfc (double x) return two - tiny; } } -weak_alias (__erfc, erfc) -#ifdef NO_LONG_DOUBLE -strong_alias (__erfc, __erfcl) -weak_alias (__erfc, erfcl) -#endif +libm_alias_double (__erfc, erfc) |