26    is scale*(1+TMP) without intermediate rounding.  The bit representation of
27    scale is in SBITS, however it has a computed exponent that may have
30    adjustment of scale, positive k here means the result may overflow and
35   double_t scale, y;  in specialcase()  local
39       /* k > 0, the exponent of scale might have overflowed by 1.  */  in specialcase()
41       scale = asdouble (sbits);  in specialcase()
42       y = 2 * (scale + scale * tmp);  in specialcase()
47   scale = asdouble (sbits);  in specialcase()
48   y = scale + scale * tmp;  in specialcase()
56       lo = scale - y + scale * tmp;  in specialcase()
83   double_t kd, r, r2, scale, tail, tmp;  in exp2()  local
114   /* 2^(k/N) ~= scale * (1 + tail).  */  in exp2()
118   /* This is only a valid scale when -1023*N < k < 1024*N.  */  in exp2()
120   /* exp2(x) = 2^(k/N) * 2^r ~= scale + scale * (tail + 2^r - 1).  */  in exp2()
134   scale = asdouble (sbits);  in exp2()
135   /* Note: tmp == 0 or |tmp| > 2^-65 and scale > 2^-928, so there  in exp2()
137   return eval_as_double (scale + scale * tmp);  in exp2()