scipy logsumexp()是否处理​​下溢挑战?

use*_*627 5 python arrays numpy scipy natural-logarithm

scipy的logsumexp() 实现是否包括通过从每个元素中减去数组中的最大找到值来防止下溢的hack?

下面解释的那个,其中m = maxval:

在此输入图像描述

unu*_*tbu 9

您可以检查logsumexp 此处定义的源代码.(请注意,doc页面上有一个指向源的链接).

你会看到的:

a_max = a.max(axis=0)
...
out = log(sum(exp(a - a_max), axis=0))
Run Code Online (Sandbox Code Playgroud)

所以是的,scipy's logsumexp减去每个元素的最大值.