python - 如何处理python中64float精度的指数溢出？

Question

我是 python 的新手，很抱歉这个简单的问题。

在下面的代码中，我想计算指数然后取对数。

Y=numpy.log(1+numpy.exp(1000))

问题是，当我取 710 或更大数字的指数时，numpy.exp() 函数会返回“inf”，即使我使用 64float 打印它也会打印“inf”。

任何有关此问题的帮助将不胜感激。

Answer 1

您可以使用函数 np.logaddexp() 来执行此类操作。它计算 logaddexp(x1, x2) == log(exp(x1) + exp(x2)) 而不显式计算中间 exp() 值。这样就避免了溢出。由于 exp(0.0) == 1，您将按预期计算 np.logaddexp(0.0, 1000.0) 并获得 1000.0 的结果.