Re: [Numpy-discussion] float128 / longdouble on PPC - is it broken?

2011-10-26 Thread David Cournapeau
On Wed, Oct 26, 2011 at 12:49 AM, Matthew Brett matthew.br...@gmail.com wrote: That reminds me of the old joke about the Irishman giving directions - If I were you, I wouldn't start from here. Sounds about accurate `1 I would rather spend some time on implementing/integrating portable quad

Re: [Numpy-discussion] float128 / longdouble on PPC - is it broken?

2011-10-26 Thread Nathaniel Smith
On Tue, Oct 25, 2011 at 4:49 PM, Matthew Brett matthew.br...@gmail.com wrote: I guess from your answer that such a warning would be complicated to implement, and if that's the case, I can imagine it would be low priority. I assume the problem is more that it would be a weirdo check that

Re: [Numpy-discussion] float128 / longdouble on PPC - is it broken?

2011-10-26 Thread Matthew Brett
Hi, On Wed, Oct 26, 2011 at 1:07 AM, Nathaniel Smith n...@pobox.com wrote: On Tue, Oct 25, 2011 at 4:49 PM, Matthew Brett matthew.br...@gmail.com wrote: I guess from your answer that such a warning would be complicated to implement, and if that's the case, I can imagine it would be low

Re: [Numpy-discussion] float128 / longdouble on PPC - is it broken?

2011-10-26 Thread Pauli Virtanen
26.10.2011 10:07, Nathaniel Smith kirjoitti: On Tue, Oct 25, 2011 at 4:49 PM, Matthew Brett matthew.br...@gmail.com wrote: I guess from your answer that such a warning would be complicated to implement, and if that's the case, I can imagine it would be low priority. I assume the problem

Re: [Numpy-discussion] float128 / longdouble on PPC - is it broken?

2011-10-25 Thread Pauli Virtanen
25.10.2011 06:59, Matthew Brett kirjoitti: res = np.longdouble(2)**64 res-1 36893488147419103231.0 Can you check if long double works properly (not a given) in C on that platform: long double x; x = powl(2, 64); x -= 1; printf(%g %Lg\n, (double)x, x); or, in

Re: [Numpy-discussion] float128 / longdouble on PPC - is it broken?

2011-10-25 Thread Charles R Harris
On Mon, Oct 24, 2011 at 10:59 PM, Matthew Brett matthew.br...@gmail.comwrote: Hi, I just ran into this on a PPC machine: In [1]: import numpy as np In [2]: np.__version__ Out[2]: '2.0.0.dev-4daf949' In [3]: res = np.longdouble(2)**64 In [4]: res Out[4]: 18446744073709551616.0 In

Re: [Numpy-discussion] float128 / longdouble on PPC - is it broken?

2011-10-25 Thread Matthew Brett
Hi, On Tue, Oct 25, 2011 at 7:31 AM, Charles R Harris charlesr.har...@gmail.com wrote: On Mon, Oct 24, 2011 at 10:59 PM, Matthew Brett matthew.br...@gmail.com wrote: Hi, I just ran into this on a PPC machine: In [1]: import numpy as np In [2]: np.__version__ Out[2]:

Re: [Numpy-discussion] float128 / longdouble on PPC - is it broken?

2011-10-25 Thread Matthew Brett
Hi, On Tue, Oct 25, 2011 at 2:43 AM, Pauli Virtanen p...@iki.fi wrote: 25.10.2011 06:59, Matthew Brett kirjoitti: res = np.longdouble(2)**64 res-1 36893488147419103231.0 Can you check if long double works properly (not a given) in C on that platform:        long double x;        x =

Re: [Numpy-discussion] float128 / longdouble on PPC - is it broken?

2011-10-25 Thread Charles R Harris
On Tue, Oct 25, 2011 at 11:45 AM, Matthew Brett matthew.br...@gmail.comwrote: Hi, On Tue, Oct 25, 2011 at 2:43 AM, Pauli Virtanen p...@iki.fi wrote: 25.10.2011 06:59, Matthew Brett kirjoitti: res = np.longdouble(2)**64 res-1 36893488147419103231.0 Can you check if long double works

Re: [Numpy-discussion] float128 / longdouble on PPC - is it broken?

2011-10-25 Thread Matthew Brett
Hi, On Tue, Oct 25, 2011 at 10:52 AM, Charles R Harris charlesr.har...@gmail.com wrote: On Tue, Oct 25, 2011 at 11:45 AM, Matthew Brett matthew.br...@gmail.com wrote: Hi, On Tue, Oct 25, 2011 at 2:43 AM, Pauli Virtanen p...@iki.fi wrote: 25.10.2011 06:59, Matthew Brett kirjoitti: res

Re: [Numpy-discussion] float128 / longdouble on PPC - is it broken?

2011-10-25 Thread Matthew Brett
On Tue, Oct 25, 2011 at 11:05 AM, Matthew Brett matthew.br...@gmail.com wrote: Hi, On Tue, Oct 25, 2011 at 10:52 AM, Charles R Harris charlesr.har...@gmail.com wrote: On Tue, Oct 25, 2011 at 11:45 AM, Matthew Brett matthew.br...@gmail.com wrote: Hi, On Tue, Oct 25, 2011 at 2:43 AM,

Re: [Numpy-discussion] float128 / longdouble on PPC - is it broken?

2011-10-25 Thread Pauli Virtanen
25.10.2011 19:45, Matthew Brett kirjoitti: [clip] or, in case the platform doesn't have powl: long double x; x = pow(2, 64); x -= 1; printf(%g %Lg\n, (double)x, x); Both the same as numpy: [mb312@jerry ~]$ gcc test.c test.c: In function 'main': test.c:5:

Re: [Numpy-discussion] float128 / longdouble on PPC - is it broken?

2011-10-25 Thread Matthew Brett
Hi, On Tue, Oct 25, 2011 at 11:14 AM, Pauli Virtanen p...@iki.fi wrote: 25.10.2011 19:45, Matthew Brett kirjoitti: [clip] or, in case the platform doesn't have powl:         long double x;         x = pow(2, 64);         x -= 1;         printf(%g %Lg\n, (double)x, x); Both the same as

Re: [Numpy-discussion] float128 / longdouble on PPC - is it broken?

2011-10-25 Thread Matthew Brett
Hi, On Tue, Oct 25, 2011 at 12:01 PM, Derek Homeier de...@astro.physik.uni-goettingen.de wrote: On 25 Oct 2011, at 20:05, Matthew Brett wrote: Both the same as numpy: [mb312@jerry ~]$ gcc test.c test.c: In function 'main': test.c:5: warning: incompatible implicit declaration of built-in

Re: [Numpy-discussion] float128 / longdouble on PPC - is it broken?

2011-10-25 Thread Pauli Virtanen
25.10.2011 20:29, Matthew Brett kirjoitti: [clip] In [7]: (res-1) / 2**32 Out[7]: 8589934591.98 In [8]: np.float((res-1) / 2**32) Out[8]: 4294967296.0 Looks like a bug in the C library installed on the machine, then. It's either in wontfix territory for us, or in the cast to doubles

Re: [Numpy-discussion] float128 / longdouble on PPC - is it broken?

2011-10-25 Thread Derek Homeier
Hi, On 25 Oct 2011, at 21:14, Pauli Virtanen wrote: 25.10.2011 20:29, Matthew Brett kirjoitti: [clip] In [7]: (res-1) / 2**32 Out[7]: 8589934591.98 In [8]: np.float((res-1) / 2**32) Out[8]: 4294967296.0 Looks like a bug in the C library installed on the machine, then. It's

Re: [Numpy-discussion] float128 / longdouble on PPC - is it broken?

2011-10-25 Thread Matthew Brett
Hi, On Tue, Oct 25, 2011 at 12:14 PM, Pauli Virtanen p...@iki.fi wrote: 25.10.2011 20:29, Matthew Brett kirjoitti: [clip] In [7]: (res-1) / 2**32 Out[7]: 8589934591.98 In [8]: np.float((res-1) / 2**32) Out[8]: 4294967296.0 Looks like a bug in the C library installed on the

Re: [Numpy-discussion] float128 / longdouble on PPC - is it broken?

2011-10-25 Thread David Cournapeau
On Tue, Oct 25, 2011 at 8:22 PM, Matthew Brett matthew.br...@gmail.com wrote: Hi, On Tue, Oct 25, 2011 at 12:14 PM, Pauli Virtanen p...@iki.fi wrote: 25.10.2011 20:29, Matthew Brett kirjoitti: [clip] In [7]: (res-1) / 2**32 Out[7]: 8589934591.98 In [8]: np.float((res-1) / 2**32)

Re: [Numpy-discussion] float128 / longdouble on PPC - is it broken?

2011-10-25 Thread Matthew Brett
Hi, On Tue, Oct 25, 2011 at 2:58 PM, David Cournapeau courn...@gmail.com wrote: On Tue, Oct 25, 2011 at 8:22 PM, Matthew Brett matthew.br...@gmail.com wrote: Hi, On Tue, Oct 25, 2011 at 12:14 PM, Pauli Virtanen p...@iki.fi wrote: 25.10.2011 20:29, Matthew Brett kirjoitti: [clip] In [7]:

[Numpy-discussion] float128 / longdouble on PPC - is it broken?

2011-10-24 Thread Matthew Brett
Hi, I just ran into this on a PPC machine: In [1]: import numpy as np In [2]: np.__version__ Out[2]: '2.0.0.dev-4daf949' In [3]: res = np.longdouble(2)**64 In [4]: res Out[4]: 18446744073709551616.0 In [5]: 2**64 Out[5]: 18446744073709551616L In [6]: res-1 Out[6]: 36893488147419103231.0