Re: [Numpy-discussion] Strange error raised by scipy.special.erf

2012-01-24 Thread Pierre Haessig
Le 22/01/2012 11:28, Nadav Horesh a écrit :
  special.erf(26.5)
 1.0
  special.erf(26.6)
 Traceback (most recent call last):
   File pyshell#7, line 1, in module
 special.erf(26.6)
 FloatingPointError: underflow encountered in erf
  special.erf(26.7)
 1.0
  
I can confirm this same behaviour with numpy 1.5.1/scipy 0.9.0
Indeed 26.5 and 26.7 works, while 26.6 raises the underflow... weird
enough !
-- 
Pierre
___
NumPy-Discussion mailing list
NumPy-Discussion@scipy.org
http://mail.scipy.org/mailman/listinfo/numpy-discussion


Re: [Numpy-discussion] Strange error raised by scipy.special.erf

2012-01-24 Thread Nadav Horesh
I filed a ticket (#1590).

 Thank you for the verification.

   Nadav.

From: numpy-discussion-boun...@scipy.org [numpy-discussion-boun...@scipy.org] 
On Behalf Of Pierre Haessig [pierre.haes...@crans.org]
Sent: 24 January 2012 16:01
To: numpy-discussion@scipy.org
Subject: Re: [Numpy-discussion] Strange error raised by scipy.special.erf

Le 22/01/2012 11:28, Nadav Horesh a écrit :
  special.erf(26.5)
 1.0
  special.erf(26.6)
 Traceback (most recent call last):
   File pyshell#7, line 1, in module
 special.erf(26.6)
 FloatingPointError: underflow encountered in erf
  special.erf(26.7)
 1.0

I can confirm this same behaviour with numpy 1.5.1/scipy 0.9.0
Indeed 26.5 and 26.7 works, while 26.6 raises the underflow... weird
enough !
--
Pierre
___
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] Strange error raised by scipy.special.erf

2012-01-22 Thread Nadav Horesh
With N.seterr(all='raise'):

 from scipy import special
 import scipy
 special.erf(26.6)
1.0
 scipy.__version__
'0.11.0.dev-81dc505'
 import numpy as N
 N.seterr(all='raise')
{'over': 'warn', 'divide': 'warn', 'invalid': 'warn', 'under': 'ignore'}
 special.erf(26.5)
1.0
 special.erf(26.6)
Traceback (most recent call last):
  File pyshell#7, line 1, in module
special.erf(26.6)
FloatingPointError: underflow encountered in erf
 special.erf(26.7)
1.0

What is so special in 26.6?
I have this error also with previous versions of scipy

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