logaritmo - log array python

Does scipy logsumexp() lidiar con el desafío de subdesbordamiento?

Question

¿La implementación de logsumexp() de logsumexp() incluye el hack que previene el underflow al restar el máximo valor encontrado en la matriz de cada elemento?

El que se explica a continuación, donde m = maxval :

Answer 1

Puede inspeccionar el código fuente que define logsumexp aquí . (Tenga en cuenta que hay un enlace a la fuente en la página del documento ).

Verás:

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

Así que sí, logsumexp de logsumexp está restando el máximo de cada elemento.