当使用两个1F1函数的和时,Tricomi汇合超几何函数的计算可能是病态的,因为它们的大小可能几乎相等,但符号相反。mpmath函数"hyperu“在内部使用任意精度,并在默认模式下生成35位有效数字的结果。这些数字中有多少是可靠的?它依赖于传递的参数吗?
import mpmath
x = mpmath.hyperu(a, b + 1, u)
发布于 2019-09-12 15:40:07
我刚刚收到mpmath的主要作者Fredrik Johansson发来的一封电子邮件,确认完整的35位数字是可用的。他写道:"hyperu在内部使用了自适应更高的精度,因此结果几乎总是准确到用户设置的全精度。“
https://stackoverflow.com/questions/57850844
复制相似问题