28    is scale*(1+TMP) without intermediate rounding.  The bit representation of
29    scale is in SBITS, however it has a computed exponent that may have
32    adjustment of scale, positive k here means the result may overflow and
37   double_t scale, y;  in specialcase()  local
41       /* k > 0, the exponent of scale might have overflowed by <= 460.  */  in specialcase()
43       scale = asdouble (sbits);  in specialcase()
44       y = 0x1p1009 * (scale + scale * tmp);  in specialcase()
49   scale = asdouble (sbits);  in specialcase()
50   y = scale + scale * tmp;  in specialcase()
58       lo = scale - y + scale * tmp;  in specialcase()
87   double_t kd, z, r, r2, scale, tail, tmp;  in exp_inline()  local
132   /* 2^(k/N) ~= scale * (1 + tail).  */  in exp_inline()
136   /* This is only a valid scale when -1023*N < k < 1024*N.  */  in exp_inline()
138   /* exp(x) = 2^(k/N) * exp(r) ~= scale + scale * (tail + exp(r) - 1).  */  in exp_inline()
152   scale = asdouble (sbits);  in exp_inline()
153   /* Note: tmp == 0 or |tmp| > 2^-200 and scale > 2^-739, so there  in exp_inline()
155   return eval_as_double (scale + scale * tmp);  in exp_inline()