3 /* k minimizes |exp(k ln2) - 2**k|, see cmath/__cexp.c */
4 static const uint32_t k = 1799;
5 static const double kln2 = 1246.97177782734161156;
7 /* exp(x)/2 for large x */
8 double __expo2(double x)
12 /* note that k is odd and scale*scale overflows */
13 INSERT_WORDS(scale, (0x3ff + k/2) << 20, 0);
14 /* exp(x - k ln2) * 2**(k-1) */
15 return exp(x - kln2) * scale * scale;