

def test_complex(self):
    c = np.array([(1 + 2j), (- 3), 4j, 5])
    b = (np.arange(8).reshape(4, 2) + 0.5j)
    x = solve_circulant(c, b)
    y = solve(circulant(c), b)
    assert_allclose(x, y)

dataset/ETHPy150Open scipy/scipy/scipy/linalg/tests/test_basic.py TestSolveCirculant.test_complex/original