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正在减去每个元素的最大值。