On Wed, Mar 31, 2010 at 3:38 PM, David Warde-Farley <d...@cs.toronto.edu> wrote: > Unfortunately there's no good way of getting around order-of- > operations-related rounding error using the reduce() machinery, that I > know of. >
That seems reasonable, but receiving a nan, in this case, does not. Are my expectations are unreasonable? Would sorting the values before reducing be an ugly(-enough) solution? It seems to fix it in this particular case. Or is the method you posted in the link below a better solution? > I'd like to see a 'logsumexp' and 'logsumexp2' in NumPy instead, using > the generalized ufunc architecture, to do it over an arbitrary > dimension of an array, rather than needing to invoke 'reduce' on > logaddexp. I tried my hand at writing one but I couldn't manage to get > it working for some reason, and I haven't had time to revisit it: > http://mail.scipy.org/pipermail/numpy-discussion/2010-January/048067.html > I saw this original post and was hoping someone would post a response. Maybe now... _______________________________________________ NumPy-Discussion mailing list NumPy-Discussion@scipy.org http://mail.scipy.org/mailman/listinfo/numpy-discussion
