long double __tanl(long double x) { long double y[2],z=0.0L; int64_t n, ix; /* High word of x. */ GET_LDOUBLE_MSW64(ix,x); /* |x| ~< pi/4 */ ix &= 0x7fffffffffffffffLL; if(ix <= 0x3ffe921fb54442d1LL) return __kernel_tanl(x,z,1); /* tanl(Inf or NaN) is NaN */ else if (ix>=0x7fff000000000000LL) { if (ix == 0x7fff000000000000LL) { GET_LDOUBLE_LSW64(n,x); if (n == 0) __set_errno (EDOM); } return x-x; /* NaN */ } /* argument reduction needed */ else { n = __ieee754_rem_pio2l(x,y); return __kernel_tanl(y[0],y[1],1-((n&1)<<1)); /* 1 -- n even -1 -- n odd */ } }
long double tanl(long double x) { union IEEEl2bits z; int i, e0, s; double xd[NX], yd[PREC]; long double hi, lo; z.e = x; s = z.bits.sign; z.bits.sign = 0; /* If x = +-0 or x is subnormal, then tan(x) = x. */ if (z.bits.exp == 0) return (x); /* If x = NaN or Inf, then tan(x) = NaN. */ if (z.bits.exp == 32767) return ((x - x) / (x - x)); /* Optimize the case where x is already within range. */ if (z.e < M_PI_4) { hi = __kernel_tanl(z.e, 0, 0); return (s ? -hi : hi); } /* Split z.e into a 24-bit representation. */ e0 = ilogbl(z.e) - 23; z.e = scalbnl(z.e, -e0); for (i = 0; i < NX; i++) { xd[i] = (double)((int32_t)z.e); z.e = (z.e - xd[i]) * two24; } /* yd contains the pieces of xd rem pi/2 such that |yd| < pi/4. */ e0 = __kernel_rem_pio2(xd, yd, e0, NX, PREC); #if PREC == 2 hi = (long double)yd[0] + yd[1]; lo = yd[1] - (hi - yd[0]); #else /* PREC == 3 */ long double t; t = (long double)yd[2] + yd[1]; hi = t + yd[0]; lo = yd[0] - (hi - t); #endif switch (e0 & 3) { case 0: case 2: hi = __kernel_tanl(hi, lo, 0); break; case 1: case 3: hi = __kernel_tanl(hi, lo, 1); break; } return (s ? -hi : hi); }
long double __tanl(long double x) { long double y[2],z=0.0; int32_t n, se, i0, i1; /* High word of x. */ GET_LDOUBLE_WORDS(se,i0,i1,x); /* |x| ~< pi/4 */ se &= 0x7fff; if(se <= 0x3ffe) return __kernel_tanl(x,z,1); /* tan(Inf or NaN) is NaN */ else if (se==0x7fff) { if (i1 == 0 && i0 == 0x80000000) __set_errno (EDOM); return x-x; } /* argument reduction needed */ else { n = __ieee754_rem_pio2l(x,y); return __kernel_tanl(y[0],y[1],1-((n&1)<<1)); /* 1 -- n even -1 -- n odd */ } }
long double tanl(long double x) { union IEEEl2bits z; int e0, s; long double y[2]; long double hi, lo; z.e = x; s = z.bits.sign; z.bits.sign = 0; /* If x = +-0 or x is subnormal, then tan(x) = x. */ if (z.bits.exp == 0) return (x); /* If x = NaN or Inf, then tan(x) = NaN. */ if (z.bits.exp == 32767) return ((x - x) / (x - x)); ENTERI(); /* Optimize the case where x is already within range. */ if (z.e < M_PI_4) { hi = __kernel_tanl(z.e, 0, 0); RETURNI(s ? -hi : hi); } e0 = __ieee754_rem_pio2l(x, y); hi = y[0]; lo = y[1]; switch (e0 & 3) { case 0: case 2: hi = __kernel_tanl(hi, lo, 0); break; case 1: case 3: hi = __kernel_tanl(hi, lo, 1); break; } RETURNI(hi); }