[Numpy-discussion] Log Arrays

[Charles R Harris]

On Thu, May 8, 2008 at 11:18 AM, David Cournapeau <cournape at gmail.com>
wrote:

> On Fri, May 9, 2008 at 2:06 AM, Nadav Horesh <nadavh at visionsense.com>
> wrote:
> > 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.

Chuck
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mail.python.org/pipermail/numpy-discussion/attachments/20080508/f69f4497/attachment.html>