我正在使用python 3和numpy,我正在尝试使用np.geomspace计算数组。 我的问题是数组不尊重我想要的边界。
例如我有:
Polygon
差异非常小,但是我需要这个数组来显示插值函数,并且它会引发错误:“ValueError:x_new中的值低于插值范围。”,因为geomspace数组在之前启动我问的限制。
当我使用logspace时,这个问题是一样的:
In : np.geomspace(2.1951999999999994e-08, 1)[0]
Out : 2.1951999999999987e-08
我不知道python浮点数是否存在精度问题,我试图寻找解决方案但没有成功。
感谢您的帮助。
感谢您告诉我这是一个用计算机定义浮动的问题。
但这并没有真正帮助我,因为事实上我作为一个例子展示的数字来自另一个带有python函数的计算,即:
In [243]: np.logspace(np.log10(2.1951999999999994e-08), 1)[0]
Out[243]: 2.1951999999999987e-08
其中a和b是两个float64 numpy数组。
所以当我这样做时:
max(a[0], b[0])
我的答案小于最大值(a [0],b [0])。 这是绘制插值的一个大问题。