[SciPy-Dev] ERROR: test_arpack.test_hermitian_modes

[Nils Wagner]

Hi Pauli,

Thank you for very much for the extension to generalized 
eigenvalue problems.
Your recent changes wrt arpack generate an error/failure.

Nils

======================================================================
ERROR: test_arpack.test_hermitian_modes(True, 
<gen-hermitian>, 'F', 2, 'SM', None, None, <function 
aslinearoperator at 0x2640a28>)
----------------------------------------------------------------------
Traceback (most recent call last):
   File 
"/home/nwagner/local/lib64/python2.6/site-packages/nose-0.11.2.dev-py2.6.egg/nose/case.py", 
line 183, in runTest
     self.test(*self.arg)
   File 
"/home/nwagner/local/lib64/python2.6/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 
"/home/nwagner/local/lib64/python2.6/site-packages/scipy/sparse/linalg/eigen/arpack/arpack.py", 
line 1425, in eigsh
     OPinv=OPinv)
   File 
"/home/nwagner/local/lib64/python2.6/site-packages/scipy/sparse/linalg/eigen/arpack/arpack.py", 
line 1251, in eigs
     params.iterate()
   File 
"/home/nwagner/local/lib64/python2.6/site-packages/scipy/sparse/linalg/eigen/arpack/arpack.py", 
line 727, in iterate
     self.workd[yslice] = self.OP(self.workd[xslice])
   File 
"/home/nwagner/local/lib64/python2.6/site-packages/scipy/sparse/linalg/eigen/arpack/arpack.py", 
line 641, in <lambda>
     self.OP = lambda x: Minv_matvec(matvec(x))
   File 
"/home/nwagner/local/lib64/python2.6/site-packages/scipy/sparse/linalg/interface.py", 
line 123, in matvec
     y = self._matvec(x)
   File 
"/home/nwagner/local/lib64/python2.6/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_hermitian_modes(True, 
<gen-hermitian>, 'D', 2, 'LM', None, None, <function 
aslinearoperator at 0x2640a28>)
----------------------------------------------------------------------
Traceback (most recent call last):
   File 
"/home/nwagner/local/lib64/python2.6/site-packages/nose-0.11.2.dev-py2.6.egg/nose/case.py", 
line 183, in runTest
     self.test(*self.arg)
   File 
"/home/nwagner/local/lib64/python2.6/site-packages/scipy/sparse/linalg/eigen/arpack/tests/test_arpack.py", 
line 196, in eval_evec
     assert_allclose(LHS, RHS, rtol=_rtol[typ], 
atol=_atol[typ], err_msg=err)
   File 
"/home/nwagner/local/lib64/python2.6/site-packages/numpy/testing/utils.py", 
line 1167, in assert_allclose
     verbose=verbose, header=header)
   File 
"/home/nwagner/local/lib64/python2.6/site-packages/numpy/testing/utils.py", 
line 635, in assert_array_compare
     raise AssertionError(msg)
AssertionError:
Not equal to tolerance rtol=2.22045e-13, atol=2.22045e-13
error for eigsh:general, typ=D, which=LM, sigma=None, 
mattype=aslinearoperator, OPpart=None, mode=normal
(mismatch 100.0%)
  x: array([[ 0.40172292+3.97785014j, 
 1.26541377-2.2227806j ],
        [ 0.67562574+2.73428778j, 
-6.32196547-3.30680082j],
        [-0.75603444+1.17747523j, 
-0.02785408-1.70286801j],...
  y: array([[ 0.40172292+3.97785014j, 
 1.26541377-2.2227806j ],
        [ 0.67562574+2.73428778j, 
-6.32196547-3.30680082j],
        [-0.75603444+1.17747523j, 
-0.02785408-1.70286801j],...