[Scipy-tickets] [SciPy] #1636: Test Failures in scipy.linalg

SciPy Trac scipy-tickets@scipy....
Thu Mar 29 13:54:04 CDT 2012


#1636: Test Failures in scipy.linalg
--------------------------+-------------------------------------------------
 Reporter:  jseabold      |       Owner:  somebody   
     Type:  defect        |      Status:  new        
 Priority:  normal        |   Milestone:  Unscheduled
Component:  scipy.linalg  |     Version:  devel      
 Keywords:  qr, lapack    |  
--------------------------+-------------------------------------------------
 I noticed the following sign errors in TestQR. Is it true that these are
 only known up to a sign? I often see the "incorrect" sign versus examples
 elsewhere (and in Matlab). Is this possibly platform/LAPACK dependent?
 Should these tests get abs()?

 {{{
 [~/]
 [1]: from scipy import linalg

 [~/]
 [2]: linalg.test()
 Running unit tests for scipy.linalg
 NumPy version 1.7.0.dev-83480b6
 NumPy is installed in /usr/local/lib/python2.7/dist-packages/numpy
 SciPy version 0.11.0.dev-0fbfdbc
 SciPy is installed in /usr/local/lib/python2.7/dist-packages/scipy
 Python version 2.7.1+ (r271:86832, Apr 11 2011, 18:13:53) [GCC 4.5.2]
 nose version 1.0.0
 ...................................................................................................................................................F..F...F.F..F.....F......F..F....F.....F..F...F........................................................................................................................................................................K................................................................................................................................................................................................
 ======================================================================
 FAIL: test_random_complex_pivoting (test_decomp.TestQR)
 ----------------------------------------------------------------------
 Traceback (most recent call last):
   File "/usr/local/lib/python2.7/dist-
 packages/scipy/linalg/tests/test_decomp.py", line 1391, in
 test_random_complex_pivoting
     assert_array_almost_equal(q,q2)
   File "/usr/local/lib/python2.7/dist-packages/numpy/testing/utils.py",
 line 846, in assert_array_almost_equal
     header=('Arrays are not almost equal to %d decimals' % decimal))
   File "/usr/local/lib/python2.7/dist-packages/numpy/testing/utils.py",
 line 677, in assert_array_compare
     raise AssertionError(msg)
 AssertionError:
 Arrays are not almost equal to 6 decimals

 (mismatch 55.0%)
  x: array([[-0.00331202-0.17922166j,  0.09483222+0.07983082j,
         -0.27964260+0.09344627j,  0.00465585-0.14319313j,
         -0.15392783-0.24820327j, -0.25569917-0.1927968j ,...
  y: array([[ 0.00331202+0.17922166j,  0.09483222+0.07983082j,
          0.27964260-0.09344627j,  0.00465585-0.14319313j,
          0.15392783+0.24820327j,  0.25569917+0.1927968j ,...

 ======================================================================
 FAIL: test_random_pivoting (test_decomp.TestQR)
 ----------------------------------------------------------------------
 Traceback (most recent call last):
   File "/usr/local/lib/python2.7/dist-
 packages/scipy/linalg/tests/test_decomp.py", line 1245, in
 test_random_pivoting
     assert_array_almost_equal(q,q2)
   File "/usr/local/lib/python2.7/dist-packages/numpy/testing/utils.py",
 line 846, in assert_array_almost_equal
     header=('Arrays are not almost equal to %d decimals' % decimal))
   File "/usr/local/lib/python2.7/dist-packages/numpy/testing/utils.py",
 line 677, in assert_array_compare
     raise AssertionError(msg)
 AssertionError:
 Arrays are not almost equal to 6 decimals

 (mismatch 50.0%)
  x: array([[-0.0043425 ,  0.33678706,  0.20818063,  0.44351143,
 0.26530978,
          0.11052461, -0.26568683,  0.05482736,  0.00597287,  0.32064401,
         -0.03433912, -0.04185174,  0.27262429,  0.03804727,
 -0.01196234,...
  y: array([[ 0.0043425 ,  0.33678706,  0.20818063,  0.44351143,
 0.26530978,
          0.11052461,  0.26568683, -0.05482736,  0.00597287,  0.32064401,
          0.03433912,  0.04185174,  0.27262429,  0.03804727,
 0.01196234,...

 ======================================================================
 FAIL: test_random_tall_e_pivoting (test_decomp.TestQR)
 ----------------------------------------------------------------------
 Traceback (most recent call last):
   File "/usr/local/lib/python2.7/dist-
 packages/scipy/linalg/tests/test_decomp.py", line 1325, in
 test_random_tall_e_pivoting
     assert_array_almost_equal(q,q2)
   File "/usr/local/lib/python2.7/dist-packages/numpy/testing/utils.py",
 line 846, in assert_array_almost_equal
     header=('Arrays are not almost equal to %d decimals' % decimal))
   File "/usr/local/lib/python2.7/dist-packages/numpy/testing/utils.py",
 line 677, in assert_array_compare
     raise AssertionError(msg)
 AssertionError:
 Arrays are not almost equal to 6 decimals

 (mismatch 52.0%)
  x: array([[-0.03448532, -0.09992597, -0.04543201, ..., -0.04599707,
          0.04235504,  0.05398889],
        [-0.08044947, -0.04159241,  0.06052022, ..., -0.00430117,...
  y: array([[ 0.03448532,  0.09992597, -0.04543201, ..., -0.04599707,
          0.04235504, -0.05398889],
        [ 0.08044947,  0.04159241,  0.06052022, ..., -0.00430117,...

 ======================================================================
 FAIL: test_random_tall_pivoting (test_decomp.TestQR)
 ----------------------------------------------------------------------
 Traceback (most recent call last):
   File "/usr/local/lib/python2.7/dist-
 packages/scipy/linalg/tests/test_decomp.py", line 1296, in
 test_random_tall_pivoting
     assert_array_almost_equal(q,q2)
   File "/usr/local/lib/python2.7/dist-packages/numpy/testing/utils.py",
 line 846, in assert_array_almost_equal
     header=('Arrays are not almost equal to %d decimals' % decimal))
   File "/usr/local/lib/python2.7/dist-packages/numpy/testing/utils.py",
 line 677, in assert_array_compare
     raise AssertionError(msg)
 AssertionError:
 Arrays are not almost equal to 6 decimals

 (mismatch 75.9975%)
  x: array([[-0.03448532, -0.09992597, -0.04543201, ..., -0.00188773,
         -0.05446374, -0.02469976],
        [-0.08044947, -0.04159241,  0.06052022, ..., -0.0211871 ,...
  y: array([[ 0.03448532,  0.09992597, -0.04543201, ...,  0.08233543,
         -0.02840087,  0.10068053],
        [ 0.08044947,  0.04159241,  0.06052022, ...,  0.10416556,...

 ======================================================================
 FAIL: test_random_trap_pivoting (test_decomp.TestQR)
 ----------------------------------------------------------------------
 Traceback (most recent call last):
   File "/usr/local/lib/python2.7/dist-
 packages/scipy/linalg/tests/test_decomp.py", line 1348, in
 test_random_trap_pivoting
     assert_array_almost_equal(q,q2)
   File "/usr/local/lib/python2.7/dist-packages/numpy/testing/utils.py",
 line 846, in assert_array_almost_equal
     header=('Arrays are not almost equal to %d decimals' % decimal))
   File "/usr/local/lib/python2.7/dist-packages/numpy/testing/utils.py",
 line 677, in assert_array_compare
     raise AssertionError(msg)
 AssertionError:
 Arrays are not almost equal to 6 decimals

 (mismatch 52.0%)
  x: array([[-0.11200166, -0.10697248,  0.04754531, ...,  0.13145621,
          0.067521  , -0.06431018],
        [-0.04720016, -0.0936498 ,  0.10080133, ...,  0.040549  ,...
  y: array([[ 0.11200166,  0.10697248, -0.04754531, ..., -0.13145621,
         -0.067521  , -0.06431018],
        [ 0.04720016,  0.0936498 , -0.10080133, ..., -0.040549  ,...

 ======================================================================
 FAIL: test_simple_complex_pivoting (test_decomp.TestQR)
 ----------------------------------------------------------------------
 Traceback (most recent call last):
   File "/usr/local/lib/python2.7/dist-
 packages/scipy/linalg/tests/test_decomp.py", line 1188, in
 test_simple_complex_pivoting
     assert_array_almost_equal(q,q2)
   File "/usr/local/lib/python2.7/dist-packages/numpy/testing/utils.py",
 line 846, in assert_array_almost_equal
     header=('Arrays are not almost equal to %d decimals' % decimal))
   File "/usr/local/lib/python2.7/dist-packages/numpy/testing/utils.py",
 line 677, in assert_array_compare
     raise AssertionError(msg)
 AssertionError:
 Arrays are not almost equal to 6 decimals

 (mismatch 100.0%)
  x: array([[-0.44367825+0.j        , -0.34615118-0.76601974j,
         -0.01408746+0.31039373j],
        [-0.17747130-0.62114956j, -0.36858690+0.38942008j,...
  y: array([[ 0.44367825+0.j        ,  0.34615118+0.76601974j,
          0.01408746-0.31039373j],
        [ 0.17747130+0.62114956j,  0.36858690-0.38942008j,...

 ======================================================================
 FAIL: test_simple_fat_e_pivoting (test_decomp.TestQR)
 ----------------------------------------------------------------------
 Traceback (most recent call last):
   File "/usr/local/lib/python2.7/dist-
 packages/scipy/linalg/tests/test_decomp.py", line 1081, in
 test_simple_fat_e_pivoting
     assert_array_almost_equal(q,q2)
   File "/usr/local/lib/python2.7/dist-packages/numpy/testing/utils.py",
 line 846, in assert_array_almost_equal
     header=('Arrays are not almost equal to %d decimals' % decimal))
   File "/usr/local/lib/python2.7/dist-packages/numpy/testing/utils.py",
 line 677, in assert_array_compare
     raise AssertionError(msg)
 AssertionError:
 Arrays are not almost equal to 6 decimals

 (mismatch 100.0%)
  x: array([[-0.21693046, -0.97618706],
        [-0.97618706,  0.21693046]])
  y: array([[ 0.21693046,  0.97618706],
        [ 0.97618706, -0.21693046]])

 ======================================================================
 FAIL: test_simple_fat_pivoting (test_decomp.TestQR)
 ----------------------------------------------------------------------
 Traceback (most recent call last):
   File "/usr/local/lib/python2.7/dist-
 packages/scipy/linalg/tests/test_decomp.py", line 1058, in
 test_simple_fat_pivoting
     assert_array_almost_equal(q,q2)
   File "/usr/local/lib/python2.7/dist-packages/numpy/testing/utils.py",
 line 846, in assert_array_almost_equal
     header=('Arrays are not almost equal to %d decimals' % decimal))
   File "/usr/local/lib/python2.7/dist-packages/numpy/testing/utils.py",
 line 677, in assert_array_compare
     raise AssertionError(msg)
 AssertionError:
 Arrays are not almost equal to 6 decimals

 (mismatch 100.0%)
  x: array([[-0.21693046, -0.97618706],
        [-0.97618706,  0.21693046]])
  y: array([[ 0.21693046,  0.97618706],
        [ 0.97618706, -0.21693046]])

 ======================================================================
 FAIL: test_simple_pivoting (test_decomp.TestQR)
 ----------------------------------------------------------------------
 Traceback (most recent call last):
   File "/usr/local/lib/python2.7/dist-
 packages/scipy/linalg/tests/test_decomp.py", line 922, in
 test_simple_pivoting
     assert_array_almost_equal(q,q2)
   File "/usr/local/lib/python2.7/dist-packages/numpy/testing/utils.py",
 line 846, in assert_array_almost_equal
     header=('Arrays are not almost equal to %d decimals' % decimal))
   File "/usr/local/lib/python2.7/dist-packages/numpy/testing/utils.py",
 line 677, in assert_array_compare
     raise AssertionError(msg)
 AssertionError:
 Arrays are not almost equal to 6 decimals

 (mismatch 33.3333333333%)
  x: array([[-0.20628425,  0.84710051, -0.48976273],
        [-0.92827912, -0.32770096, -0.17581226],
        [-0.30942637,  0.41836921,  0.85394527]])
  y: array([[ 0.20628425,  0.84710051, -0.48976273],
        [ 0.92827912, -0.32770096, -0.17581226],
        [ 0.30942637,  0.41836921,  0.85394527]])

 ======================================================================
 FAIL: test_simple_tall_e_pivoting (test_decomp.TestQR)
 ----------------------------------------------------------------------
 Traceback (most recent call last):
   File "/usr/local/lib/python2.7/dist-
 packages/scipy/linalg/tests/test_decomp.py", line 993, in
 test_simple_tall_e_pivoting
     assert_array_almost_equal(q,q2)
   File "/usr/local/lib/python2.7/dist-packages/numpy/testing/utils.py",
 line 846, in assert_array_almost_equal
     header=('Arrays are not almost equal to %d decimals' % decimal))
   File "/usr/local/lib/python2.7/dist-packages/numpy/testing/utils.py",
 line 677, in assert_array_compare
     raise AssertionError(msg)
 AssertionError:
 Arrays are not almost equal to 6 decimals

 (mismatch 50.0%)
  x: array([[-0.20628425,  0.84710051],
        [-0.92827912, -0.32770096],
        [-0.30942637,  0.41836921]])
  y: array([[ 0.20628425,  0.84710051],
        [ 0.92827912, -0.32770096],
        [ 0.30942637,  0.41836921]])

 ======================================================================
 FAIL: test_simple_tall_pivoting (test_decomp.TestQR)
 ----------------------------------------------------------------------
 Traceback (most recent call last):
   File "/usr/local/lib/python2.7/dist-
 packages/scipy/linalg/tests/test_decomp.py", line 972, in
 test_simple_tall_pivoting
     assert_array_almost_equal(q,q2)
   File "/usr/local/lib/python2.7/dist-packages/numpy/testing/utils.py",
 line 846, in assert_array_almost_equal
     header=('Arrays are not almost equal to %d decimals' % decimal))
   File "/usr/local/lib/python2.7/dist-packages/numpy/testing/utils.py",
 line 677, in assert_array_compare
     raise AssertionError(msg)
 AssertionError:
 Arrays are not almost equal to 6 decimals

 (mismatch 66.6666666667%)
  x: array([[-0.20628425,  0.84710051, -0.48976273],
        [-0.92827912, -0.32770096, -0.17581226],
        [-0.30942637,  0.41836921,  0.85394527]])
  y: array([[ 0.20628425,  0.84710051,  0.48976273],
        [ 0.92827912, -0.32770096,  0.17581226],
        [ 0.30942637,  0.41836921, -0.85394527]])

 ======================================================================
 FAIL: test_simple_trap_pivoting (test_decomp.TestQR)
 ----------------------------------------------------------------------
 Traceback (most recent call last):
   File "/usr/local/lib/python2.7/dist-
 packages/scipy/linalg/tests/test_decomp.py", line 953, in
 test_simple_trap_pivoting
     assert_array_almost_equal(q,q2)
   File "/usr/local/lib/python2.7/dist-packages/numpy/testing/utils.py",
 line 846, in assert_array_almost_equal
     header=('Arrays are not almost equal to %d decimals' % decimal))
   File "/usr/local/lib/python2.7/dist-packages/numpy/testing/utils.py",
 line 677, in assert_array_compare
     raise AssertionError(msg)
 AssertionError:
 Arrays are not almost equal to 6 decimals

 (mismatch 100.0%)
  x: array([[-0.21693046, -0.97618706],
        [-0.97618706,  0.21693046]])
  y: array([[ 0.21693046,  0.97618706],
        [ 0.97618706, -0.21693046]])

 ----------------------------------------------------------------------
 Ran 555 tests in 2.460s

 FAILED (KNOWNFAIL=1, failures=12)
 [2]: <nose.result.TextTestResult run=555 errors=0 failures=12>

 [~/]
 [3]: np.show_config()
 atlas_threads_info:
     libraries = ['lapack', 'ptf77blas', 'ptcblas', 'atlas']
     library_dirs = ['/home/skipper/built/atlas_files/lib']
     define_macros = [('ATLAS_INFO', '"\\"3.9.32\\""')]
     language = f77
     include_dirs = ['/home/skipper/built/atlas_files/include']
 blas_opt_info:
     libraries = ['ptf77blas', 'ptcblas', 'atlas']
     library_dirs = ['/home/skipper/built/atlas_files/lib']
     define_macros = [('ATLAS_INFO', '"\\"3.9.32\\""')]
     language = c
     include_dirs = ['/home/skipper/built/atlas_files/include']
 atlas_blas_threads_info:
     libraries = ['ptf77blas', 'ptcblas', 'atlas']
     library_dirs = ['/home/skipper/built/atlas_files/lib']
     define_macros = [('ATLAS_INFO', '"\\"3.9.32\\""')]
     language = c
     include_dirs = ['/home/skipper/built/atlas_files/include']
 lapack_opt_info:
     libraries = ['lapack', 'ptf77blas', 'ptcblas', 'atlas']
     library_dirs = ['/home/skipper/built/atlas_files/lib']
     define_macros = [('ATLAS_INFO', '"\\"3.9.32\\""')]
     language = f77
     include_dirs = ['/home/skipper/built/atlas_files/include']
 lapack_mkl_info:
   NOT AVAILABLE
 blas_mkl_info:
   NOT AVAILABLE
 mkl_info:
   NOT AVAILABLE

 [~/]
 [4]: from scipy import version

 [~/]
 [5]: version.full_version
 [5]: '0.11.0.dev-0fbfdbc'
 }}}

-- 
Ticket URL: <http://projects.scipy.org/scipy/ticket/1636>
SciPy <http://www.scipy.org>
SciPy is open-source software for mathematics, science, and engineering.


More information about the Scipy-tickets mailing list