aboutsummaryrefslogtreecommitdiff
diff options
context:
space:
mode:
authorAndreas Jaeger <aj@suse.de>2001-06-19 12:41:02 +0000
committerAndreas Jaeger <aj@suse.de>2001-06-19 12:41:02 +0000
commit20f421e1a1bdc4eb203198825120975a52dd1e30 (patch)
treeea998677a162a1ef2bdbe815047feccfcc5ce884 /sysdeps
parent128 bit implementation of tanh. (diff)
downloadglibc-20f421e1a1bdc4eb203198825120975a52dd1e30.tar.gz
glibc-20f421e1a1bdc4eb203198825120975a52dd1e30.tar.bz2
glibc-20f421e1a1bdc4eb203198825120975a52dd1e30.zip
Update
2001-06-19 Andreas Jaeger <aj@suse.de> * sysdeps/ieee754/ldbl-128/s_tanhl.c: New file. * math/libm-test.inc (tanh_test): Test for 2^-56. * sysdeps/ieee754/ldbl-96/s_tanhl.c (__tanhl): Make sure result equals argument when x < 2^-55. Patches by Stephen L. Moshier <moshier@na-net.ornl.gov>.
Diffstat (limited to 'sysdeps')
-rw-r--r--sysdeps/ieee754/ldbl-96/s_tanhl.c2
1 files changed, 1 insertions, 1 deletions
diff --git a/sysdeps/ieee754/ldbl-96/s_tanhl.c b/sysdeps/ieee754/ldbl-96/s_tanhl.c
index 5c0ac9a95d..45cea903e3 100644
--- a/sysdeps/ieee754/ldbl-96/s_tanhl.c
+++ b/sysdeps/ieee754/ldbl-96/s_tanhl.c
@@ -78,7 +78,7 @@ static long double one=1.0, two=2.0, tiny = 1.0e-4900L;
if ((ix|j0|j1) == 0)
return x; /* x == +- 0 */
if (ix<0x3fc8) /* |x|<2**-55 */
- return x*(one+x); /* tanh(small) = small */
+ return x*(one+tiny); /* tanh(small) = small */
if (ix>=0x3fff) { /* |x|>=1 */
t = __expm1l(two*fabsl(x));
z = one - two/(t+two);