python numpy machine epsilon

Question

Próbuję zrozumieć, czym jest epsilon maszyny. Według Wikipedii można ją obliczyć w następujący sposób:

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 

Jednak nadaje się tylko do liczb podwójnej precyzji. Interesuje mnie modyfikowanie go tak, aby obsługiwał również liczby pojedynczej precyzji. Czytałem, że można używać numpy, szczególnie klasy numpy.float32. Czy ktokolwiek może pomóc w modyfikacji funkcji?

Answer 1

Łatwiejszy sposób, aby uzyskać epsilon maszynowy dla danego typu float jest użycie np.finfo():

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

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

Answer 2

Już będzie działać, jak zauważył Dawid!

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

Answer 3

Kolejny łatwy sposób uzyskać Epsilon to:

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