Yes, aim was to be acurate down to 1 lsb while significantly faster. Feel free to drop terms from the polynomial if you can live with less accurate results!
The coefficients were generated by a package called Sollya, I've used it a few times to develop accurate chebyshev approximations for functions.
Please, Would you mind one of these days updating your blog post with the instructions you gave to sollya? I'm trying something stupid with log1p and can't get sollya to help, mostly because I'm not putting enough time to read all the docs...
Comments
Your result is significantly slower than the versions presented in the article, though yours has more terms and so may be more accurate.
Yes, aim was to be acurate down to 1 lsb while significantly faster. Feel free to drop terms from the polynomial if you can live with less accurate results!
The coefficients were generated by a package called Sollya, I've used it a few times to develop accurate chebyshev approximations for functions.
Abramowitz & Stegun is another good reference.
Please, Would you mind one of these days updating your blog post with the instructions you gave to sollya? I'm trying something stupid with log1p and can't get sollya to help, mostly because I'm not putting enough time to read all the docs...