mpmath与Mathematica高精度计算结果存在差异，求排查原因

问题根源 #

你遇到的结果差异，核心原因是指数部分的表示精度问题：

• 在你的mpmath代码中，使用mp.mpf(str(.2))来表示0.2，但0.2是十进制分数1/5，无法精确转换为二进制浮点数（mpmath的mpf类型基于二进制浮点实现），因此这个值实际上是一个近似值，引入了微小误差。
• 这个误差经过指数函数放大后，最终导致计算结果从第20位小数开始出现偏差。

而Mathematica的计算中，SetPrecision[Exp[-64*.2], 300]会自动以高精度解析并处理-64*.2，将其识别为精确分数-64/5，从而得到更接近真实值的结果。

解决方案 #

在mpmath中，改用精确的分数形式表示指数，避免使用十进制小数的近似值：

from mpmath import mp
mp.dps = 300

# 用精确分数计算指数：-64/5 = -12.8
result = mp.exp(mp.mpf(-64) / mp.mpf(5))
print(result)

修改后的代码会生成与MathematicaN[Exp[-64/5], 300]完全一致的高精度结果，因为指数部分被精确表示，消除了近似误差。

补充说明 #

• 对于无法用有限二进制位表示的十进制小数（如0.2、0.1等），永远不要直接用字符串或浮点数初始化mpf，改用分数形式或精确的分子分母表示，比如mp.mpf(1)/5是精确的，而mp.mpf("0.2")是近似值。
• 处理复数高精度计算时，同样要确保复数的实部和虚部都用精确形式表示，避免引入不必要的近似误差。

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