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

Question

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

下面解释的那个，m = maxval:

Answer 1

您可以检查定义logsumexp here 的源代码。 （请注意，the doc page 上有指向源的链接）。

你会看到：

a_max = a.max(axis=0)
...
out = log(sum(exp(a - a_max), axis=0))

所以是的，scipy 的 logsumexp 正在从每个元素中减去最大值。