g?dZddlZddlZddlmZmZmZmZddlZ ddl m Z ddl m Z mZmZmZmZmZmZddlmZddlmZddlmZmZmZdd lmZdd lmZd d l m!Z!d d l"m#Z#d dl$m%Z%ej&e'Z(de#dede#fdZ)de!de!fdZ*de#de#fdZ+dddddZ,dddddZ-e dgfgedgfgedgfgedgfggZ.e dgfgedgfgedgfgedgfggZ/de.DZ0e dgfgZ1e dgfedgfgZ2e dgfge1e2gZ3e dgfgZ4e dgfedgfgZ5e dgfge4e5gZ6de3DZ7e fefefefefeefeefeefeefeefeefeeefeeefeeefeeefeeefeeefeeefeeefeeefefefeefeefgZ8e.dde0ddde3dde7ddddde9ddDidZ:d eeee!fd!eee!fdee;fd"ZdS)%a Noise transformation module The goal of this module is to transform one 1-qubit noise channel (given by the QuantumError class) into another, built from specified "building blocks" (given as Kraus matrices) such that the new channel is as close as possible to the original one in the Hilber-Schmidt metric. For a typical use case, consider a simulator for circuits built from the Clifford group. Computations on such circuits can be simulated at polynomial time and space, but not all noise channels can be used in such a simulation. To enable noisy Clifford simulation one can transform the given noise channel into the closest one, Hilbert-Schmidt wise, that can be used in a Clifford simulator. N)SequenceListUnionCallable)Reset)IGateXGateYGateZGateHGateSGateSdgGate) transpile)MissingOptionalLibraryError)KrausSuperOpChi)QuantumChannel)TranspilerError) QuantumError) NoiseModel) NoiseError noise_modelfuncreturnc6tj|}|jD]\}}|||j|<|jD]8\}}|D]\}}|||j||<9|S)aReturn a new noise model by applyign a function to all quantum errors. This returns a new noise model containing the resulting errors from applying the supplied function to all QuantumErrors in the noise model. This function should have singature `func(error: QuantumError) -> QuantumError` where the returned quantum error is defined on the same number of qubits as the original error. Args: noise_model: the noise model to be transformed. func: function for transforming QuantumErrors. Returns: The transpiled noise model. Raises: NoiseError: if the transformation failed. )copydeepcopy_default_quantum_errorsitems_local_quantum_errors)rr new_noisekeyerror inst_name noise_dicqubitss [/granary/edev/quantum/lib/python3.11/site-packages/qiskit_aer/utils/noise_transformation.pytransform_noise_modelr*:s( k**I7==??== U15e )#.. ) ? E E G GMM 9&__.. M MMFEAEeI +I 6v > > M r%c t|jfi|}n(#t$r}td|d||d}~wwxYwt t ||jS)aReturn a new quantum error containin transpiled circuits. This returns a new QuantumError containing the circuits resulting from transpiling all error circuits using :func:`qiskit.transpile` with the passed keyword agruments. Args: error: the quantum error to be transformed. transpile_kwargs: kwargs for passing to qiskit transpile function. Returns: The transformed quantum error. Raises: NoiseError: if the transformation failed. z Failed to transpile circuits in z with kwargs N)rcircuitsrrrzip probabilities)r%transpile_kwargstranspiled_circserrs r)transpile_quantum_errorr3Zs"$U^HH7GHH  Uu U UCS U U    ,e.ABB C CCs :5:c .fd}t||S)aReturn a new noise model with transpiled QuantumErrors. This returns a new noise model containing the resulting errors from transpiling all QuantumErrors in the noise model using :func:`transpile_quantum_error` function with the passed keyword agruments. Args: noise_model: the noise model to be transformed. transpile_kwargs: kwargs for passing to qiskit transpile function. Returns: The transpiled noise model. Raises: NoiseError: if the transformation failed. ct|fiSN)r3)r%r0s r)rz#transpile_noise_model..funcs&uAA0@AAAr+r*)rr0rs ` r)transpile_noise_modelr8ts3&BBBBB !d 3 33r+operator_string operator_dict operator_listcJt|ttfstd|jj|jdkrtd|jd|t}| }||vrt|d| t||j}n.#t$r!}td|d|jd |d}~wwxYw|t| \}}|t|tstd | d |D}n%#t$r}td ||d}~wwxYwt||d d}tjd t#|z d} | dks| d krtdd | z t%dgf| fg} t||D]\} } | | | ft| Std)a Return a ``QuantumError`` object that approximates an error as a mixture of specified operators (channels). The approximation is done by minimizing the Hilbert-Schmidt distance between the process matrix of the target error channel (:math:`T`) and the process matrix of the output channel (:math:`S = \sum_i{p_i S_i}`), i.e. :math:`Tr[(T-S)^\dagger (T-S)]`, where :math:`[p_1, p_2, ..., p_n]` denote probabilities and :math:`[S_1, S_2, ..., S_n]` denote basis operators (channels). See `arXiv:1207.0046 `_ for the details. Args: error (QuantumError or QuantumChannel): the error to be approximated. The number of qubits must be 1 or 2. operator_string (string): a name for a pre-made set of building blocks for the output channel (Default: None). Possible values are ``'pauli'``, ``'reset'``, ``'clifford'``. operator_dict (dict): a dictionary whose values are the building blocks for the output channel (Default: None). E.g. {"x": XGate(), "y": YGate()}, keys "x" and "y" are not used in transformation. operator_list (list): list of building block operators for the output channel (Default: None). E.g. [XGate(), YGate()] Returns: QuantumError: the approximate quantum error. Raises: NoiseError: if any invalid argument is specified or approximation failed. MissingOptionalLibraryError: if cvxpy is not installed. Note: The operator input precedence is: ``list`` < ``dict`` < ``string``. If a ``string`` is given, ``dict`` is overwritten; if a ``dict`` is given, ``list`` is overwritten. The ``string`` supports only 1- or 2-qubit errors and its possible values are ``'pauli'``, ``'reset'``, ``'clifford'``. The ``'clifford'`` does not support 2-qubit errors. zInvalid input error type: rz2Only 1-qubit and 2-qubit noises can be converted, z-qubit noise found in modelNz3 is not a valid operator_string. It must be one of zPreset 'z<' operators do not support the approximation of errors with z qubitsz!operator_list is not a sequence: c`g|]+}t|tr|nt|dfg,S)) isinstancerr.0ops r) z-approximate_quantum_error..sI!^44Q,Qy:Q:Qr+z%Invalid type found in operator list: r? rzChannel probabilities sum to zHQuantum error approximation failed - no approximating operators detected)r@rrr __class____name__ num_qubits_PRESET_OPERATOR_TABLEkeyslowerKeyErrorr.r!r_transform_by_operator_listnproundsumrappend) r%r:r;r<valid_operator_stringsr2_ channel_listr/ identity_prob noise_opsoperator probabilitys r)approximate_quantum_errorrYsZ elN; < <RPeo6NPPQQQ ! = = = =   "!7!!>)//11 "8 8 8">>%;>>  2?CEDTUMM   J?JJ050@JJJ     3 3 5 56= -22 RPPPQQ Q _'LL _ _ _T]TTUU[^ ^ _4L%HHL S%7%7!7;; 1    1 1PQ=NPPQQ Qwwnm45 %( %F%F 6 6 !Hk   h 4 5 5 5 5I&&& _ ` ``s0$B== C(C##C(6 E E% E  E%c6fd}t||S)a Replace all noises in a noise model with ones approximated by a mixture of operators (channels). Args: model (NoiseModel): the noise model to be approximated. All noises in the model must be 1- or 2-qubit noises. operator_string (string): a name for a pre-made set of building blocks for the output channel (Default: None). Possible values are ``'pauli'``, ``'reset'``, ``'clifford'``. operator_dict (dict): a dictionary whose values are the building blocks for the output channel (Default: None). E.g. {"x": XGate(), "y": YGate()}, keys "x" and "y" are not used in transformation. operator_list (list): list of building block operators for the output channel (Default: None). E.g. [XGate(), YGate()] Returns: NoiseModel: the approximate noise model. Raises: NoiseError: if any invalid argument is specified or approximation failed. MissingOptionalLibraryError: if cvxpy is not installed. Note: The operator input precedence is: ``list`` < ``dict`` < ``string``. If a ``string`` is given, ``dict`` is overwritten; if a ``dict`` is given, ``list`` is overwritten. The ``string`` supports only 1- or 2-qubit errors and its possible values are ``'pauli'``, ``'reset'``, ``'clifford'``. The ``'clifford'`` does not support 2-qubit errors. c*t|S)Nr9)rY)noiser;r<r:s r) approximatez,approximate_noise_model..approximate s&( +''     r+r7)modelr:r;r<r]s ``` r)approximate_noise_modelr_s;D        ! 4 44r+r?c,g|]}tD]}||zS) _PAULIS_Q1rBop_q0op_q1s r)rDrDs)NNN%:NN% NNNNr+c,g|]}tD]}||zSra) _RESET_Q1rcs r)rDrD#s)KKKKKuu}KKKKr+)r?rc<g|]}dt|DS)cg|]}|dgfS)rra)rBgates r)rDz.Ls 8 8 8TdQC[ 8 8 8r+)_CLIFFORD_GATES)rBjs r)rDrDLs, O O OQ 8 8_Q%7 8 8 8 O O Or+)pauliresetclifford basis_opstargetcrd|djz}ttj|g|z}d|D}t |jt |}tj||f}t|D]p\}}t|D][\}} ||krt|| |||<%||kr||||||<Ct||||<\qdtj fd|Dz} dtj fd|D} |} ddl } n@#t$r3}td t!d d d d |d}~wwxYw| |}| | | ||| j|zz| |dk|dk| j|z| kg}||j}|S)aD Transform (or approximate) the target quantum channel into a mixture of basis operators (channels) and return the mixing probabilities. The approximate channel is found by minimizing the Hilbert-Schmidt distance between the process matrix of the target channel (:math:`T`) and the process matrix of the output channel (:math:`S = \sum_i{p_i S_i}`), i.e. :math:`Tr[(T-S)^\dagger (T-S)]`, where :math:`[p_1, p_2, ..., p_n]` denote probabilities and :math:`[S_1, S_2, ..., S_n]` denote basis operators (channels). Such an optimization can represented as a quadratic program: Minimize :math:`p^T A p + b^T p`, subject to :math:`p \ge 0`, `\sum_i{p_i} = 1`, `\sum_i{p_i} = 1`. The last constraint is for making the approximation conservative by forcing the output channel have more error than the target channel in the sense that a "fidelity" against identity channel should be worse. See `arXiv:1207.0046 `_ for the details. Args: target: Quantum channel to be transformed. basis_ops: a list of channel objects constructing the output channel. Returns: list: A list of amplitudes (probabilities) of basis that define the output channel. Raises: MissingOptionalLibraryError: if cvxpy is not installed. rrc6g|]}t|jSra)rdatarAs r)rDz/_transform_by_operator_list..zs 777rc"ggl777r+c0g|]}t|Sra)_hs_inner_prod_real)rBS_iTs r)rDz/_transform_by_operator_list..s$MMMs*1c22MMMr+c ttjtjt |Sr6)floatrNabstracer)channels r)fidelityz-_transform_by_operator_list..fidelitys,RVBHWW%5%56677888r+c&g|] }|Srara)rBrCrs r)rDz/_transform_by_operator_list..s!===RHHRLL===r+Nz7cvxpy module needs to be installed to use this feature.cvxpyz%Transformation/Approximation of noisezpip install cvxpyzTCVXPY is required to solve an optimization problem of approximating a noise channel.)libnamename pip_installmsgr?) constraints)rHrrNeyerrulenzeros enumeraterx_hs_normarrayr ImportErrorloggerr%rVariableProblemMinimize quad_formrzrPsolvevalue)rqrrNbasis_ops_matsnAiryrlS_jb source_fids target_fidrr2xprobr/rzrs @@r)rMrMQs}H Yq\ $$Arvayy!!"Y.I87Y777N F A NA !QAN++((3// ( (FAs1uu-c377!QQA$q'!Q"3--!Q ( RXMMMMnMMM N NNA999(====9===>>K&!!J   NOOO)8+.       qA == uq!,,qsQw677YYq\\Q&Q 0AZ0OP   D  JJLLLGM sE F.F  Fchtjtj|j|zjSr6rNr~conjrzreal)rs r)rrs$ 8BGAJJL1$ % % **r+chtjtj|j|zjSr6r)rBs r)rxrxs$ 8BGACLL1$ % % **r+)?__doc__rloggingtypingrrrrnumpyrNqiskit.circuitr%qiskit.circuit.library.standard_gatesrr r r r r rqiskit.compilerrqiskit.exceptionsr%qiskit.quantum_info.operators.channelrrr5qiskit.quantum_info.operators.channel.quantum_channelrqiskit.transpiler.exceptionsr noise.errorsrnoise.noise_modelrnoise.noiseerrorr getLoggerrGrr*r3r8rYr_ _PAULIS_Q0rb _PAULIS_Q0Q1_RESET_Q0_TO_0_RESET_Q0_TO_1 _RESET_Q0_RESET_Q1_TO_0_RESET_Q1_TO_1rg _RESET_Q0Q1rkrangerIr|rMrrxrar+r)rsx" 222222222222 &%%%%%999999EEEEEEEEEEPPPPPP888888''''''******))))))  8 $ $zj@D<D DDDD44z4*44444#$d[a[a[a[a[a|7;$^b*5*5*5*5*5\!~%%''A3 0EEGGaS>2BeeggPQs^DT U !~%%''A3 0EEGGaS>2BeeggPQs^DT U NN:NNN 577QC.!577QC.5577QC.1uwwn ~~ > 577QC.!577QC.5577QC.1uwwn ~~ > KK)KKK  UWWJ UWWJ WYYL UWWJ UWWJ UWWeegg UWWeegg UWWeegg UWWeegg UWWggii WYY UWWeegguuww WYY! UWWeegguuww UWWeeggwwyy! WYY# UWWeeggwwyy! UWWeegguuww WYY! UWWeegguuww UWWJ UWWJ UWWeegg WYY5> abb>   QRR= qrr? O O%%2,, O O O  Tnl:;<T .,. /T %[TTTTn+++ +++++r+