I would like to ask if the maintainers would be in favour of adding
support to np.logaddexp in order to support complex numbers.
scipy.special.logsumexp already does, and it would be relatively
straightforward to support them in numpy too.
Complex numbers are very common when working in quantum physics, where
we have complex-valued probability amplitudes.