diff options
Diffstat (limited to 'sysdeps/libm-ieee754/s_asinhf.c')
| -rw-r--r-- | sysdeps/libm-ieee754/s_asinhf.c | 20 |
1 files changed, 10 insertions, 10 deletions
diff --git a/sysdeps/libm-ieee754/s_asinhf.c b/sysdeps/libm-ieee754/s_asinhf.c index d57508de9f..fac256d37a 100644 --- a/sysdeps/libm-ieee754/s_asinhf.c +++ b/sysdeps/libm-ieee754/s_asinhf.c @@ -8,7 +8,7 @@ * * Developed at SunPro, a Sun Microsystems, Inc. business. * Permission to use, copy, modify, and distribute this - * software is freely granted, provided that this notice + * software is freely granted, provided that this notice * is preserved. * ==================================================== */ @@ -21,13 +21,13 @@ static char rcsid[] = "$NetBSD: s_asinhf.c,v 1.5 1995/05/12 04:57:39 jtc Exp $"; #include "math_private.h" #ifdef __STDC__ -static const float +static const float #else -static float +static float #endif one = 1.0000000000e+00, /* 0x3F800000 */ ln2 = 6.9314718246e-01, /* 0x3f317218 */ -huge= 1.0000000000e+30; +huge= 1.0000000000e+30; #ifdef __STDC__ float __asinhf(float x) @@ -35,21 +35,21 @@ huge= 1.0000000000e+30; float __asinhf(x) float x; #endif -{ +{ float t,w; int32_t hx,ix; GET_FLOAT_WORD(hx,x); ix = hx&0x7fffffff; if(ix>=0x7f800000) return x+x; /* x is inf or NaN */ - if(ix< 0x31800000) { /* |x|<2**-28 */ + if(ix< 0x38000000) { /* |x|<2**-14 */ if(huge+x>one) return x; /* return x inexact except 0 */ - } - if(ix>0x4d800000) { /* |x| > 2**28 */ + } + if(ix>0x47000000) { /* |x| > 2**14 */ w = __ieee754_logf(fabsf(x))+ln2; - } else if (ix>0x40000000) { /* 2**28 > |x| > 2.0 */ + } else if (ix>0x40000000) { /* 2**14 > |x| > 2.0 */ t = fabsf(x); w = __ieee754_logf((float)2.0*t+one/(__ieee754_sqrtf(x*x+one)+t)); - } else { /* 2.0 > |x| > 2**-28 */ + } else { /* 2.0 > |x| > 2**-14 */ t = x*x; w =__log1pf(fabsf(x)+t/(one+__ieee754_sqrtf(one+t))); } |