[Scipy-tickets] [SciPy] #1515: Arpack complex_nonsymmetric_modes failure on 64-bit Linux

SciPy Trac scipy-tickets@scipy....
Mon Sep 26 09:58:50 CDT 2011


#1515: Arpack complex_nonsymmetric_modes failure on 64-bit Linux
---------------------------------+------------------------------------------
 Reporter:  rgommers             |       Owner:  wnbell
     Type:  defect               |      Status:  new   
 Priority:  normal               |   Milestone:  0.10.0
Component:  scipy.sparse.linalg  |     Version:  0.9.0 
 Keywords:                       |  
---------------------------------+------------------------------------------

Comment(by bsouthey):

 I do not get that error with Fedora provided Atlas 3.8.3 installed but
 other similar errors. But if I build numpy (dev version) and scipy using
 "ATLAS=None python setup.py install" I do not get any arpack errors. So I
 think this is an Atlas related issue but I do not have Atlas 3.8.2
 available.

 I do get an error for "test_arpack.test_symmetric_modes" under Python2.5,
 2.6 and 2.7. Curiously I also get an error for
 "test_arpack.test_hermitian_modes" using Fedora 14 host with Python2.7 and
 numpy dev 2.0.0.dev-aded70c but not on another Fedora 15 host with
 Python2.7 numpy 1.6.1rc2.

 {{{
 $ python -c "import scipy; scipy.test(); scipy.show_config()"
 Running unit tests for scipy
 NumPy version 2.0.0.dev-aded70c
 NumPy is installed in /usr/lib64/python2.7/site-packages/numpy
 SciPy version 0.10.0b2
 SciPy is installed in /usr/lib64/python2.7/site-packages/scipy
 Python version 2.7 (r27:82500, Sep 16 2010, 18:02:00) [GCC 4.5.1 20100907
 (Red Hat 4.5.1-3)]
 [snip]

 ======================================================================
 ERROR: test_arpack.test_hermitian_modes(True, <gen-hermitian>, 'F', 2,
 'SA', None, None, <function aslinearoperator at 0x16c32a8>)
 ----------------------------------------------------------------------
 Traceback (most recent call last):
   File "/usr/lib/python2.7/site-
 packages/nose-1.0.0-py2.7.egg/nose/case.py", line 187, in runTest
     self.test(*self.arg)
   File "/usr/lib64/python2.7/site-
 packages/scipy/sparse/linalg/eigen/arpack/tests/test_arpack.py", line 168,
 in eval_evec
     eval, evec = eigs_func(ac, k, bc, **kwargs)
   File "/usr/lib64/python2.7/site-
 packages/scipy/sparse/linalg/eigen/arpack/arpack.py", line 1425, in eigsh
     OPinv=OPinv)
   File "/usr/lib64/python2.7/site-
 packages/scipy/sparse/linalg/eigen/arpack/arpack.py", line 1251, in eigs
     params.iterate()
   File "/usr/lib64/python2.7/site-
 packages/scipy/sparse/linalg/eigen/arpack/arpack.py", line 727, in iterate
     self.workd[yslice] = self.OP(self.workd[xslice])
   File "/usr/lib64/python2.7/site-
 packages/scipy/sparse/linalg/eigen/arpack/arpack.py", line 641, in
 <lambda>
     self.OP = lambda x: Minv_matvec(matvec(x))
   File "/usr/lib64/python2.7/site-
 packages/scipy/sparse/linalg/interface.py", line 123, in matvec
     y = self._matvec(x)
   File "/usr/lib64/python2.7/site-
 packages/scipy/sparse/linalg/eigen/arpack/arpack.py", line 945, in _matvec
     % (self.ifunc.__name__, info))
 ValueError: Error in inverting M: function gmres did not converge (info =
 60).

 ======================================================================
 FAIL: test_arpack.test_symmetric_modes(True, <gen-symmetric>, 'f', 2,
 'SM', None, 0.5, <function aslinearoperator at 0x16c32a8>, None, 'normal')
 ----------------------------------------------------------------------
 Traceback (most recent call last):
   File "/usr/lib/python2.7/site-
 packages/nose-1.0.0-py2.7.egg/nose/case.py", line 187, in runTest
     self.test(*self.arg)
   File "/usr/lib64/python2.7/site-
 packages/scipy/sparse/linalg/eigen/arpack/tests/test_arpack.py", line 187,
 in eval_evec
     err_msg=err)
   File "/usr/lib64/python2.7/site-
 packages/scipy/sparse/linalg/eigen/arpack/tests/test_arpack.py", line 74,
 in assert_allclose_cc
     assert_allclose(actual, conj(desired), **kw)
   File "/usr/lib64/python2.7/site-packages/numpy/testing/utils.py", line
 1213, in assert_allclose
     verbose=verbose, header=header)
   File "/usr/lib64/python2.7/site-packages/numpy/testing/utils.py", line
 677, in assert_array_compare
     raise AssertionError(msg)
 AssertionError:
 Not equal to tolerance rtol=0.000238419, atol=0.000238419
 error for eigsh:general, typ=f, which=SM, sigma=0.5,
 mattype=aslinearoperator, OPpart=None, mode=normal
 (mismatch 100.0%)
  x: array([ 34.53577423,   5.23088455], dtype=float32)
  y: array([ 34.52712250-0.j,   5.23087454-0.j], dtype=complex64)
 }}}

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


More information about the Scipy-tickets mailing list