diff options
Diffstat (limited to 'sysdeps/ieee754/flt-32/s_erff.c')
-rw-r--r-- | sysdeps/ieee754/flt-32/s_erff.c | 11 |
1 files changed, 7 insertions, 4 deletions
diff --git a/sysdeps/ieee754/flt-32/s_erff.c b/sysdeps/ieee754/flt-32/s_erff.c index c8b6287503..693934d59a 100644 --- a/sysdeps/ieee754/flt-32/s_erff.c +++ b/sysdeps/ieee754/flt-32/s_erff.c @@ -20,7 +20,10 @@ static char rcsid[] = "$NetBSD: s_erff.c,v 1.4 1995/05/10 20:47:07 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-float.h> #include <fix-int-fp-convert-zero.h> static const float @@ -104,7 +107,7 @@ float __erff(float x) GET_FLOAT_WORD(hx,x); ix = hx&0x7fffffff; if(ix>=0x7f800000) { /* erf(nan)=nan */ - i = ((u_int32_t)hx>>31)<<1; + i = ((uint32_t)hx>>31)<<1; return (float)(1-i)+one/x; /* erf(+-inf)=+-1 */ } @@ -152,7 +155,7 @@ float __erff(float x) r = __ieee754_expf(-z*z-(float)0.5625)*__ieee754_expf((z-x)*(z+x)+R/S); if(hx>=0) return one-r/x; else return r/x-one; } -weak_alias (__erff, erff) +libm_alias_float (__erf, erf) float __erfcf(float x) { @@ -162,7 +165,7 @@ float __erfcf(float x) ix = hx&0x7fffffff; if(ix>=0x7f800000) { /* erfc(nan)=nan */ /* erfc(+-inf)=0,2 */ - float ret = (float)(((u_int32_t)hx>>31)<<1)+one/x; + float ret = (float)(((uint32_t)hx>>31)<<1)+one/x; if (FIX_INT_FP_CONVERT_ZERO && ret == 0.0f) return 0.0f; return ret; @@ -227,4 +230,4 @@ float __erfcf(float x) return two-tiny; } } -weak_alias (__erfcf, erfcf) +libm_alias_float (__erfc, erfc) |