[SciPy-Dev] Many errors/failures running scipy.test() using latest numpy and scipy

[Benjamin Root]

On Tue, Jun 21, 2011 at 1:43 PM, Pauli Virtanen <pav at iki.fi> wrote:

> On Tue, 21 Jun 2011 12:23:54 -0500, Benjamin Root wrote:
> > Using the example given in the other mail thread:
> >
> > import numpy as np
> > import scipy.sparse.linalg
> > A = np.random.rand(6, 6) + 1j*np.random.rand(6, 6)
> > b = np.random.rand(6) + 1j*np.random.rand(6)
> > x, info = scipy.sparse.linalg.gmres(A, b)
> > y = np.linalg.solve(A, b)
> > print abs(x - y).max()
> >
> > produces:
> >
> > 0.427293469893
>
> Yep, that seems broken. I'll add some test cases for the iterative solvers
> so that we catch this if it occurs.
>
>        Pauli
>
>
I performed the tests on another machine of mine with completely different
specs (32-bit, gcc 4.5.2, Ubuntu 11.04), and the arpack errors were much
fewer.  I still had the other errors (like the netcdf one and the Anderson
one), but I attached the relevant ones here.

I hope this is helpful,
Ben Root
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mail.python.org/pipermail/scipy-dev/attachments/20110622/4e557018/attachment.html>
-------------- next part --------------
Running unit tests for scipy
NumPy version 2.0.0.dev-adab903
NumPy is installed in /home/ben/Programs/numpy/numpy
SciPy version 0.10.0.dev
SciPy is installed in /home/ben/Programs/scipy/scipy
Python version 2.7.1+ (r271:86832, Apr 11 2011, 18:05:24) [GCC 4.5.2]
nose version 0.11.4


======================================================================
ERROR: test_arpack.test_hermitian_modes(True, <gen-hermitian>, 'F', 2, 'SM', None, None, <function aslinearoperator at 0x971c5a4>)
----------------------------------------------------------------------
Traceback (most recent call last):
  File "/usr/lib/pymodules/python2.7/nose/case.py", line 186, in runTest
    self.test(*self.arg)
  File "/home/ben/Programs/scipy/scipy/sparse/linalg/eigen/arpack/tests/test_arpack.py", line 168, in eval_evec
    eval, evec = eigs_func(ac, k, bc, **kwargs)
  File "/home/ben/Programs/scipy/scipy/sparse/linalg/eigen/arpack/arpack.py", line 1425, in eigsh
    OPinv=OPinv)
  File "/home/ben/Programs/scipy/scipy/sparse/linalg/eigen/arpack/arpack.py", line 1251, in eigs
    params.iterate()
  File "/home/ben/Programs/scipy/scipy/sparse/linalg/eigen/arpack/arpack.py", line 727, in iterate
    self.workd[yslice] = self.OP(self.workd[xslice])
  File "/home/ben/Programs/scipy/scipy/sparse/linalg/eigen/arpack/arpack.py", line 641, in <lambda>
    self.OP = lambda x: Minv_matvec(matvec(x))
  File "/home/ben/Programs/scipy/scipy/sparse/linalg/interface.py", line 123, in matvec
    y = self._matvec(x)
  File "/home/ben/Programs/scipy/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_lu (test_decomp.TestLUSolve)
----------------------------------------------------------------------
Traceback (most recent call last):
  File "/home/ben/Programs/scipy/scipy/linalg/tests/test_decomp.py", line 732, in test_lu
    assert_array_equal(x1,x2)
  File "/home/ben/Programs/numpy/numpy/testing/utils.py", line 706, in assert_array_equal
    verbose=verbose, header='Arrays are not equal')
  File "/home/ben/Programs/numpy/numpy/testing/utils.py", line 635, in assert_array_compare
    raise AssertionError(msg)
AssertionError: 
Arrays are not equal

(mismatch 90.0%)
 x: array([-0.36278861,  0.75920432,  0.19844271,  0.65479681,  0.40510514,
        0.79604687, -1.05344321, -0.25597019, -0.34997286,  0.07080989])
 y: array([-0.36278861,  0.75920432,  0.19844271,  0.65479681,  0.40510514,
        0.79604687, -1.05344321, -0.25597019, -0.34997286,  0.07080989])