def test_decompose_gray_random(self):
for matrix_size in [2, 4, 8, 16]:
for i in range(4):
A = np.array(unitary_group.rvs(matrix_size))
self.check_correct_product(A, qd.two_level_decompose(A))
def check_two_level_decompose(self, matrix):
matrix = np.array(matrix)
self.check_correct_product(matrix, qd.two_level_decompose(matrix))