Lines Matching +full:39 +full:x
12 /* atan(x)
14 * 1. Reduce x to positive by atan(x) = -atan(-x).
15 * 2. According to the integer k=4t+0.25 chopped, t=x, the argument
19 * [0,7/16] atan(x) = t-t^3*(a1+t^2*(a2+...(a10+t^2*a11)...)
20 * [7/16,11/16] atan(x) = atan(1/2) + atan( (t-0.5)/(1+t/2) )
21 * [11/16.19/16] atan(x) = atan( 1 ) + atan( (t-1)/(1+t) )
22 * [19/16,39/16] atan(x) = atan(3/2) + atan( (t-1.5)/(1+1.5t) )
23 * [39/16,INF] atan(x) = atan(INF) + atan( -1/t )
70 atan(double x) in atan() argument
75 GET_HIGH_WORD(hx,x); in atan()
77 if(ix>=0x44100000) { /* if |x| >= 2^66 */ in atan()
79 GET_LOW_WORD(low,x); in atan()
82 return x+x; /* NaN */ in atan()
85 } if (ix < 0x3fdc0000) { /* |x| < 0.4375 */ in atan()
86 if (ix < 0x3e400000) { /* |x| < 2^-27 */ in atan()
87 if(huge+x>one) return x; /* raise inexact */ in atan()
91 x = fabs(x); in atan()
92 if (ix < 0x3ff30000) { /* |x| < 1.1875 */ in atan()
93 if (ix < 0x3fe60000) { /* 7/16 <=|x|<11/16 */ in atan()
94 id = 0; x = (2.0*x-one)/(2.0+x); in atan()
95 } else { /* 11/16<=|x|< 19/16 */ in atan()
96 id = 1; x = (x-one)/(x+one); in atan()
99 if (ix < 0x40038000) { /* |x| < 2.4375 */ in atan()
100 id = 2; x = (x-1.5)/(one+1.5*x); in atan()
101 } else { /* 2.4375 <= |x| < 2^66 */ in atan()
102 id = 3; x = -1.0/x; in atan()
106 z = x*x; in atan()
111 if (id<0) return x - x*(s1+s2); in atan()
113 z = atanhi[id] - ((x*(s1+s2) - atanlo[id]) - x); in atan()