我要求一个简短的解释,指出所涉及的各种加速方法。这仅仅是为了好奇。
例如mpmath website告诉指数函数公式用于复平面中的操作,但是现在我想将它保留在普通实数上使用的操作中。
可选地,有一个链接到源代码部分的链接,其中实际计算以mpmath执行,我没有运气搜索可下载的源。
答案 0 :(得分:4)
最多1000个十进制数字的精度基本函数通常由Taylor series计算,收敛性很好。如果你明确知道你需要什么样的精度,Chebyshev polynomials以最少的努力获得最好的质量。
为了获得非凡的精确度,您使用的是Arithmetic-geometric mean,它以二次方式收敛到解决方案,并允许您快速计算ln(x)。逆函数exp(x)通过Newton-Raphson迭代计算。从那以后,有一些算法能够计算arctan,tan,sin,cos ......
有关算法的闲暇,请参阅Jean-Michel Muller的Elementary Functions: Algorithms and Implementation。