[SciPy-user] floating point exeption

Steve Schmerler elcorto at gmx.net
Sun Dec 11 09:06:15 CST 2005


Andrew Straw wrote:
> Steve Schmerler wrote:
> 
> 
>>	scipy.test()
>>
>>I get a floating point exception:
>>
>>......................................................................................
>>.........................................................................................................................
>>.............................Gleitkomma-Ausnahme
>> 
>>
> 
> 
> Are you running Debian sarge on a machine with SSE?
> 
> GNU libc version 2.3.2 has a bug[1] "feclearexcept() error on CPUs with
> SSE" (fixed in 2.3.3) which has been  submitted to Debian[2] but not
> fixed in sarge. You can see my page about the issue, complete with patch[3].
> 
> [1] http://sources.redhat.com/bugzilla/show_bug.cgi?id=10
> [2] http://bugs.debian.org/cgi-bin/bugreport.cgi?bug=279294
> [3]
> http://www.its.caltech.edu/~astraw/coding.html#libc-patched-for-debian-sarge-to-fix-floating-point-exceptions-on-sse2
> 

As mentioned in a previous mail the exception occures when calling some 
fumctions in scipy.special._cephes.

I'm running Debian sarge stable. Some days ago I pretended that my CPU 
doesn't have SSE, but a look in /proc/cpuinfo reveals

-----------------------------------------------------------------------------------
processor       : 0
vendor_id       : AuthenticAMD
cpu family      : 6
model           : 8
model name      : AMD Athlon(TM) XP 2400+
stepping        : 1
cpu MHz         : 1990.656
cache size      : 256 KB
fdiv_bug        : no
hlt_bug         : no
f00f_bug        : no
coma_bug        : no
fpu             : yes
fpu_exception   : yes
cpuid level     : 1
wp              : yes
flags           : fpu vme de pse tsc msr pae mce cx8 apic sep mtrr pge 
mca cmov pat pse36 mmx fxsr sse syscall mmxext 3dnowext 3dnow
bogomips        : 3984.77
-----------------------------------------------------------------------------------

Is the patch [3] for SSE and/or SSE2?

I'm relatively new to Linux and don't really feel like fiddling with my 
libc :) Is it advisable to install libc6 from testing (version 2.3.5)?

cheers,
steve

-- 
Should array indices start at 0 or 1? My compromise of 0.5 was rejected 
without, I thought, proper consideration.
-- Stan Kelly-Bootle



More information about the SciPy-user mailing list