f8a349522c5d07b09d4825fa835c23ef456f991d,scipy/sandbox/multigrid/tests/test_relaxation.py,TestRelaxation,check_jacobi,#TestRelaxation#,40

Before Change

        assert_almost_equal(x,array([0.5,1.0,0.5]))

        N = 1
        A = spdiags([2*ones(N),-ones(N),-ones(N)],[0,-1,1],N,N).T
        x = arange(N).astype(numpy.float64)
        b = array([10])
        jacobi(A,x,b)
        assert_almost_equal(x,array([5]))

After Change

        assert_almost_equal(x,array([0.5,1.0,0.5]))

        N = 1
        A = spdiags([2*ones(N),-ones(N),-ones(N)],[0,-1,1],N,N,format="csr")
        x = arange(N).astype(numpy.float64)
        b = array([10])
        jacobi(A,x,b)
        assert_almost_equal(x,array([5]))

Instances