45{
46
47
48
49 static double piby2 = 1.5707963267948966e+00;
50 double chi, clo,
v,
s,
q,
z;
51
52
53
54 unsigned long long ux,
aux, xneg;
56 aux = ux & ~SIGNBIT_DP64;
58
61
62
63
64 if (
aux > 0x4003800000000000)
65 {
66
68 {
69
72 }
73 else if (
v > 0x4370000000000000)
74 {
75
76 if (xneg)
78 else
80 }
81
83
84 chi = 1.57079632679489655800e+00;
85 clo = 6.12323399573676480327e-17;
86 }
87 else if (
aux > 0x3ff3000000000000)
88 {
89 x = (
v-1.5)/(1.0+1.5*
v);
90
91 chi = 9.82793723247329054082e-01;
92 clo = 1.39033110312309953701e-17;
93 }
94 else if (
aux > 0x3fe6000000000000)
95 {
97
98 chi = 7.85398163397448278999e-01;
99 clo = 3.06161699786838240164e-17;
100 }
101 else if (
aux > 0x3fdc000000000000)
102 {
103 x = (2.0*
v-1.0)/(2.0+
v);
104
105 chi = 4.63647609000806093515e-01;
106 clo = 2.26987774529616809294e-17;
107 }
108 else
109 {
111 chi = 0.0;
112 clo = 0.0;
113 }
114
115
116
119 (0.268297920532545909e0 +
120 (0.447677206805497472e0 +
121 (0.220638780716667420e0 +
122 (0.304455919504853031e-1 +
123 0.142316903342317766e-3*
s)*
s)*
s)*
s)/
124 (0.804893761597637733e0 +
125 (0.182596787737507063e1 +
126 (0.141254259931958921e1 +
127 (0.424602594203847109e0 +
128 0.389525873944742195e-1*
s)*
s)*
s)*
s);
129
130 z = chi - ((
q - clo) -
x);
131
134}
double __cdecl _handle_error(char *fname, int opcode, unsigned long long value, int type, int flags, int error, double arg1, double arg2, int nargs)
Handles an error condition.
GLint GLint GLint GLint GLint x
GLdouble GLdouble GLdouble GLdouble q
#define GET_BITS_DP64(x, ux)