g5dZddlZddlZddlZddlmZmZmZddl m Z m Z m Z ddl mZmZddlmZddlmZddlmZdd lmZdd lmZdd lmZdd lmZGd deZdS)z$Dynamical Decoupling insertion pass.N)GateDelayReset)IGateUGateU3Gate) DAGOpNode DAGInNode) matrix_equal)OneQubitEulerDecomposer)TransformationPass)TranspilerError)InstructionDurations)Optimize1qGates)deprecate_funccleZdZdZedd dfd ZdZd Zd ed e d e fd Z xZ S)DynamicalDecouplinga Dynamical decoupling insertion pass. This pass works on a scheduled, physical circuit. It scans the circuit for idle periods of time (i.e. those containing delay instructions) and inserts a DD sequence of gates in those spots. These gates amount to the identity, so do not alter the logical action of the circuit, but have the effect of mitigating decoherence in those idle periods. As a special case, the pass allows a length-1 sequence (e.g. [XGate()]). In this case the DD insertion happens only when the gate inverse can be absorbed into a neighboring gate in the circuit (so we would still be replacing Delay with something that is equivalent to the identity). This can be used, for instance, as a Hahn echo. This pass ensures that the inserted sequence preserves the circuit exactly (including global phase). .. plot:: :alt: Output from the previous code. :include-source: import numpy as np from qiskit.circuit import QuantumCircuit from qiskit.circuit.library import XGate from qiskit.transpiler import PassManager, InstructionDurations from qiskit.transpiler.passes import ALAPSchedule, DynamicalDecoupling from qiskit.visualization import timeline_drawer # Because the legacy passes do not propagate the scheduling information correctly, it is # necessary to run a no-op "re-schedule" before the output circuits can be drawn. def draw(circuit): from qiskit import transpile scheduled = transpile( circuit, optimization_level=0, instruction_durations=InstructionDurations(), scheduling_method="alap", ) return timeline_drawer(scheduled) circ = QuantumCircuit(4) circ.h(0) circ.cx(0, 1) circ.cx(1, 2) circ.cx(2, 3) circ.measure_all() durations = InstructionDurations( [("h", 0, 50), ("cx", [0, 1], 700), ("reset", None, 10), ("cx", [1, 2], 200), ("cx", [2, 3], 300), ("x", None, 50), ("measure", None, 1000)] ) # balanced X-X sequence on all qubits dd_sequence = [XGate(), XGate()] pm = PassManager([ALAPSchedule(durations), DynamicalDecoupling(durations, dd_sequence)]) circ_dd = pm.run(circ) draw(circ_dd) # Uhrig sequence on qubit 0 n = 8 dd_sequence = [XGate()] * n def uhrig_pulse_location(k): return np.sin(np.pi * (k + 1) / (2 * n + 2)) ** 2 spacing = [] for k in range(n): spacing.append(uhrig_pulse_location(k) - sum(spacing)) spacing.append(1 - sum(spacing)) pm = PassManager( [ ALAPSchedule(durations), DynamicalDecoupling(durations, dd_sequence, qubits=[0], spacing=spacing), ] ) circ_dd = pm.run(circ) draw(circ_dd) zInstead, use :class:`~.PadDynamicalDecoupling`, which performs the same function but requires scheduling and alignment analysis passes to run prior to it.z1.1.0)additional_msgsinceNTc*t||_||_||_||_||_||_|C||_|D])}|j |j vrt|j d(dSdS)aDynamical decoupling initializer. Args: durations (InstructionDurations): Durations of instructions to be used in scheduling. dd_sequence (list[Gate]): sequence of gates to apply in idle spots. qubits (list[int]): physical qubits on which to apply DD. If None, all qubits will undergo DD (when possible). spacing (list[float]): a list of spacings between the DD gates. The available slack will be divided according to this. The list length must be one more than the length of dd_sequence, and the elements must sum to 1. If None, a balanced spacing will be used [d/2, d, d, ..., d, d, d/2]. skip_reset_qubits (bool): if True, does not insert DD on idle periods that immediately follow initialized/reset qubits (as qubits in the ground state are less susceptible to decoherence). target (Target): The :class:`~.Target` representing the target backend, if both ``durations`` and this are specified then this argument will take precedence and ``durations`` will be ignored. Nz. in dd_sequence is not supported in the target) super__init__ _durations _dd_sequence_qubits_spacing_skip_reset_qubits_target durationsnameoperation_namesr) selfr dd_sequencequbitsspacingskip_reset_qubitstargetgate __class__s n/granary/edev/quantum/lib/python3.11/site-packages/qiskit/transpiler/passes/scheduling/dynamical_decoupling.pyrzDynamicalDecoupling.__init__ns< #'  "3  $..00DO#  9F$:::)9TTT;    c h "t|jdks|jddtd|jtd||}t|j}d}|dkr|dzdkrtdtjd}|jD])}| | }*t|t d std tj |dd}|j4tt!||_nt|j|_|jrFt'|jdkst)d |jDrtd n d|z }|dz }|g|g|dz zz|gz|_t+|jD]>} |jD]4}||| s|j| n5?i} |jD]f} d} t1|jD]H\} }|}||j| <||| |_| |jz } I| | | <g|}|D]}t9|jt<s)||j|j |j!dF|j d}|"|j#} | |jvr)||j|j |j!dtI|%|}tI|&|}|j'rYt9|tPst9|jtRr*||j|j |j!dC| | } |jj| z ""dkr*||j|j |j!d|dkr_|jd* }tW,|\}}}}t9|tZrZt9|jt\t^fr9|jj0\}}}tcj2|||||||j_0||z }nt9|tZrZt9|jt\t^fr9|jj0\}}}tcj2|||||||j_0||z }n*||j|j |j!d"fd|jD}"t'|z }tgt|dz dz } || xx|z cc<tij5||jD]L\}!}|!dkr&|t=|!|gd||||gdM|j6|z|_6|S)aJRun the DynamicalDecoupling pass on dag. Args: dag (DAGCircuit): a scheduled DAG. Returns: DAGCircuit: equivalent circuit with delays interrupted by DD, where possible. Raises: TranspilerError: if the circuit is not mapped on physical qubits. qNz"DD runs on physical circuits only.z#DD runs after circuit is scheduled.rz8DD sequence must contain an even number of gates (or 1).T) ignore_phasez4The DD sequence does not make an identity operation.c3"K|] }|dkV dS)rN).0as r* z*DynamicalDecoupling.run..s&-K-Ka!e-K-K-K-K-K-Kr+zRThe spacings must be given in terms of fractions of the slack period and sum to 1.F)checkc4g|]}t|zSr2)int)r3r4slacks r* z+DynamicalDecoupling.run..s#:::qC NN:::r+)7lenqregsgetrduration_update_inst_durationsrnpeyedot to_matrixr ranglersetrange num_qubitsrsumanylist$_DynamicalDecoupling__gate_supporteddiscard enumerate to_mutablecopy_empty_liketopological_op_nodes isinstanceoprapply_operation_backqargscargsfind_bitindexnext predecessors successorsrr rinverser angles_and_phaser rrparamsr compose_u3r8 itertools zip_longest global_phase)#r"dagr num_pulsessequence_gphasenoopr(midendqargindex_sequence_duration_mapphysical_qubitdd_sequence_durationrWnew_dagnd dag_qubitpredsuccu_invthetaphilamphasetheta_rphi_rlam_rtheta_lphi_llam_ltaus unused_slack middle_indextaur9s# @r*runzDynamicalDecoupling.runs sy>>Q  #)--T":":"B!"FGG G < !"GHH H//44 *++  ??A~""%&`aaa6!99D) 2 2xx 0 011egg&7&7&9&9MMM ^%&\]]] htAwqz22O < uS^^%5%56677DLLt|,,DL = E4=!!Q&&#-K-KT]-K-K-K*K*K&%8' j.C'C ESEZ!^$<<uDDM&&  D)  ,,T488L((...E')#"l O ON#$ ():;; 6 6 t((+/!%( ) dN C C $ 5$$:N ' 7 7%%''**,,: J: JBbeU++ ,,RUBHbhe,TTT I \\)44:NT\11,,RUBHbhe,TTT((,,--Dr**++D& dI..*TWe2L2L00"(RW0XXX#>~#N EN%99Ezz,,RUBHbhe,TTTQ)!,4466@@BB)@)B)B)S)STY)Z)Z&sCdI..:dgv3W3W,0GN)GUE%4%?uc3&&DGN$u,OOi00 Z%QW5Y5Y ,0GN)GUE%4%?sC%&&DGN$u,OO00"(RW0XXX;:::DM:::D 3t99,LD A 233L    , .   &249JKK Q Q T7700si[PU0VVV#00 {%0PPP#*#7/#IG r+c t}|jrg}tj5tjdt |jD]@\}}|D]&\\}}}|||||jf'A dddn #1swxYwY| ||j |j /| |j t|j dd|S)aUUpdate instruction durations with circuit information. If the dag contains gate calibrations and no instruction durations were provided through the target or as a standalone input, the circuit calibration durations will be used. The priority order for instruction durations is: target > standalone > circuit. ignore)actioncategoryNdt) r_calibrations_propwarningscatch_warnings simplefilterDeprecationWarningitemsappendr>updaterrgetattr) r"rbcirc_durations cal_durationsr( gate_calsr$ parametersschedules r*r?z*DynamicalDecoupling._update_inst_durationssv .//  ! DM(** \ \%X@RSSSS(+'='C'C'E'E\\OD):C//:K:K\\6,h%,,dFJHY-Z[[[[\\  \ \ \ \ \ \ \ \ \ \ \ \ \ \ \  ! !-1B C C C ? &  ! !$/74?DRV3W3W X X XsA6B..B25B2r(rhreturnc\|j"|j|j|frdSdS)z/A gate is supported on the qubit (qarg) or not.N)rTTF)rinstruction_supportedr )r"r(rhs r*__gate_supportedz$DynamicalDecoupling.__gate_supported5s4 < 4<#E#EdiX\W^#E#_#_ 4ur+)NNTN) __name__ __module__ __qualname____doc__rrrr?rr8boolrK __classcell__)r)s@r*rrsLL\^ a ae$$$$$$L@@@D0Tr+r)rr_rnumpyr@qiskit.circuitrrr%qiskit.circuit.library.standard_gatesrrrqiskit.dagcircuitr r (qiskit.quantum_info.operators.predicatesr qiskit.synthesis.one_qubitr qiskit.transpiler.basepassesr qiskit.transpiler.exceptionsr'qiskit.transpiler.instruction_durationsr%qiskit.transpiler.passes.optimizationrqiskit.utils.deprecationrrr2r+r*rs?+*----------FFFFFFFFFF22222222AAAAAA>>>>>>;;;;;;888888HHHHHHAAAAAA333333ZZZZZ,ZZZZZr+