gpt4 book ai didi

python - 如何设置浮点精度以避免 Python 中的浮点运算错误?

转载 作者:行者123 更新时间:2023-12-04 03:57:14 24 4
gpt4 key购买 nike

我一直在 Python 中遇到一个我似乎无法弄清楚的浮点运算错误。

问题:我需要创建一个权重,使所有权重总和为 1,而不是,例如:0.99999999999999。

例如,如下代码:

values = numpy.array([9626.40000000034,      0. ,      0. ,      0. ,      0. ,      0. ,
0. , 0. , 36907.300000000000054])
weights = values/values.sum()
weights.sum()

产量:

0.99999999999999989

而不是 1。我试过乘以 1000,转换为字符串(以切断精度),然后转换回 float 并除以 1000。它不起作用。我也尝试过使用 Decimal。

from decimal import *
string_weight = []
float_weight = []
getcontext().prec = 3
for number in weights:
string_weight.append(Decimal(str(number)))
for string in string_weight:
float_weight.append(float(string))
fuel_weights = numpy.array(fuel_weights_float)
fuel_weights.sum()

答案是:

1.0009999999999999

这不是我想要的。我只想要一个简单的“1.0”。

sys.version 报告给出:

3.6.8 |Anaconda, Inc.| (default, Dec 29 2018, 19:04:46) 
[GCC 4.2.1 Compatible Clang 4.0.1 (tags/RELEASE_401/final)]

我正在使用 Mac OS X Catalina。

最佳答案

这个二进制算术问题的解决方案是使用 Decimal,这里是如何正确使用它。

首先,让我分享一个更清晰的问题示例。

import numpy
from decimal import *

# creating dummy values and weights
values = 1e-10 * numpy.ones(5)
weights = values/values.sum()
weights.sum()

产量:

0.9999999999999999

当我应用 Decimal 来解决这个错误时,我遇到了一个对我来说不直观的细微差别:通过 Decimal 传递除法运算的结果不会对二进制错误做任何事情,即:

getcontext().prec = 5
Decimal(values[0])

产量:

Decimal('1.0000000000000000364321973154977415791655470655996396089904010295867919921875E-10')

为了让 Decimal 纠正二进制算术错误,Decimal 必须包含在除法运算中,如下所示。

getcontext().prec = 5
Decimal(1)/Decimal(7)

产生与 Python docs website 中所述相同的答案:

Decimal('0.14286')

在我的例子中,正确的应用程序看起来像这样。

weights_list = []
values = 1e-10 * numpy.ones(5)
sum_values = values.sum()
for value in values:
getcontext().prec = 5
weight = Decimal(value)/Decimal(sum_values)
weights_list.append(weight)
weights = numpy.array(weights_list)
weights.sum()

结果是正确的数学答案,而不是二进制算术答案。

Decimal('1.0000')

可以使用 numpy.float() 将其转换为数字。

关于python - 如何设置浮点精度以避免 Python 中的浮点运算错误?,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/63654072/

24 4 0
Copyright 2021 - 2024 cfsdn All Rights Reserved 蜀ICP备2022000587号
广告合作:1813099741@qq.com 6ren.com