[SciPy-dev] LAPACK issue

[Nils Wagner]

Hi all,

scipy.test() has some failures which seem to be related to the new
LAPACK version
http://www.netlib.org/lapack/lapack-3.1.1.changes


======================================================================
FAIL: check_syevr (scipy.lib.tests.test_lapack.test_flapack_float)
----------------------------------------------------------------------
Traceback (most recent call last):
  File
"/usr/local/lib64/python2.5/site-packages/scipy/lib/lapack/tests/esv_tests.py",
line 41, in check_syevr
    assert_array_almost_equal(w,exact_w)
  File
"/usr/local/lib64/python2.5/site-packages/numpy/testing/utils.py", line
230, in assert_array_almost_equal
    header='Arrays are not almost equal')
  File
"/usr/local/lib64/python2.5/site-packages/numpy/testing/utils.py", line
215, in assert_array_compare
    assert cond, msg
AssertionError:
Arrays are not almost equal

(mismatch 33.3333333333%)
 x: array([-0.66992444,  0.48769474,  9.18222618], dtype=float32)
 y: array([-0.66992434,  0.48769389,  9.18223045])

======================================================================
FAIL: check_syevr_irange (scipy.lib.tests.test_lapack.test_flapack_float)
----------------------------------------------------------------------
Traceback (most recent call last):
  File
"/usr/local/lib64/python2.5/site-packages/scipy/lib/lapack/tests/esv_tests.py",
line 66, in check_syevr_irange
    assert_array_almost_equal(w,exact_w[rslice])
  File
"/usr/local/lib64/python2.5/site-packages/numpy/testing/utils.py", line
230, in assert_array_almost_equal
    header='Arrays are not almost equal')
  File
"/usr/local/lib64/python2.5/site-packages/numpy/testing/utils.py", line
215, in assert_array_compare
    assert cond, msg
AssertionError:
Arrays are not almost equal

(mismatch 33.3333333333%)
 x: array([-0.66992444,  0.48769474,  9.18222618], dtype=float32)
 y: array([-0.66992434,  0.48769389,  9.18223045])

======================================================================
FAIL: Solve: single precision complex
----------------------------------------------------------------------
Traceback (most recent call last):
  File
"/usr/local/lib64/python2.5/site-packages/scipy/linsolve/umfpack/tests/test_umfpack.py",
line 32, in check_solve_complex_without_umfpack
    assert_array_almost_equal(a*x, b)
  File
"/usr/local/lib64/python2.5/site-packages/numpy/testing/utils.py", line
230, in assert_array_almost_equal
    header='Arrays are not almost equal')
  File
"/usr/local/lib64/python2.5/site-packages/numpy/testing/utils.py", line
215, in assert_array_compare
    assert cond, msg
AssertionError:
Arrays are not almost equal

(mismatch 20.0%)
 x: array([ 1.        +0.j,  1.99999809+0.j,  3.        +0.j, 
4.00000048+0.j,
        5.        +0.j], dtype=complex64)
 y: array([1, 2, 3, 4, 5])

----------------------------------------------------------------------

These failures are not present in case of LAPACK 3.0.

Any comment ?

Nils