[SciPy-Dev] Many errors/failures running scipy.test() using latest numpy and scipy
Pauli Virtanen
pav at iki.fi
Tue Jun 21 14:43:33 EDT 2011
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
More information about the SciPy-Dev
mailing list