[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
OT: negative zeros with ATLAS 3.8.4.x
From: |
Marco Caliari |
Subject: |
OT: negative zeros with ATLAS 3.8.4.x |
Date: |
Mon, 11 Apr 2011 17:07:38 +0200 (CEST) |
User-agent: |
Alpine 1.00 (DEB 882 2007-12-20) |
Dear maintainers,
the following code
clear all
m = 45;
y = linspace(0,1,m+1)';
y = y(1:m);
F = [0:m-1].'*2*pi*y.';
F(:,1)
produces some strange negative zeros (in random positions) with Octave
3.4.0 and ATLAS 3.8.4.x (x=1 or x=2, Core232SSE3), whereas everything is
fine with ATLAS 3.8.3. So, I suspect a bug in ATLAS (uninitialized
values?). Maybe someone is able to write a neat C code to submit to ATLAS
developers and/or to check that there is surely no fault in Octave.
Best regards,
Marco