python : Geomspace and Logspace do not return true boundaries

Question

我正在尝试使用 python 3 和 numpy 创建一个对数刻度的数组。

我的数组边界有问题，返回的数组给出的下限低于我想要的限制。

例如:

In : a = np.array(379/19)
In : x = np.geomspace(a, 20)
In : x[0] >= a
Out: False

而使用 linspace 它可以工作:

In : y = np.linspace(a, 20)
In : y[0] >= a
Out: True

这是 logspace 和 geomspace 的已知问题吗？

感谢您的帮助!

Answer 1

你本质上是testing floats for equality ，这是一个坏主意。

a    is 19.94736842105263
x[0] is 19.947368421052627
y[0] is 19.94736842105263

所以它们非常接近，但并不完全相同。这就是 float 学通常的工作方式，由于舍入误差，您通常无法保证一个值是否大于或小于另一个值。

您到底想用这些边界做什么？