14 #pragma warning ( disable : 4723 )
74 two54 = 1.80143985094819840000e+16,
75 Lg1 = 6.666666666666735130e-01,
76 Lg2 = 3.999999999940941908e-01,
77 Lg3 = 2.857142874366239149e-01,
78 Lg4 = 2.222219843214978396e-01,
79 Lg5 = 1.818357216161805012e-01,
80 Lg6 = 1.531383769920937332e-01,
81 Lg7 = 1.479819860511658591e-01;
83 static const double zero = 0.0;
87 double hfsq,
f,
s,
z,R,
w,
t1,t2,dk;
94 if (hx < 0x00100000) {
95 if (((hx&0x7fffffff)|lx)==0)
97 if (hx<0)
return (
x-
x)/
zero;
101 if (hx >= 0x7ff00000)
return x+
x;
104 i = (hx+0x95f64)&0x100000;
108 if((0x000fffff&(2+hx))<3) {
109 if(
f==
zero) {
if(
k==0)
return zero;
else {dk=(double)
k;
112 R =
f*
f*(0.5-0.33333333333333333*
f);
113 if(
k==0)
return f-R;
else {dk=(double)
k;
128 if(
k==0)
return f-(hfsq-
s*(hfsq+R));
else
131 if(
k==0)
return f-
s*(
f-R);
else
143 if (_LIB_VERSION == _IEEE_ || isnan(
x) ||
x > 0.0)
146 return __kernel_standard(
x,
x, 16);
147 return __kernel_standard(
x,
x, 17);