def METHOD_NAME(
    q1: int,
    q2: int,
    args: sim.CliffordTableauSimulationState,
    operations: List[ops.Operation],
    qubits: List['cirq.Qid'],
):
    protocols.act_on(ops.CNOT, args, qubits=[qubits[q1], qubits[q2]], allow_decompose=False)
    operations.append(ops.CNOT(qubits[q1], qubits[q2]))