diff options
author | Dan Albert <danalbert@google.com> | 2016-02-24 13:48:45 -0800 |
---|---|---|
committer | Dan Albert <danalbert@google.com> | 2016-02-24 13:51:18 -0800 |
commit | b9de1157289455b0ca26daff519d4a0ddcd1fa13 (patch) | |
tree | 4c56cc0a34b91f17033a40a455f26652304f7b8d /gcc-4.8.3/libquadmath/math/sinq.c | |
parent | 098157a754787181cfa10e71325832448ddcea98 (diff) | |
download | toolchain_gcc-b9de1157289455b0ca26daff519d4a0ddcd1fa13.tar.gz toolchain_gcc-b9de1157289455b0ca26daff519d4a0ddcd1fa13.tar.bz2 toolchain_gcc-b9de1157289455b0ca26daff519d4a0ddcd1fa13.zip |
Update 4.8.1 to 4.8.3.
My previous drop was the wrong version. The platform mingw is
currently using 4.8.3, not 4.8.1 (not sure how I got that wrong).
From ftp://ftp.gnu.org/gnu/gcc/gcc-4.8.3/gcc-4.8.3.tar.bz2.
Bug: http://b/26523949
Change-Id: Id85f1bdcbbaf78c7d0b5a69e74c798a08f341c35
Diffstat (limited to 'gcc-4.8.3/libquadmath/math/sinq.c')
-rw-r--r-- | gcc-4.8.3/libquadmath/math/sinq.c | 82 |
1 files changed, 82 insertions, 0 deletions
diff --git a/gcc-4.8.3/libquadmath/math/sinq.c b/gcc-4.8.3/libquadmath/math/sinq.c new file mode 100644 index 000000000..989f679d6 --- /dev/null +++ b/gcc-4.8.3/libquadmath/math/sinq.c @@ -0,0 +1,82 @@ +/* sinq.c -- __float128 version of s_sin.c. + * Conversion to long double by Jakub Jelinek, jj@ultra.linux.cz. + */ + +/* + * ==================================================== + * Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved. + * + * Developed at SunPro, a Sun Microsystems, Inc. business. + * Permission to use, copy, modify, and distribute this + * software is freely granted, provided that this notice + * is preserved. + * ==================================================== + */ + +/* sinq(x) + * Return sine function of x. + * + * kernel function: + * __quadmath_kernel_sinq ... sine function on [-pi/4,pi/4] + * __quadmath_kernel_cosq ... cose function on [-pi/4,pi/4] + * __quadmath_rem_pio2q ... argument reduction routine + * + * Method. + * Let S,C and T denote the sin, cos and tan respectively on + * [-PI/4, +PI/4]. Reduce the argument x to y1+y2 = x-k*pi/2 + * in [-pi/4 , +pi/4], and let n = k mod 4. + * We have + * + * n sin(x) cos(x) tan(x) + * ---------------------------------------------------------- + * 0 S C T + * 1 C -S -1/T + * 2 -S -C T + * 3 -C S -1/T + * ---------------------------------------------------------- + * + * Special cases: + * Let trig be any of sin, cos, or tan. + * trig(+-INF) is NaN, with signals; + * trig(NaN) is that NaN; + * + * Accuracy: + * TRIG(x) returns trig(x) nearly rounded + */ + +#include "quadmath-imp.h" + +__float128 +sinq (__float128 x) +{ + __float128 y[2],z=0.0Q; + int64_t n, ix; + + /* High word of x. */ + GET_FLT128_MSW64(ix,x); + + /* |x| ~< pi/4 */ + ix &= 0x7fffffffffffffffLL; + if(ix <= 0x3ffe921fb54442d1LL) + return __quadmath_kernel_sinq(x,z,0); + + /* sin(Inf or NaN) is NaN */ + else if (ix>=0x7fff000000000000LL) { + if (ix == 0x7fff000000000000LL) { + GET_FLT128_LSW64(n,x); + } + return x-x; + } + + /* argument reduction needed */ + else { + n = __quadmath_rem_pio2q(x,y); + switch(n&3) { + case 0: return __quadmath_kernel_sinq(y[0],y[1],1); + case 1: return __quadmath_kernel_cosq(y[0],y[1]); + case 2: return -__quadmath_kernel_sinq(y[0],y[1],1); + default: + return -__quadmath_kernel_cosq(y[0],y[1]); + } + } +} |