

def test_random(self):
    n = 20
    a = random([n, n])
    for i in range(n):
        a[(i, i)] = (20 * (0.1 + a[(i, i)]))
    for i in range(4):
        b = random([n, 3])
        x = solve(a, b)
        assert_array_almost_equal(dot(a, x), b)

dataset/ETHPy150Open scipy/scipy/scipy/linalg/tests/test_basic.py TestSolve.test_random/original