SNG_EXPBITS
exp = (i >> (32 - 1 - SNG_EXPBITS)) & mask(SNG_EXPBITS);
s = e >> SNG_EXPBITS;
result = hs << (DBL_EXPBITS - SNG_EXPBITS);