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

Pauli Virtanen pav@iki...
Tue Jun 21 13:43:33 CDT 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