[SciPy-dev] float exception, fblas cgemv - newscipy

Andrew Straw strawman at astraw.com
Tue Nov 1 13:04:38 CST 2005


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 at scipy.net
> http://www.scipy.net/mailman/listinfo/scipy-dev




More information about the Scipy-dev mailing list