如何在Python中获取Fortran级别的浮点数精度？

嘿，我明白你想要在Python里复刻Fortran那种更高精度的π值的需求——毕竟numpy默认的双精度π（numpy.pi，值为3.141592653589793）确实比Fortran里高精度版本（3.141592653589793238462643...）短几位。咱们来聊聊可行的解决方案：

先搞清楚精度差异的根源 #

Fortran里你用到的高精度π，对应的是128位扩展精度浮点数（通常用real(kind=16)声明），而numpy.pi是标准的64位双精度浮点数。直接用Python内置的float或者普通的numpy类型，都达不到那个精度级别。

解决方案1：用mpmath实现任意精度（最通用靠谱） #

mpmath是Python专门用于任意精度数学计算的库，完全可以生成和Fortran高精度π一致的值，甚至能设置更高的精度。

步骤如下：

1. 先安装库：

   pip install mpmath
   

2. 设置精度并获取高精度π：

   import mpmath as mp
   # 设置小数点后的精度位数，这里设为50位，足够覆盖你提到的Fortran版本的π
   mp.mp.dps = 50
   high_precision_pi = mp.pi
   print(high_precision_pi)
   # 输出会是：3.14159265358979323846264338327950288419716939937510
   

   你可以根据需求调整mp.mp.dps的值，设置得越大，精度越高。

解决方案2：正确使用numpy的float128（注意平台限制） #

如果你想用numpy的float128，要注意一个坑：numpy的float128在不同平台上的实现不一样。比如在x86架构的Linux系统上，它对应的是硬件支持的80位扩展精度；而在一些其他平台（比如macOS）可能没有这个类型，或者是模拟的128位。

关键是：不能直接把numpy.pi转成float128（那样只是把双精度值扩展，不会提升精度），而是要从高精度的字符串转换：

import numpy as np
# 用你提到的Fortran级别的π字符串
pi_high_precision_str = "3.14159265358979323846264338327950288419716939937510"
pi_float128 = np.float128(pi_high_precision_str)
print(pi_float128)

这样就能得到和Fortran一致精度的π值了，但要注意平台兼容性问题。

为什么你之前试float128没成功？ #

大概率是你直接把numpy.pi转成了float128，但numpy.pi本身只有64位精度，转成128位也只是补零，不会增加有效数字。必须从更高精度的数据源（比如字符串或者专门的高精度库生成）转换才行。

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