Is there really no vectorized way to do this? We have a utility `sklearn.utils.extmath.logsumexp` that does a `logaddexp` over an entire array, and that's more stable than a series of `logaddexp`s.