33 if((uint64_t)v+v > 0xFFEULL<<52)
35 return ldexp(((v&((1LL<<52)-1)) + (1LL<<52)) * (v>>63|1), (v>>52&0x7FF)-1075);
39 if((uint32_t)v+v > 0xFF000000U)
41 return ldexp(((v&0x7FFFFF) + (1<<23)) * (v>>31|1), (v>>23&0xFF)-150);
48 for (i = 0; i < 8; i++)
64 else if(d-d)
return 0x7FF0000000000000LL + ((int64_t)(d<0)<<63) + (d!=d);
66 return (int64_t)(d<0)<<63 | (e+1022LL)<<52 | (int64_t)((fabs(d)-0.5)*(1LL<<53));
72 else if(d-d)
return 0x7F800000 + ((d<0)<<31) + (d!=d);
74 return (d<0)<<31 | (e+126)<<23 | (int64_t)((fabs(d)-0.5)*(1<<24));
79 int e, i;
double f; uint64_t m;
81 f = fabs(frexp(d, &e));
82 if (f >= 0.5 && f < 1) {
86 m = (uint64_t)ldexp(f, 64);
89 }
else if (f != 0.0) {