3 /* k minimizes |exp(k ln2) - 2**k|, see cmath/__cexpf.c */
4 static const uint32_t k = 235;
5 static const float kln2 = 162.88958740f;
7 /* expf(x)/2 for large x */
8 float __expo2f(float x)
12 /* note that k is odd and scale*scale overflows */
13 SET_FLOAT_WORD(scale, (0x7f + k/2) << 23);
14 /* exp(x - k ln2) * 2**(k-1) */
15 return expf(x - kln2) * scale * scale;