/* Exponential base 2 function.
- Copyright (C) 2011-2012 Free Software Foundation, Inc.
+ Copyright (C) 2011-2013 Free Software Foundation, Inc.
This program is free software: you can redistribute it and/or modify
it under the terms of the GNU General Public License as published by
+ 21844/6081075 * z^13
- 929569/638512875 * z^15
+ ...
- Since |z| <= log(2)/1024 < 0.0007, the relative error of the z^13 term
- is < 0.0007^12 < 2^-120 <= 2^-LDBL_MANT_DIG, therefore we can truncate
- the series after the z^11 term. */
+ Since |z| <= log(2)/1024 < 0.0007, the relative contribution of the
+ z^13 term is < 0.0007^12 < 2^-120 <= 2^-LDBL_MANT_DIG, therefore we
+ can truncate the series after the z^11 term. */
{
long double nm = roundl (x * 256.0L); /* = 256 * n + m */