diff options
| author | Rich Felker <dalias@aerifal.cx> | 2012-03-19 22:07:43 -0400 |
|---|---|---|
| committer | Rich Felker <dalias@aerifal.cx> | 2012-03-19 22:07:43 -0400 |
| commit | 97721a5508415a2f10eb068e022093811c9ff8be (patch) | |
| tree | 88e9ce153895ad949576fa7ce1eeee4b02286479 /src/math/sinhl.c | |
| parent | acb744921b73f5a73803e533e5e4a4896d164a26 (diff) | |
| parent | 0cbb65479147ecdaa664e88cc2a5a925f3de502f (diff) | |
| download | musl-97721a5508415a2f10eb068e022093811c9ff8be.tar.gz | |
Merge remote branch 'nsz/master'
Diffstat (limited to 'src/math/sinhl.c')
| -rw-r--r-- | src/math/sinhl.c | 8 |
1 files changed, 4 insertions, 4 deletions
diff --git a/src/math/sinhl.c b/src/math/sinhl.c index 2252dec9..8a54677e 100644 --- a/src/math/sinhl.c +++ b/src/math/sinhl.c @@ -35,7 +35,7 @@ long double sinhl(long double x) return sinh(x); } #elif LDBL_MANT_DIG == 64 && LDBL_MAX_EXP == 16384 -static const long double one = 1.0, huge = 1.0e4931L; +static const long double huge = 1.0e4931L; long double sinhl(long double x) { @@ -55,12 +55,12 @@ long double sinhl(long double x) /* |x| in [0,25], return sign(x)*0.5*(E+E/(E+1))) */ if (ix < 0x4003 || (ix == 0x4003 && i0 <= 0xc8000000)) { /* |x| < 25 */ if (ix < 0x3fdf) /* |x|<2**-32 */ - if (huge + x > one) + if (huge + x > 1.0) return x;/* sinh(tiny) = tiny with inexact */ t = expm1l(fabsl(x)); if (ix < 0x3fff) - return h*(2.0*t - t*t/(t + one)); - return h*(t + t/(t + one)); + return h*(2.0*t - t*t/(t + 1.0)); + return h*(t + t/(t + 1.0)); } /* |x| in [25, log(maxdouble)] return 0.5*exp(|x|) */ |