diff options
author | Anssi Hannula <anssi.hannula@bitwise.fi> | 2020-01-27 12:45:11 +0200 |
---|---|---|
committer | Siddhesh Poyarekar <siddhesh@sourceware.org> | 2020-12-18 12:10:31 +0530 |
commit | 69a7ca7705626177e09d2fec6cf571f2015b76fc (patch) | |
tree | e37bd9f3764b8e028e8b196e6425b602254aa336 /sysdeps/ieee754 | |
parent | f67f9c9af228f6b84579cb8c86312d3a7a206a55 (diff) | |
download | glibc-69a7ca7705626177e09d2fec6cf571f2015b76fc.tar.gz glibc-69a7ca7705626177e09d2fec6cf571f2015b76fc.tar.xz glibc-69a7ca7705626177e09d2fec6cf571f2015b76fc.zip |
ieee754: Remove unused __sin32 and __cos32
The __sin32 and __cos32 functions were only used in the now removed slow path of asin and acos.
Diffstat (limited to 'sysdeps/ieee754')
-rw-r--r-- | sysdeps/ieee754/dbl-64/sincos32.c | 62 |
1 files changed, 0 insertions, 62 deletions
diff --git a/sysdeps/ieee754/dbl-64/sincos32.c b/sysdeps/ieee754/dbl-64/sincos32.c index c00e5d1a0f..a28932dffe 100644 --- a/sysdeps/ieee754/dbl-64/sincos32.c +++ b/sysdeps/ieee754/dbl-64/sincos32.c @@ -127,68 +127,6 @@ __c32 (mp_no *x, mp_no *y, mp_no *z, int p) __cpy (&s, z, p); } -/* Receive double x and two double results of sin(x) and return result which is - more accurate, computing sin(x) with multi precision routine c32. */ -double -SECTION -__sin32 (double x, double res, double res1) -{ - int p; - mp_no a, b, c; - p = 32; - __dbl_mp (res, &a, p); - __dbl_mp (0.5 * (res1 - res), &b, p); - __add (&a, &b, &c, p); - if (x > 0.8) - { - __sub (&hp, &c, &a, p); - __c32 (&a, &b, &c, p); - } - else - __c32 (&c, &a, &b, p); /* b=sin(0.5*(res+res1)) */ - __dbl_mp (x, &c, p); /* c = x */ - __sub (&b, &c, &a, p); - /* if a > 0 return min (res, res1), otherwise return max (res, res1). */ - if ((a.d[0] > 0 && res >= res1) || (a.d[0] <= 0 && res <= res1)) - res = res1; - LIBC_PROBE (slowasin, 2, &res, &x); - return res; -} - -/* Receive double x and two double results of cos(x) and return result which is - more accurate, computing cos(x) with multi precision routine c32. */ -double -SECTION -__cos32 (double x, double res, double res1) -{ - int p; - mp_no a, b, c; - p = 32; - __dbl_mp (res, &a, p); - __dbl_mp (0.5 * (res1 - res), &b, p); - __add (&a, &b, &c, p); - if (x > 2.4) - { - __sub (&pi, &c, &a, p); - __c32 (&a, &b, &c, p); - b.d[0] = -b.d[0]; - } - else if (x > 0.8) - { - __sub (&hp, &c, &a, p); - __c32 (&a, &c, &b, p); - } - else - __c32 (&c, &b, &a, p); /* b=cos(0.5*(res+res1)) */ - __dbl_mp (x, &c, p); /* c = x */ - __sub (&b, &c, &a, p); - /* if a > 0 return max (res, res1), otherwise return min (res, res1). */ - if ((a.d[0] > 0 && res <= res1) || (a.d[0] <= 0 && res >= res1)) - res = res1; - LIBC_PROBE (slowacos, 2, &res, &x); - return res; -} - /* Compute sin() of double-length number (X + DX) as Multi Precision number and return result as double. If REDUCE_RANGE is true, X is assumed to be the original input and DX is ignored. */ |