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.

Best,

Filippo Vicentini