python mesin numpy epsilon

Saya mencoba memahami apa itu mesin epsilon. Menurut Wikipedia, itu dapat dihitung sebagai berikut:

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

Namun, ini hanya cocok untuk angka presisi ganda. Saya tertarik untuk memodifikasinya untuk mendukung juga angka presisi tunggal. Saya membaca bahwa numpy dapat digunakan, terutama di numpy.float32kelas. Adakah yang bisa membantu dengan memodifikasi fungsi?

Cara yang lebih mudah untuk mendapatkan mesin epsilon untuk tipe float tertentu adalah dengan menggunakan np.finfo():

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

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

Cara lain yang mudah untuk mendapatkan epsilon adalah:

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

Itu sudah akan berhasil, seperti yang ditunjukkan David!

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