'solve_triangular on a simple 2x2 complex matrix'

def test_simple_complex(self):

    A = array([[(1 + 1j), 0], [1j, 2]])
    b = identity(2)
    sol = solve_triangular(A, b, lower=True, trans=1)
    assert_array_almost_equal(sol, [[(0.5 - 0.5j), ((-0.25) - 0.25j)], [0, 0.5]])
