gsdZddlmZddlmZddlmZddlZddl m Z ddl m Z ddl mZdd lmZmZmZdd lmZdd lmZdd lmZmZerdd lmZGddeZdS)z Stabilizer state class. ) annotations) Collection) TYPE_CHECKINGN QiskitError)OpShape)Operator)CliffordPauli PauliList) _append_x) QuantumState)QuantumCircuit Instruction)circuitceZdZdZ dFdGfd Ze dHdIdZdZdZe dZ dJdZ dZ dZ dKdZdKdZdLdZdZdMdZdMdZ dNdOd"ZdNdPd'ZdQd(ZdJdRd,Z dJdSd0Z dJdTd1ZdNdUd2ZdNdVd4ZdNdWd7Zd8Zed9Zed:Z ed;Z!ed<Z" dNdXdCZ# dYdZdEZ$xZ%S)[StabilizerStateaStabilizerState class. Stabilizer simulator using the convention from reference [1]. Based on the internal class :class:`~qiskit.quantum_info.Clifford`. .. code-block:: from qiskit import QuantumCircuit from qiskit.quantum_info import StabilizerState, Pauli # Bell state generation circuit qc = QuantumCircuit(2) qc.h(0) qc.cx(0, 1) stab = StabilizerState(qc) # Print the StabilizerState print(stab) # Calculate the StabilizerState measurement probabilities dictionary print (stab.probabilities_dict()) # Calculate expectation value of the StabilizerState print (stab.expectation_value(Pauli('ZZ'))) .. code-block:: text StabilizerState(StabilizerTable: ['+XX', '+ZZ']) {'00': 0.5, '11': 0.5} 1 Given a list of stabilizers, :meth:`qiskit.quantum_info.StabilizerState.from_stabilizer_list` returns a state stabilized by the list .. code-block:: python from qiskit.quantum_info import StabilizerState stabilizer_list = ["ZXX", "-XYX", "+ZYY"] stab = StabilizerState.from_stabilizer_list(stabilizer_list) References: 1. S. Aaronson, D. Gottesman, *Improved Simulation of Stabilizer Circuits*, Phys. Rev. A 70, 052328 (2004). `arXiv:quant-ph/0406196 `_ TdataUStabilizerState | Clifford | Pauli | QuantumCircuit | circuit.instruction.Instructionvalidateboolcnt|tr |j|_nQt|tr't ||_nt |||_t tj |jj ddS)zInitialize a StabilizerState object. Args: data: Data from which the stabilizer state can be constructed. validate: validate that the stabilizer state data is a valid Clifford. r) num_qubits_r num_qubits_l)op_shapeN) isinstancer_datar r to_instructionsuper__init__rauto num_qubits)selfrr __class__s `/granary/edev/quantum/lib/python3.11/site-packages/qiskit/quantum_info/states/stabilizerstate.pyr zStabilizerState.__init__Ts dO , , 2DJJ e $ $ 2!$"5"5"7"788DJJ"$11DJ ',DJstz))Q #F++ 5 5JC!VC[EGEN!VC[EGEN UWU^egen4 4EEz""  AH (4=3G3J)JKK KH (4=3G3J)JKK KH!|q  qq! 7..""z"" 8 8AH (4=3I!3L)LMM MH (4=3I!3L)LMM MH!|q   Q1!44 4E R%dm&:1&= @TUV@W&WXX XE Q)'DM4H4K*KLLL LE 4Q 77GG 19>>< r&c*t|ts" t|}n#t$rYdSwxYw|j}|j|krdSt j|jj|jjd|jj z}t j|jj|jjd|jj ztj fd|DsdSt|D]&}| |}|dkrdS'dS)aReturn True if the two generating sets generate the same stabilizer group. Args: other (StabilizerState): another StabilizerState. Returns: bool: True if other has a generating set that generates the same StabilizerState. Frnc:g|]}|S)commutes).0r} pauli_others r% z)StabilizerState.equiv.._s%KKKuu~~k22KKKr&rmT)rrrr"r from_symplecticrrxrwr{rur4rpr)r#r6r" pauli_origiexp_valrs @r%equivzStabilizerState.equivCs?%11  '..   uu _  z ) )5. J tz0!dj6K2K   / K  2A 8N4N  vKKKK KKKLL 5 z""  A,,[^<** H HJqM%(HqG||!R\H\B]B]11(AwGGG4Q3RRVWW 199 QJE A::YZZ Z r&c d|zd|zz}t|jD]N}|t|j||j||j||j|z }O|dz}|dk|dkzrt dt|dk}|xj|jzc_|xj|jzc_||fS)z)Aaronson-Gottesman rowsum helper functionrnrorz*Invalid rowsum in measurement calculation.)rpr"rrrsrtrr) accum_pauli accum_phase row_pauli row_phasenewrrs r%_rowsumzStabilizerState._rowsum_s9}q;.9/00  E O33 E"IK$6 e8Lkm\aNb DD   AI$!) $ LJKK K$!)nn  $  $ K''r&cF|j|}|j|}|j}|j}t||||f}t||||f}t||||\}}||j|<|j||<|j||<dS)zUpdating StabilizerState Clifford in the non-deterministic rowsum calculation. row and accum are rows in the StabilizerState Clifford.N)rqrtrsr rr) r@accumrowrrrtrsrrs r%rz(StabilizerState._rowsum_nondeterministicrs N3' nU+ J J1S61S6*++ QuXqx011 #2#:#: i$ $  [!,u=%=%r&c|j|}|j}|}t|j||j|f}t||||\}}|}||_|S)zUpdating an auxiliary Pauli aux_pauli in the deterministic rowsum calculation. The StabilizerState itself is not updated.)rqr rtrsrr)r@rrrrrrs r%rz%StabilizerState._rowsum_deterministicsq N3' o  8:c?HJsO<== #2#:#: i$ $  [ % r&r|rpr list[str] outcome_probrcjd}|}tt|D]}||dkr|t||z dz } t|jjdd| fs_|| d} |t||| krt| ||<t||||<d}|}|dkrd |} ||| <dS|tddnt||gD]} |} t| | |<|} | |t||z dz | | || d|z||dS) aRecursive helper function for calculating the probabilities Args: qubits (range): range of qubits outcome (list[str]): outcome being built outcome_prob (float): probability of the outcome probs (dict[str, float]): holds the outcomes and probability results outcome_bitstring (str): target outcome to measure which reduces measurements, None if not targeting a specific target XrmNrrrng?) rYrprrr@rwrrrjoin_get_probabilities)r#r|rrrrqubit_for_branchingr[rrsingle_qubit_outcome str_outcome new_outcomestab_cpys r%rz"StabilizerState._get_probabilitiess$$&#yy{{s6{{## , ,AqzS  F a! 353<.qqq%x899,585L5LUTU5V5V((0-a0115III &))=%>%> &)):1)=%>%> '( *+' " $ $'''**K!-E+  F") !QKKK'(;<==>   ",,..K/23G/H/HK+ ,xxzzH  ( (F &99A=?AU     ' ' cL&85BS      r& None | strc|t|jj}n|}dgt|z}d}i}|||||||-|D]\}} t | |||<|S)aIHelper Function for calculating the subsystem measurement probability dictionary. When the targeted outcome_bitstring value is set, then only the single outcome_bitstring probability will be calculated. Args: outcome_bitstring (None or str): targeted outcome bitstring to perform a measurement calculation for, this will significantly reduce the number of calculation performed (Default: None) qargs (None or list): subsystems to return probabilities for, if None return for all subsystems (Default: None). decimals (None or int): the number of decimal places to round values. If None no rounding is done (Default: None). Returns: dict: The measurement probabilities in dict (key) form. NrrP)rpr@r"rrrround) r#rrcrr|rrrrrs r%rz'StabilizerState._get_probabilities_dicts, =4=344FFF%#f++% "$ uFWXXX  #kkmm 4 4 U"5(33c  r&)T)rrrr)FF)r'r(r)rr*rr+r)NN)r+rM)r+r )r6rr+rr2)r6rbrcrdr+r)rhr rcrir+rj)r6rr+r)rcrirrr+r)rrrcrirrr+r)rcrirrr+r)rcrdr+r)rcrdr+r)rrrcrir+r) r|rprrrrMrrrr)NNN)rrrcrirrr+r)&__name__ __module__ __qualname____doc__r classmethodr0r7r>propertyr@rErJrLrRrUrWrZr]rargrrrrrrrrr staticmethodrrrrrr __classcell__)r$s@r%rr$s,--dddddddd2!&', [:;;;DDDX''''MMMZZZ        5555 $&TX0?????B''''R:"# " " " " " JAEddddd0#####J<:...`  \ ((\($!!\!*\8"& @@@@@H)-!# %%%%%%%%%r&r)r __future__rcollections.abcrtypingrnumpyruqiskit.exceptionsr&qiskit.quantum_info.operators.op_shaper&qiskit.quantum_info.operators.operatorr (qiskit.quantum_info.operators.symplecticr r r :qiskit.quantum_info.operators.symplectic.clifford_circuitsr (qiskit.quantum_info.states.quantum_staterqiskit.circuitrrqiskitrrrr&r%rs;#"""""&&&&&& ))))))::::::;;;;;;OOOOOOOOOOPPPPPPAAAAAA66666666b b b b b lb b b b b r&