Strange error raised by scipy.special.erf

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.

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

I filed a ticket (#1590). Thank you for the verification. Nadav. ________________________________________ From: numpy-discussion-bounces@scipy.org [numpy-discussion-bounces@scipy.org] On Behalf Of Pierre Haessig [pierre.haessig@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
participants (2)
-
Nadav Horesh
-
Pierre Haessig