Python:Geomspace(类似np.logspace)不尊重边界

时间:2017-05-23 13:55:13

标签: python arrays python-3.x numpy

我正在使用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])。 这是绘制插值的一个大问题。

0 个答案:

没有答案