def test_solve_discrete_lyapunov_complex():
 A = np.array([[(0.5 + 0.3j), (0.1 + 0.1j)], [1, 0]])
 B = np.eye(2)
 X = qme.solve_discrete_lyapunov(A, B)
 assert_allclose((np.dot(np.dot(A, X), A.conj().transpose()) - X), (- B), atol=1e-15)

def test_solve_discrete_lyapunov_complex ( ) : A = np . array ( [ [ ( 0.5 + 0.3j ) , ( 0.1 + 0.1j ) ] , [ 1 , 0 ] ] ) B = np . eye ( 2 ) X = qme . solve_discrete_lyapunov ( A , B ) assert_allclose ( ( np . dot ( np . dot ( A , X ) , A . conj ( ) . transpose ( ) ) - X ) , ( - B ) , atol = 1e-15 )
complex test .