如何通过subprocess获取脚本计算数值的原精度结果？

问题分析 #

你用subprocess调用脚本获取变量值时，因为被调用脚本默认输出的是简化的数值字符串（比如Python里print(2/3)只会输出有限位数的十进制表示），转成float后无法还原原变量的精确浮点值，导致精度不一致。

解决方案核心思路 #

让被调用脚本输出能精确还原原数值的字符串表示，而非简化的打印结果，这样调用方解析后就能得到和原变量完全一致的数值。

情况1：被调用脚本是Python #

修改script.py，输出浮点数的完整精度表示：

# script.py
x = 2/3
print(repr(x))  # 输出能精确还原该浮点数的字符串

或者用固定精度格式化（双精度浮点数需要17位有效数字才能精确表示）：

# script.py
import sys
x = 2/3
print("{0:.{1}g}".format(x, sys.float_info.dig))

调用方解析代码：

import subprocess
import ast

cmd_line = ["python", "script.py"]
# 去除输出中的换行和空格
s = subprocess.check_output(cmd_line, text=True).strip()
# 用ast.literal_eval解析更安全，直接float(s)也可
x = ast.literal_eval(s)

情况2：被调用脚本是R（你的实际场景） #

在R脚本中输出17位有效数字的字符串（R的双精度浮点数同样需要17位来精确还原）：

# script.R
x <- 2/3
cat(sprintf("%.17g", x))

调用方解析代码和Python场景一致，获取字符串后转float即可得到和R中x完全相同的数值。

为什么原方法不行？ #

默认的print（Python）或print（R）会输出简化的十进制字符串，位数不足以精确还原双精度浮点数的二进制存储值。当你把这个简化字符串转成float时，得到的是一个近似值，和原变量存在微小差异。

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