float exception, fblas cgemv - newscipy
Hi, on another debian sarge based 32Bit machine (where I made a couple of links so that ATLAS is found) I get the following check_default_beta_y (scipy.linalg.fblas.test_fblas.test_cgemv) ... zsh: 8527 floating point exception The shortest example to reproduce this is gdb file /usr/bin/python run from scipy import * from scipy.basic.random import normal a = normal(0.,1.,(3,3)) x = arange(shape(a)[0]) print a print x linalg.fblas.cgemv(1,a,x) Giving
print a [[-0.70685647 0.65494695 0.9179007 ] [-0.30862233 0.58137533 2.01614468] [ 0.18585164 -0.8134587 1.23613689]] print x [0 1 2] linalg.fblas.cgemv(1,a,x)
Program received signal SIGFPE, Arithmetic exception. [Switching to Thread 16384 (LWP 26791)] 0x404167a0 in ATL_cgemvC_a1_x1_bXi0_y1 () from /usr/lib/sse2/libatlas.so.3 Does anyone else see this? Any ideas what to do with this? Best, Arnd System details: ipython import scipy scipy.__core_version__ scipy.__scipy_version__ scipy.__scipy_config__.show() gives Python 2.3.5 (#2, Sep 4 2005, 22:01:42) In [1]:import scipy Importing io to scipy Importing fftpack to scipy Importing special to scipy Importing cluster to scipy Importing sparse to scipy Importing utils to scipy Importing interpolate to scipy Importing integrate to scipy Importing signal to scipy Importing optimize to scipy Importing linalg to scipy Importing stats to scipy In [2]:scipy.__core_version__ Out[2]:'0.4.3.1376' In [3]:scipy.__scipy_version__ Out[3]:'0.4.2_1400' In [4]:scipy.__scipy_config__.show() dfftw_info: NOT AVAILABLE blas_opt_info: libraries = ['f77blas', 'cblas', 'atlas'] library_dirs = ['/usr/lib/sse2'] define_macros = [('ATLAS_INFO', '"\\"3.6.0\\""')] language = c djbfft_info: NOT AVAILABLE atlas_blas_threads_info: NOT AVAILABLE lapack_opt_info: libraries = ['lapack', 'f77blas', 'cblas', 'atlas'] library_dirs = ['/usr/lib/atlas', '/usr/lib/'] define_macros = [('ATLAS_INFO', '"\\"3.6.0\\""')] language = f77 atlas_info: libraries = ['lapack', 'f77blas', 'cblas', 'atlas'] library_dirs = ['/usr/lib/atlas', '/usr/lib/'] language = f77 define_macros = [('ATLAS_INFO', '"\\"3.6.0\\""')] fftw_info: NOT AVAILABLE atlas_blas_info: libraries = ['f77blas', 'cblas', 'atlas'] library_dirs = ['/usr/lib/sse2'] language = c define_macros = [('ATLAS_INFO', '"\\"3.6.0\\""')] atlas_threads_info: NOT AVAILABLE Atlas: ls -l /usr/lib/atlas total 17168 lrwxrwxrwx 1 root root 12 Oct 25 00:32 libblas.so -> libblas.so.3 lrwxrwxrwx 1 root root 14 Mar 23 2005 libblas.so.2 -> libblas.so.2.3 -rw-r--r-- 1 root root 3051888 Mar 4 2005 libblas.so.2.3 lrwxrwxrwx 1 root root 14 Nov 18 2004 libblas.so.3 -> libblas.so.3.0 -rw-r--r-- 1 root root 3443424 Oct 27 2004 libblas.so.3.0 lrwxrwxrwx 1 root root 14 Oct 25 00:32 liblapack.so -> liblapack.so.3 lrwxrwxrwx 1 root root 16 Mar 23 2005 liblapack.so.2 -> liblapack.so.2.3 -rw-r--r-- 1 root root 5500688 Mar 4 2005 liblapack.so.2.3 lrwxrwxrwx 1 root root 16 Nov 18 2004 liblapack.so.3 -> liblapack.so.3.0 -rw-r--r-- 1 root root 5537840 Oct 27 2004 liblapack.so.3.0 drwxr-xr-x 2 root root 4096 Nov 18 2004 sse2 ls -l /usr/lib/ | grep atlas drwxr-xr-x 3 root root 4096 Oct 25 00:32 atlas lrwxrwxrwx 1 root root 16 Mar 22 2004 libatlas.so -> sse2/libatlas.so lrwxrwxrwx 1 root root 15 Mar 23 2005 libatlas.so.2 -> libatlas.so.2.3 -rw-r--r-- 1 root root 2835272 Mar 4 2005 libatlas.so.2.3 lrwxrwxrwx 1 root root 15 Nov 18 2004 libatlas.so.3 -> libatlas.so.3.0 -rw-r--r-- 1 root root 3234288 Oct 27 2004 libatlas.so.3.0 lrwxrwxrwx 1 root root 34 Mar 22 2004 liblapack.so -> /usr/lib/atlas/sse2/liblapack.so.3 lrwxrwxrwx 1 root root 25 Mar 22 2004 liblapack_atlas.so -> sse2/liblapack_atlas.so.3 lrwxrwxrwx 1 root root 22 Mar 23 2005 liblapack_atlas.so.2 -> liblapack_atlas.so.2.3 -rw-r--r-- 1 root root 60568 Mar 4 2005 liblapack_atlas.so.2.3 lrwxrwxrwx 1 root root 22 Nov 18 2004 liblapack_atlas.so.3 -> liblapack_atlas.so.3.0 -rw-r--r-- 1 root root 131224 Oct 27 2004 liblapack_atlas.so.3.0
Hi Arnd, Make sure you're not being bitten by this bug in glibc < 2.3.3: http://sources.redhat.com/bugzilla/show_bug.cgi?id=10 http://bugs.debian.org/cgi-bin/bugreport.cgi?bug=279294 Arnd Baecker wrote:
Hi,
on another debian sarge based 32Bit machine (where I made a couple of links so that ATLAS is found) I get the following
check_default_beta_y (scipy.linalg.fblas.test_fblas.test_cgemv) ... zsh: 8527 floating point exception
The shortest example to reproduce this is
gdb file /usr/bin/python run
from scipy import * from scipy.basic.random import normal a = normal(0.,1.,(3,3)) x = arange(shape(a)[0]) print a print x linalg.fblas.cgemv(1,a,x)
Giving
print a
[[-0.70685647 0.65494695 0.9179007 ] [-0.30862233 0.58137533 2.01614468] [ 0.18585164 -0.8134587 1.23613689]]
print x
[0 1 2]
linalg.fblas.cgemv(1,a,x)
Program received signal SIGFPE, Arithmetic exception. [Switching to Thread 16384 (LWP 26791)] 0x404167a0 in ATL_cgemvC_a1_x1_bXi0_y1 () from /usr/lib/sse2/libatlas.so.3
Does anyone else see this? Any ideas what to do with this?
Best, Arnd
System details:
ipython import scipy scipy.__core_version__ scipy.__scipy_version__ scipy.__scipy_config__.show()
gives
Python 2.3.5 (#2, Sep 4 2005, 22:01:42)
In [1]:import scipy Importing io to scipy Importing fftpack to scipy Importing special to scipy Importing cluster to scipy Importing sparse to scipy Importing utils to scipy Importing interpolate to scipy Importing integrate to scipy Importing signal to scipy Importing optimize to scipy Importing linalg to scipy Importing stats to scipy In [2]:scipy.__core_version__ Out[2]:'0.4.3.1376' In [3]:scipy.__scipy_version__ Out[3]:'0.4.2_1400' In [4]:scipy.__scipy_config__.show() dfftw_info: NOT AVAILABLE
blas_opt_info: libraries = ['f77blas', 'cblas', 'atlas'] library_dirs = ['/usr/lib/sse2'] define_macros = [('ATLAS_INFO', '"\\"3.6.0\\""')] language = c
djbfft_info: NOT AVAILABLE
atlas_blas_threads_info: NOT AVAILABLE
lapack_opt_info: libraries = ['lapack', 'f77blas', 'cblas', 'atlas'] library_dirs = ['/usr/lib/atlas', '/usr/lib/'] define_macros = [('ATLAS_INFO', '"\\"3.6.0\\""')] language = f77
atlas_info: libraries = ['lapack', 'f77blas', 'cblas', 'atlas'] library_dirs = ['/usr/lib/atlas', '/usr/lib/'] language = f77 define_macros = [('ATLAS_INFO', '"\\"3.6.0\\""')]
fftw_info: NOT AVAILABLE
atlas_blas_info: libraries = ['f77blas', 'cblas', 'atlas'] library_dirs = ['/usr/lib/sse2'] language = c define_macros = [('ATLAS_INFO', '"\\"3.6.0\\""')]
atlas_threads_info: NOT AVAILABLE
Atlas:
ls -l /usr/lib/atlas total 17168 lrwxrwxrwx 1 root root 12 Oct 25 00:32 libblas.so -> libblas.so.3 lrwxrwxrwx 1 root root 14 Mar 23 2005 libblas.so.2 -> libblas.so.2.3 -rw-r--r-- 1 root root 3051888 Mar 4 2005 libblas.so.2.3 lrwxrwxrwx 1 root root 14 Nov 18 2004 libblas.so.3 -> libblas.so.3.0 -rw-r--r-- 1 root root 3443424 Oct 27 2004 libblas.so.3.0 lrwxrwxrwx 1 root root 14 Oct 25 00:32 liblapack.so -> liblapack.so.3 lrwxrwxrwx 1 root root 16 Mar 23 2005 liblapack.so.2 -> liblapack.so.2.3 -rw-r--r-- 1 root root 5500688 Mar 4 2005 liblapack.so.2.3 lrwxrwxrwx 1 root root 16 Nov 18 2004 liblapack.so.3 -> liblapack.so.3.0 -rw-r--r-- 1 root root 5537840 Oct 27 2004 liblapack.so.3.0 drwxr-xr-x 2 root root 4096 Nov 18 2004 sse2
ls -l /usr/lib/ | grep atlas drwxr-xr-x 3 root root 4096 Oct 25 00:32 atlas lrwxrwxrwx 1 root root 16 Mar 22 2004 libatlas.so -> sse2/libatlas.so lrwxrwxrwx 1 root root 15 Mar 23 2005 libatlas.so.2 -> libatlas.so.2.3 -rw-r--r-- 1 root root 2835272 Mar 4 2005 libatlas.so.2.3 lrwxrwxrwx 1 root root 15 Nov 18 2004 libatlas.so.3 -> libatlas.so.3.0 -rw-r--r-- 1 root root 3234288 Oct 27 2004 libatlas.so.3.0 lrwxrwxrwx 1 root root 34 Mar 22 2004 liblapack.so -> /usr/lib/atlas/sse2/liblapack.so.3 lrwxrwxrwx 1 root root 25 Mar 22 2004 liblapack_atlas.so -> sse2/liblapack_atlas.so.3 lrwxrwxrwx 1 root root 22 Mar 23 2005 liblapack_atlas.so.2 -> liblapack_atlas.so.2.3 -rw-r--r-- 1 root root 60568 Mar 4 2005 liblapack_atlas.so.2.3 lrwxrwxrwx 1 root root 22 Nov 18 2004 liblapack_atlas.so.3 -> liblapack_atlas.so.3.0 -rw-r--r-- 1 root root 131224 Oct 27 2004 liblapack_atlas.so.3.0
_______________________________________________ Scipy-dev mailing list Scipy-dev@scipy.net http://www.scipy.net/mailman/listinfo/scipy-dev
participants (2)
-
Andrew Straw -
Arnd Baecker