def test_solve_discrete_lyapunov_B():
 A = (np.ones((2, 2)) * 0.5)
 B = np.array([[0.5, (-0.5)], [(-0.5), 0.5]])
 X = qme.solve_discrete_lyapunov(A, B)
 assert_allclose(B, X)

def test_solve_discrete_lyapunov_B ( ) : A = ( np . ones ( ( 2 , 2 ) ) * 0.5 ) B = np . array ( [ [ 0.5 , ( - 0.5 ) ] , [ ( - 0.5 ) , 0.5 ] ] ) X = qme . solve_discrete_lyapunov ( A , B ) assert_allclose ( B , X )
simple test where x is same as b .