Имеет ли scipy logsumexp() дело с проблемой underflow?

Question

Реализована ли реализация scipy logsumexp() взлома, которая предотвращает переполнение, вычитая максимальное найденное значение в массиве из каждого элемента?

Одно объяснено здесь ниже, где m = maxval:

Answer 1

Вы можете проверить исходный код, определяющий logsumexphere. (Обратите внимание, что есть ссылка на источник по адресу the doc page).

Вы увидите:

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

Так что да, SciPy-х logsumexp вычитает максимум из каждого элемента.