[Numpy-discussion] Log Arrays
Thu May 8 12:57:19 CDT 2008
On Thu, May 8, 2008 at 12:53 PM, Charles R Harris
> On Thu, May 8, 2008 at 11:18 AM, David Cournapeau <email@example.com>
>> On Fri, May 9, 2008 at 2:06 AM, Nadav Horesh <firstname.lastname@example.org>
>> > Is the 80 bits float (float96 on IA32, float128 on AMD64) isn't enough?
>> > It has a 64 bits mantissa and can represent numbers up to nearly 1E(+-)5000.
>> It only make the problem happen later, I think. If you have a GMM with
>> million of samples of high dimension with many clusters, any "linear"
>> representation will fail I think. In a sense, the IEEE format is not
>> adequate for that kind of computation.
> David, what you are using is a log(log(x)) representation internally. IEEE
> is *not* linear, it is logarithmic.
*YES*. That is precisely the point. I want 53 bits devoted to the "x"
part of "exp(-x)". The straight IEEE representation is not logarithmic
"I have come to believe that the whole world is an enigma, a harmless
enigma that is made terrible by our own mad attempt to interpret it as
though it had an underlying truth."
-- Umberto Eco
More information about the Numpy-discussion