trăn numpy máy epsilon

Tôi đang cố gắng hiểu epsilon máy là gì. Theo Wikipedia, nó có thể được tính như sau:

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

Tuy nhiên, nó chỉ phù hợp với các số chính xác gấp đôi. Tôi quan tâm đến việc sửa đổi nó để hỗ trợ các số chính xác đơn lẻ. Tôi đọc rằng numpy có thể được sử dụng, đặc biệt là numpy.float32lớp học. Ai có thể giúp sửa đổi chức năng?

Một cách dễ dàng hơn để có được epsilon của máy cho một loại phao nhất định là sử dụng np.finfo():

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

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

Một cách dễ dàng khác để có được epsilon là:

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

Nó sẽ hoạt động, như David đã chỉ ra!

>>> 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