aboutsummaryrefslogtreecommitdiff
path: root/sysdeps/ieee754/ldbl-128/e_j0l.c
diff options
context:
space:
mode:
authorAndreas Jaeger <aj@suse.de>2001-08-27 13:55:13 +0000
committerAndreas Jaeger <aj@suse.de>2001-08-27 13:55:13 +0000
commitfa9ced58f6e12a594e5250df52875cf1d67a50f0 (patch)
treee3ac76bb393e0e98fdedf9bbb55485d67983dcf5 /sysdeps/ieee754/ldbl-128/e_j0l.c
parent82ad8034aa1f812bf221ab0d2fc2fc55e2d44654 (diff)
downloadglibc-fa9ced58f6e12a594e5250df52875cf1d67a50f0.tar
glibc-fa9ced58f6e12a594e5250df52875cf1d67a50f0.tar.gz
glibc-fa9ced58f6e12a594e5250df52875cf1d67a50f0.tar.bz2
glibc-fa9ced58f6e12a594e5250df52875cf1d67a50f0.zip
* sysdeps/i386/fpu/libm-test-ulps (float): Add ulps for new tests.
* math/libm-test.inc (j0_test): Add extra tests. * sysdeps/ieee754/ldbl-128/e_j0l.c (__ieee754_j0l): Use |x|, not x, in Hankel approximation. Patches by Stephen L Moshier.
Diffstat (limited to 'sysdeps/ieee754/ldbl-128/e_j0l.c')
-rw-r--r--sysdeps/ieee754/ldbl-128/e_j0l.c8
1 files changed, 4 insertions, 4 deletions
diff --git a/sysdeps/ieee754/ldbl-128/e_j0l.c b/sysdeps/ieee754/ldbl-128/e_j0l.c
index 22a43848b2..e0a61f0658 100644
--- a/sysdeps/ieee754/ldbl-128/e_j0l.c
+++ b/sysdeps/ieee754/ldbl-128/e_j0l.c
@@ -754,16 +754,16 @@ __ieee754_j0l (long double x)
= 1/sqrt(2) * (sin(x) - cos(x))
sin(x) +- cos(x) = -cos(2x)/(sin(x) -+ cos(x))
cf. Fdlibm. */
- c = cosl (x);
- s = sinl (x);
+ c = cosl (xx);
+ s = sinl (xx);
ss = s - c;
cc = s + c;
- z = -cosl (x + x);
+ z = -cosl (xx + xx);
if ((s * c) < 0)
cc = z / ss;
else
ss = z / cc;
- z = ONEOSQPI * (p * cc - q * ss) / sqrtl (x);
+ z = ONEOSQPI * (p * cc - q * ss) / sqrtl (xx);
return z;
}