29    is scale*(1+TMP) without intermediate rounding.  The bit representation of
30    scale is in SBITS, however it has a computed exponent that may have
33    adjustment of scale, positive k here means the result may overflow and
38   double_t scale, y;  in exp_inline_special_case()  local
42       /* k > 0, the exponent of scale might have overflowed by <= 460.  */  in exp_inline_special_case()
44       scale = asdouble (sbits);  in exp_inline_special_case()
45       y = 0x1p1009 * (scale + scale * tmp);  in exp_inline_special_case()
50   scale = asdouble (sbits);  in exp_inline_special_case()
51   y = scale + scale * tmp;  in exp_inline_special_case()
59       lo = scale - y + scale * tmp;  in exp_inline_special_case()
88   double_t kd, z, r, r2, scale, tail, tmp;  in exp_inline()  local
133   /* 2^(k/N) ~= scale * (1 + tail).  */  in exp_inline()
137   /* This is only a valid scale when -1023*N < k < 1024*N.  */  in exp_inline()
139   /* exp(x) = 2^(k/N) * exp(r) ~= scale + scale * (tail + exp(r) - 1).  */  in exp_inline()
153   scale = asdouble (sbits);  in exp_inline()
154   /* Note: tmp == 0 or |tmp| > 2^-200 and scale > 2^-739, so there  in exp_inline()
156   return eval_as_double (scale + scale * tmp);  in exp_inline()