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);
}
