关于mpmath库hyperu函数默认精度及可靠位数的技术咨询

针对你的两个问题，我来逐一拆解说明：

1. hyperu函数的默认精度 #

mpmath库的全局默认精度是35位有效数字，hyperu函数作为库内的任意精度函数，会直接遵循这个默认设置。如果你需要调整精度，可以通过修改mpmath.mp.dps的值来实现，比如mpmath.mp.dps = 50就会把全局精度设为50位有效数字，后续的hyperu计算都会使用这个精度。

2. 可靠有效数字与参数的关系 #

首先要明确：直接用两个₁F₁函数求和计算Tricomi合流超几何函数时的病态问题，根源是两个函数值大小接近但符号相反，导致严重的数值抵消，丢失有效数字。但mpmath的hyperu函数内部并没有采用这种直接求和的方式，它使用了更稳定的算法（比如针对不同参数范围选择渐近展开、级数展开或者其他避免抵消的计算路径），从根本上规避了这个病态问题。

回到你的问题：默认35位有效数字的结果中，可靠位数确实依赖于传入的参数：

• 如果参数处于极端场景（比如某些会触发数值不稳定的a、b、u组合），可靠位数会有所减少，但mpmath的任意精度机制会自动预留足够的精度冗余，尽可能保证结果的可靠性；
• 在绝大多数常规参数下，35位有效数字里的绝大多数都是可靠的，足以满足大部分科学计算的需求。

代码示例 #

import mpmath

# 定义参数并计算hyperu
a = 1.5
b = 2.0
u = 3.0
result = mpmath.hyperu(a, b + 1, u)
print(f"计算结果：{result}")

内容的提问来源于stack exchange，提问作者tcp