python - 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 是从每个元素中减去最大值。