+ 21844/6081075 * z^13
- 929569/638512875 * z^15
+ ...
- Since |z| <= log(2)/1024 < 0.0007, the relative error of the z^7 term
- is < 0.0007^6 < 2^-60 <= 2^-DBL_MANT_DIG, therefore we can truncate
- the series after the z^5 term. */
+ Since |z| <= log(2)/1024 < 0.0007, the relative contribution of the
+ z^7 term is < 0.0007^6 < 2^-60 <= 2^-DBL_MANT_DIG, therefore we can
+ truncate the series after the z^5 term. */
{
double nm = round (x * 256.0); /* = 256 * n + m */