From 855a7fdbfa6c48d987203c09a9416b4d42a6f959 Mon Sep 17 00:00:00 2001 From: Andreas Jaeger Date: Tue, 29 May 2001 15:41:11 +0000 Subject: (__ieee754_sinhl): sinhl(x) = x when x < 2^-32. --- sysdeps/ieee754/ldbl-96/e_sinhl.c | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/sysdeps/ieee754/ldbl-96/e_sinhl.c b/sysdeps/ieee754/ldbl-96/e_sinhl.c index 4f9cfe2c38..646d4fde82 100644 --- a/sysdeps/ieee754/ldbl-96/e_sinhl.c +++ b/sysdeps/ieee754/ldbl-96/e_sinhl.c @@ -66,7 +66,7 @@ static long double one = 1.0, shuge = 1.0e4931L; if (jx & 0x8000) h = -h; /* |x| in [0,25], return sign(x)*0.5*(E+E/(E+1))) */ if (ix < 0x4003 || (ix == 0x4003 && i0 <= 0xc8000000)) { /* |x|<25 */ - if (ix<0x3fe3) /* |x|<2**-28 */ + if (ix<0x3fdf) /* |x|<2**-32 */ if(shuge+x>one) return x;/* sinh(tiny) = tiny with inexact */ t = __expm1l(fabsl(x)); if(ix<0x3fff) return h*(2.0*t-t*t/(t+one)); -- cgit 1.4.1