diff options
| author | Roland McGrath <roland@gnu.org> | 2006-01-31 18:56:42 +0000 |
|---|---|---|
| committer | Roland McGrath <roland@gnu.org> | 2006-01-31 18:56:42 +0000 |
| commit | 16f0ecedb5726fa46fac6a73b633aabcf469962d (patch) | |
| tree | b58989963f24a39fb71398dbf90c8d6f67f9f421 /sysdeps/ieee754/ldbl-128ibm/k_sincosl.c | |
| parent | 27d7e3b54c07b8eff35d2ff07ff9eb8ea03c853e (diff) | |
| download | glibc-16f0ecedb5726fa46fac6a73b633aabcf469962d.zip glibc-16f0ecedb5726fa46fac6a73b633aabcf469962d.tar.gz glibc-16f0ecedb5726fa46fac6a73b633aabcf469962d.tar.bz2 | |
* sysdeps/ieee754/ldbl-128ibm/k_cosl.c (__kernel_cosl): Correct index
for __sincosl_table.
* sysdeps/ieee754/ldbl-128ibm/k_sincosl.c (__kernel_sincosl): Likewise.
* sysdeps/ieee754/ldbl-128ibm/k_sinl.c (__kernel_sinl): Likewise.
* sysdeps/ieee754/ldbl-128ibm/s_ceill.c: Correct sign of 0.0.
* sysdeps/ieee754/ldbl-128ibm/s_floorl.c: Likewise.
* sysdeps/ieee754/ldbl-128ibm/s_cprojl.c: New file.
* sysdeps/ieee754/ldbl-128ibm/s_ctanhl.c: New file.
* sysdeps/ieee754/ldbl-128ibm/s_ctanl.c: New file.
Diffstat (limited to 'sysdeps/ieee754/ldbl-128ibm/k_sincosl.c')
| -rw-r--r-- | sysdeps/ieee754/ldbl-128ibm/k_sincosl.c | 21 |
1 files changed, 20 insertions, 1 deletions
diff --git a/sysdeps/ieee754/ldbl-128ibm/k_sincosl.c b/sysdeps/ieee754/ldbl-128ibm/k_sincosl.c index 1cea6fe..cd2ce7a 100644 --- a/sysdeps/ieee754/ldbl-128ibm/k_sincosl.c +++ b/sysdeps/ieee754/ldbl-128ibm/k_sincosl.c @@ -131,6 +131,25 @@ __kernel_sincosl(long double x, long double y, long double *sinx, long double *c Chebyshev polynomial of degree 10(11) and compute sinl(h+l) = sinl(h)cosl(l) + cosl(h)sinl(l) and cosl(h+l) = cosl(h)cosl(l) - sinl(h)sinl(l). */ + int six = tix; + tix = ((six - 0x3ff00000) >> 4) + 0x3fff0000; + index = 0x3ffe - (tix >> 16); + hix = (tix + (0x200 << index)) & (0xfffffc00 << index); + x = fabsl (x); + switch (index) + { + case 0: index = ((45 << 10) + hix - 0x3ffe0000) >> 8; break; + case 1: index = ((13 << 11) + hix - 0x3ffd0000) >> 9; break; + default: + case 2: index = (hix - 0x3ffc3000) >> 10; break; + } + hix = (hix << 4) & 0x3fffffff; +/* + The following should work for double but generates the wrong index. + For now the code above converts double to ieee extended to compute + the index back to double for the h value. + + index = 0x3fe - (tix >> 20); hix = (tix + (0x2000 << index)) & (0xffffc000 << index); x = fabsl (x); @@ -141,7 +160,7 @@ __kernel_sincosl(long double x, long double y, long double *sinx, long double *c default: case 2: index = (hix - 0x3fc30000) >> 14; break; } - +*/ SET_LDOUBLE_WORDS64(h, ((u_int64_t)hix) << 32, 0); if (iy) l = y - (h - x); |