38 static const double one = 1.0, half=0.5, huge = 1.0e300;
40 static double one = 1.0, half=0.5, huge = 1.0e300;
59 if(ix>=0x7ff00000)
return x*x;
65 if (ix<0x3c800000)
return w;
66 return one+(t*t)/(w+w);
70 if (ix < 0x40360000) {
79 lx = *( (((*(
unsigned*)&
one)>>29)) + (
unsigned*)&x);
80 if ((ix<0x408633CE) ||
81 ((ix==0x408633ce)&&(lx<=(
unsigned)0x8fb9f87d))) {