On Thu, Nov 6, 2008 at 1:23 PM, T J
On Wed, Nov 5, 2008 at 2:09 PM, Charles R Harris
wrote: I'm inclined to go with logaddexp and add logsumexp as an alias for logaddexp.reduce. But I'll wait until tomorrow to see if there are more comments.
When working in other bases, it seems like it would be good to avoid having to convert to base e and then back to base 2 with each function call (for example). Is there any desire add similar functions for other standard bases?
I suppose that depends on who you ask ;) What is your particular interest in these other bases and why would they be better than working in base e and converting at the end? The only one I could see really having a fast implementation is log2. In fact, I think the standard log starts in log2 by pulling in the floating point exponent and then using some sort of rational approximation of log2 over the range [1,2) on the mantissa. Chuck