蟒蛇麻木的机器

我试图理解什么是机器 ε。根据维基百科，它可以计算如下:

def machineEpsilon(func=float):
machine_epsilon = func(1)
while func(1)+func(machine_epsilon) != func(1):
machine_epsilon_last = machine_epsilon
machine_epsilon = func(machine_epsilon) / func(2)
return machine_epsilon_last

但是，它只适用于双精度数。我有兴趣修改它，以支持也单精度数字。我读到可以使用 numpy，特别是 numpy.float32类。有人能帮忙修改这个函数吗？

184465 次浏览

小开

就像大卫指出的那样，它已经起作用了！

>>> def machineEpsilon(func=float):
...     machine_epsilon = func(1)
...     while func(1)+func(machine_epsilon) != func(1):
...         machine_epsilon_last = machine_epsilon
...         machine_epsilon = func(machine_epsilon) / func(2)
...     return machine_epsilon_last
...
>>> machineEpsilon(float)
2.220446049250313e-16
>>> import numpy
>>> machineEpsilon(numpy.float64)
2.2204460492503131e-16
>>> machineEpsilon(numpy.float32)
1.1920929e-07

小开

最佳答案

为给定的 float 类型获得机器 ε 的一种更简单的方法是使用 np.finfo():

print(np.finfo(float).eps)
# 2.22044604925e-16


print(np.finfo(np.float32).eps)
# 1.19209e-07

小开

Another easy way to get epsilon is:

In [1]: 7./3 - 4./3 -1
Out[1]: 2.220446049250313e-16