HFRAC_BITS
while(hx<(1LL<<HFRAC_BITS)) { /* normalize x */
while(hx<(int64_t)(1ULL<<HFRAC_BITS)) { /* normalize x */