Re: [Numpy-discussion] Warnings not raised by np.log in 32 bit build on Windows

2013-08-23 Thread Charles R Harris
On Fri, Aug 23, 2013 at 6:29 AM, Nathaniel Smith  wrote:

> Probably the thing to do for reliable behaviour is to decide on the
> behaviour we want and then implement it "by hand". I.e., either clear the
> FP flags inside the ufunc loop (if we decide that log shouldn't raise a
> warning), or else check for nan and set the invalid flag ourselves.
> (Checking for nan should be much cheaper than computing log, I think, so
> this should be okay speed-wise?)
> On 23 Aug 2013 13:16, "Charles R Harris" 
> wrote:
>
>> These things may depend on how the compiler implements various calls.
>> Some errors went the other way with Julian's SIMD work, i.e., errors
>> getting set that were not set before. I'm not sure what can be done about
>> it.
>>
>>
>> On Thu, Aug 22, 2013 at 8:32 PM, Warren Weckesser <
>> warren.weckes...@gmail.com> wrote:
>>
>>> I'm investigating a test error in scipy 0.13.0 beta 1 that was
>>> reported by Christoph Gohlke.  The scipy issue is here:
>>> https://github.com/scipy/scipy/issues/2771
>>>
>>> I don't have a Windows environment to test it myself, but Christoph
>>> reported that this code:
>>>
>>> ```
>>> import numpy as np
>>>
>>> data = np.array([-0.375, -0.25, 0.0])
>>> s = np.log(data)
>>> ```
>>>
>>> does not generate two RuntimeWarnings when it is run with numpy 1.7.1
>>> in a 32 bit Windows 8 environment (numpy 1.7.1 compiled with Visual
>>> Studio compilers and Intel's MKL).  In 64 bit Windows, and in 64 bit
>>> linux, it generates two RuntimeWarnings.
>>>
>>> The inconsistency seems like a bug, possibly this one:
>>> https://github.com/numpy/numpy/issues/1958.
>>>
>>> Can anyone check if this also occurs in the development branch?
>>>
>>>
ISTR that is can also depend on instruction reordering in the hardware.
There was a bug opened on gcc for this a couple of years ago, but I did not
have the impression that it was going to get fixed any time soon.
Explicitly checking may add noticeable overhead and probably needs to be
profiled if we go that way.

Chuck
___
NumPy-Discussion mailing list
NumPy-Discussion@scipy.org
http://mail.scipy.org/mailman/listinfo/numpy-discussion


Re: [Numpy-discussion] Warnings not raised by np.log in 32 bit build on Windows

2013-08-23 Thread Nathaniel Smith
Probably the thing to do for reliable behaviour is to decide on the
behaviour we want and then implement it "by hand". I.e., either clear the
FP flags inside the ufunc loop (if we decide that log shouldn't raise a
warning), or else check for nan and set the invalid flag ourselves.
(Checking for nan should be much cheaper than computing log, I think, so
this should be okay speed-wise?)
On 23 Aug 2013 13:16, "Charles R Harris"  wrote:

> These things may depend on how the compiler implements various calls. Some
> errors went the other way with Julian's SIMD work, i.e., errors getting set
> that were not set before. I'm not sure what can be done about it.
>
>
> On Thu, Aug 22, 2013 at 8:32 PM, Warren Weckesser <
> warren.weckes...@gmail.com> wrote:
>
>> I'm investigating a test error in scipy 0.13.0 beta 1 that was
>> reported by Christoph Gohlke.  The scipy issue is here:
>> https://github.com/scipy/scipy/issues/2771
>>
>> I don't have a Windows environment to test it myself, but Christoph
>> reported that this code:
>>
>> ```
>> import numpy as np
>>
>> data = np.array([-0.375, -0.25, 0.0])
>> s = np.log(data)
>> ```
>>
>> does not generate two RuntimeWarnings when it is run with numpy 1.7.1
>> in a 32 bit Windows 8 environment (numpy 1.7.1 compiled with Visual
>> Studio compilers and Intel's MKL).  In 64 bit Windows, and in 64 bit
>> linux, it generates two RuntimeWarnings.
>>
>> The inconsistency seems like a bug, possibly this one:
>> https://github.com/numpy/numpy/issues/1958.
>>
>> Can anyone check if this also occurs in the development branch?
>>
>> Warren
>> ___
>> NumPy-Discussion mailing list
>> NumPy-Discussion@scipy.org
>> http://mail.scipy.org/mailman/listinfo/numpy-discussion
>>
>
>
> ___
> NumPy-Discussion mailing list
> NumPy-Discussion@scipy.org
> http://mail.scipy.org/mailman/listinfo/numpy-discussion
>
>
___
NumPy-Discussion mailing list
NumPy-Discussion@scipy.org
http://mail.scipy.org/mailman/listinfo/numpy-discussion


Re: [Numpy-discussion] Warnings not raised by np.log in 32 bit build on Windows

2013-08-23 Thread Alan G Isaac
On 8/22/2013 10:32 PM, Warren Weckesser wrote:
> Christoph
> reported that this code:
>
> ```
> import numpy as np
>
> data = np.array([-0.375, -0.25, 0.0])
> s = np.log(data)
> ```
>
> does not generate two RuntimeWarnings when it is run with numpy 1.7.1
> in a 32 bit Windows 8 environment (numpy 1.7.1 compiled with Visual
> Studio compilers and Intel's MKL).


Not sure if you want other (not Win 8) reports related to this,
but ... I'm appending (no runtime errors) below.  The OS is
Win 7 (64bit).

Alan Isaac


Enthought Python Distribution -- www.enthought.com
Version: 7.3-2 (32-bit)

Python 2.7.3 |EPD 7.3-2 (32-bit)| (default, Apr 12 2012, 14:30:37) [MSC v.1500 
3 2 bit (Intel)] on win32
Type "credits", "demo" or "enthought" for more information.
 >>> import numpy as np
 >>> np.__version__
'1.7.1'
 >>> data = np.array([-0.375, -0.25, 0.0])
 >>> s = np.log(data)
 >>>


___
NumPy-Discussion mailing list
NumPy-Discussion@scipy.org
http://mail.scipy.org/mailman/listinfo/numpy-discussion


Re: [Numpy-discussion] Warnings not raised by np.log in 32 bit build on Windows

2013-08-23 Thread Charles R Harris
These things may depend on how the compiler implements various calls. Some
errors went the other way with Julian's SIMD work, i.e., errors getting set
that were not set before. I'm not sure what can be done about it.


On Thu, Aug 22, 2013 at 8:32 PM, Warren Weckesser <
warren.weckes...@gmail.com> wrote:

> I'm investigating a test error in scipy 0.13.0 beta 1 that was
> reported by Christoph Gohlke.  The scipy issue is here:
> https://github.com/scipy/scipy/issues/2771
>
> I don't have a Windows environment to test it myself, but Christoph
> reported that this code:
>
> ```
> import numpy as np
>
> data = np.array([-0.375, -0.25, 0.0])
> s = np.log(data)
> ```
>
> does not generate two RuntimeWarnings when it is run with numpy 1.7.1
> in a 32 bit Windows 8 environment (numpy 1.7.1 compiled with Visual
> Studio compilers and Intel's MKL).  In 64 bit Windows, and in 64 bit
> linux, it generates two RuntimeWarnings.
>
> The inconsistency seems like a bug, possibly this one:
> https://github.com/numpy/numpy/issues/1958.
>
> Can anyone check if this also occurs in the development branch?
>
> Warren
> ___
> NumPy-Discussion mailing list
> NumPy-Discussion@scipy.org
> http://mail.scipy.org/mailman/listinfo/numpy-discussion
>
___
NumPy-Discussion mailing list
NumPy-Discussion@scipy.org
http://mail.scipy.org/mailman/listinfo/numpy-discussion