

def test_random_sym(self):
    n = 20
    a = random([n, n])
    for i in range(n):
        a[(i, i)] = abs((20 * (0.1 + a[(i, i)])))
        for j in range(i):
            a[(i, j)] = a[(j, i)]
    for i in range(4):
        b = random([n])
        x = solve(a, b, sym_pos=1)
        assert_array_almost_equal(dot(a, x), b)

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