web-dev-qa-db-ja.com

python numpy machine epsilon

機械イプシロンとは何かを理解しようとしています。ウィキペディアによると、次のように計算できます。

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クラスを使用できることを読みました。誰でも関数の変更を手伝うことができますか？

pythonnumpyepsilon
84
Bob

特定の浮動小数点型のマシンイプシロンを取得する簡単な方法は、 np.finfo() を使用することです。

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

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

イプシロンを取得する別の簡単な方法は次のとおりです。

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

デビッドが指摘したように、すでに機能します！

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