36 4.63647609000806093515e-01,
37 7.85398163397448278999e-01,
38 9.82793723247329054082e-01,
39 1.57079632679489655800e+00,
43 2.26987774529616870924e-17,
44 3.06161699786838301793e-17,
45 1.39033110312309984516e-17,
46 6.12323399573676603587e-17,
49 static const double aT[] = {
50 3.33333333333329318027e-01,
51 -1.99999999998764832476e-01,
52 1.42857142725034663711e-01,
53 -1.11111104054623557880e-01,
54 9.09088713343650656196e-02,
55 -7.69187620504482999495e-02,
56 6.66107313738753120669e-02,
57 -5.83357013379057348645e-02,
58 4.97687799461593236017e-02,
59 -3.65315727442169155270e-02,
60 1.62858201153657823623e-02,
78 (ix==0x7ff00000&&(low!=0)))
81 else return -
atanhi[3]-atanlo[3];
82 }
if (ix < 0x3fdc0000) {
83 if (ix < 0x3e200000) {
89 if (ix < 0x3ff30000) {
90 if (ix < 0x3fe60000) {
91 id = 0; x = (2.0*x-
one)/(2.0+x);
93 id = 1; x = (x-
one)/(x+
one);
96 if (ix < 0x40038000) {
97 id = 2; x = (x-1.5)/(
one+1.5*
x);
106 s1 = z*(
aT[0]+w*(
aT[2]+w*(
aT[4]+w*(
aT[6]+w*(
aT[8]+w*
aT[10])))));
107 s2 = w*(
aT[1]+w*(
aT[3]+w*(
aT[5]+w*(
aT[7]+w*
aT[9]))));
108 if (
id<0)
return x - x*(s1+s2);
#define GET_HIGH_WORD(i, d)
static const double atanhi[]
#define libm_hidden_def(x)
GLuint GLfloat GLfloat GLfloat GLfloat GLfloat GLfloat GLfloat GLfloat s1
GLint GLint GLint GLint z
GLint GLint GLint GLint GLint w
static const double atanlo[]