def METHOD_NAME(run_number):
    backend = NumpyBackend()
    unitary = random_unitary(2)
    c = Circuit(2)
    c.add(two_qubit_decomposition(0, 1, unitary))
    final_matrix = c.unitary(backend)
    np.testing.assert_allclose(final_matrix, unitary, atol=ATOL)