gNKdZddlmZddlZddlZddlmZddlmZddl m Z ddl Z ddl mZmZmZddlmZdd lmZmZdd lmZmZdd lmZdd lmZmZdd lmZddl m!Z!m"Z"m#Z#ddl$m%Z%ddl&m'Z'm(Z(m)Z)m*Z*ddl+m,Z,ddl-m.Z.ddl/m0Z0e GddZ1e GddZ2Gdde%Z3d dZ4dS)!zdZ d?d$Z d@d,ZdAd1ZdBd6ZdS)CBackendEstimatorV2a Evaluates expectation values for provided quantum circuit and observable combinations. The :class:`~.BackendEstimatorV2` class is a generic implementation of the :class:`~.BaseEstimatorV2` interface that is used to wrap a :class:`~.BackendV2` (or :class:`~.BackendV1`) object in the :class:`~.BaseEstimatorV2` API. It facilitates using backends that do not provide a native :class:`~.BaseEstimatorV2` implementation in places that work with :class:`~.BaseEstimatorV2`. However, if you're using a provider that has a native implementation of :class:`~.BaseEstimatorV2`, it is a better choice to leverage that native implementation as it will likely include additional optimizations and be a more efficient implementation. The generic nature of this class precludes doing any provider- or backend-specific optimizations. This class does not perform any measurement or gate mitigation, and, presently, is only compatible with Pauli-based observables. More formally, given an observable of the type :math:`O=\sum_{i=1}^Na_iP_i`, where :math:`a_i` is a complex number and :math:`P_i` is a Pauli operator, the estimator calculates the expectation :math:`\mathbb{E}(P_i)` of each :math:`P_i` and finally calculates the expectation value of :math:`O` as :math:`\mathbb{E}(O)=\sum_{i=1}^Na_i\mathbb{E}(P_i)`. The reported ``std`` is calculated as .. math:: \frac{\sum_{i=1}^{n}|a_i|\sqrt{\textrm{Var}\big(P_i\big)}}{\sqrt{N}}\:, where :math:`\textrm{Var}(P_i)` is the variance of :math:`P_i`, :math:`N=O(\epsilon^{-2})` is the number of shots, and :math:`\epsilon` is the target precision [1]. Each tuple of ``(circuit, observables, parameter values, precision)``, called an estimator primitive unified bloc (PUB), produces its own array-based result. The :meth:`~.BackendEstimatorV2.run` method can be given a sequence of pubs to run in one call. The options for :class:`~.BackendEstimatorV2` consist of the following items. * ``default_precision``: The default precision to use if none are specified in :meth:`~run`. Default: 0.015625 (1 / sqrt(4096)). * ``abelian_grouping``: Whether the observables should be grouped into sets of qubit-wise commuting observables. Default: True. * ``seed_simulator``: The seed to use in the simulator. If None, a random seed will be used. Default: None. **Reference:** [1] O. Crawford, B. van Straaten, D. Wang, T. Parks, E. Campbell, St. Brierley, Efficient quantum measurement of Pauli operators in the presence of finite sampling error. `Quantum 5, 385 `_ N)optionsbackendBackendV1 | BackendV2r6 dict | Nonect|ts tjd|dtd||_|r t di|n t |_tj |j }t|trt||j }nt|}t|g|_dS) a .. deprecated:: 1.4 The method ``BackendEstimatorV2.__init__`` will stop supporting inputs of type :class:`.BackendV1` in the `backend` parameter in a future release no earlier than 2.0. :class:`.BackendV1` is deprecated and implementations should move to :class:`.BackendV2`. Args: backend: The backend to run the primitive on. options: The options to control the default precision (``default_precision``), the operator grouping (``abelian_grouping``), and the random seed for the simulator (``seed_simulator``). z[The method `BackendEstimatorV2.__init__` will stop supporting inputs of type `BackendV1` ( z ) in the `backend` parameter in a future release no earlier than 2.0. `BackendV1` is deprecated and implementations should move to `BackendV2`.)category stacklevel)basistarget)r>Nr+) isinstancer warningswarnDeprecationWarning_backendr _optionsr from_backend basis_gatesrr?r _passmanager)selfr7r6r>opt1qs r-__init__zBackendEstimatorV2.__init__s('9--  M.&-...,       .5D**'***799 !.w77C gy ) ) >0uW^TTTEE0u===E'00r,returnr c|jS)zReturn the options)rErIs r-r6zBackendEstimatorV2.options }r,c|jS)z7Returns the backend which this sampler object based on.)rDrNs r-r7zBackendEstimatorV2.backendrOr,) precisionpubsIterable[EstimatorPubLike]rQ float | None(PrimitiveJob[PrimitiveResult[PubResult]]c |jjfd|D}||t|j|}||S)Nc:g|]}tj|Sr+)rcoerce).0pubrQs r- z*BackendEstimatorV2.run..s&LLL +C;;LLLr,)rEr"_validate_pubsr_run_submit)rIrRrQ coerced_pubsjobs ` r-runzBackendEstimatorV2.runse   7ILLLLtLLL  L)))49l33  r,list[EstimatorPub]c~t|D],\}}|jdkrtd|d|jdd-dS)NgzThe z/-th pub has precision less than or equal to 0 (z). z&But precision should be larger than 0.) enumeraterQ ValueError)rIrRirZs r-r\z!BackendEstimatorV2._validate_pubssdoo  FAs}## _1__S]___<$  r,PrimitiveResult[PubResult]ctt}tD]L\}}tt jd|jdzz }|||Mdgtz}| D]B\}}| fd|D|}t||D] \}} | ||< Ct|ddiS)Ng?r;c g|] }| Sr+r+)rYrfrRs r-r[z+BackendEstimatorV2._run..s)?)?)?a$q')?)?)?r,versionmetadata) rlistrdintmathceilrQappendlenitems _run_pubszipr) rIrRpub_dictrfrZshotsresultslst pub_results pub_results ` r-r]zBackendEstimatorV2._runst$$oo & &FAs # q(8"899::E UO " "1 % % % %&3t99$"..** ( (JE3..)?)?)?)?3)?)?)?GGK!$S+!6!6 ( ( :'  (w)Q@@@@r,rwrnlist[PubResult]c g}g}|D]F}||}||||jGt ||j||jj\}}t|} g} d} t||D]o\}}| t|jz} | | | | || | } | } | | || ||p| S)zFCompute results for pubs that all require the same value of ``shots``.)rwr%r) _preprocess_pubrqextendr1rrDrEr%rrurr_calc_expval_map_postprocess_pub)rIrRrwpreprocessed_data flat_circuitsrZdata run_resultrlcountsrxstartend expval_maps r-rtzBackendEstimatorV2._run_pubss-  0 0C'',,D  $ $T * * *   / / / /, 4=dmFb    H!,,T#455 P PIC#dm,,,C..veCi/@(5QT9BUVVJE NN400j$NN O O O Or,rZrr/c|j}|j}|j}|j}t jt j|t|}t j ||\}}tt} t j|jD]+} || } | |  || ,| ||| } t| ||S)aeConverts a pub into a list of bound circuits necessary to estimate all its observables. The circuits contain metadata explaining which bindings array index they are with respect to, and which measurement basis they are measuring. Args: pub: The pub to preprocess. Returns: The values ``(circuits, bc_param_ind, bc_obs)`` where ``circuits`` are the circuits to execute on the backend, ``bc_param_ind`` are indices of the pub's bindings array and ``bc_obs`` is the observables array, both broadcast to the shape of the pub. dtype)circuitr3parameter_valuesshapenpfromiterndindexobjectreshapebroadcast_arraysrsetupdate_bind_and_add_measurementsr/) rIrZrr3r param_shape param_indices bc_param_indbc_obs param_obs_mapindex param_indexbound_circuitss r-r~z"BackendEstimatorV2._preprocess_pubs+o /',  BJ{$;$;6JJJRRS^__ !2=+NN f#C(( Z!34 = =E&u-K + & - -fUm < < < <88BRTabb vFFFr,rdictrrc:|j}|j}tj|t}tj|t}tj|jD]r} || } || D]M\} } || | f\} }|| xx| | zz cc<|| xxtj| |dzzz cc<Ns|tj |z }t|||j}t||j ||j jdS)axComputes expectation values (evs) and standard errors (stds). The values are stored in arrays broadcast to the shape of the pub. Args: pub: The pub to postprocess. expval_map: The map data: The result data of the preprocessing. shots: The number of shots. Returns: The pub result. rg?)evsstdsr)target_precisionrwcircuit_metadatark)r2r3r zeros_liker!rrrsabssqrtrrrQrrl)rIrZrrrwrrr variancesrrpaulicoeffexpvalvariancerdata_bins r-rz#BackendEstimatorV2._postprocess_pub sH - !mL666M,e<<< Z!34 B BE&u-K &u 3 3 5 5 B B u#-k5.@#A E fun, %   BF5MMHcM$AA     B275>>)sSY??? $'M$'K$8    r,rrrrrdict[tuple[int, ...], set[str]]r0cg}|D]c\}}|||}tt|}||||} || d|S)aBind the given circuit against each parameter value set, and add necessary measurements to each. Args: circuit: The (possibly parametric) circuit of interest. parameter_values: An array of parameter value sets that can be applied to the circuit. param_obs_map: A mapping from locations in ``parameter_values`` to a sets of Pauli terms whose expectation values are required in those locations. Returns: A flat list of circuits sufficient to measure all Pauli terms in the ``param_obs_map`` values at the corresponding ``parameter_values`` location, where requisite book-keeping is stored as circuit metadata. )rsbindrsorted_create_measurement_circuitsr) rIrrrr1r pauli_strings bound_circuit meas_paulis new_circuitss r-rz-BackendEstimatorV2._bind_and_add_measurements/s(*7*=*=*?*? * * &K,11';GGM#F=$9$9::K<<{KL OOL ) ) ) )r,r list[Counts]rl6dict[tuple[tuple[int, ...], str], tuple[float, float]]ci}t||D]c\}}|d}|d}|d}t||\} } t|| | D]!\} } } | | f||| f<"d|S)aComputes the map of expectation values. Args: counts: The counts data. metadata: The metadata. Returns: The map of expectation values takes a pair of an index of the bindings array and a pauli string as a key and returns the expectation value of the pauli string with the the pub's circuit bound against the parameter value set in the index of the bindings array. orig_paulisrr)rurto_label)rIrrlrcountmetarrrexpvalsrrrrs r-rz#BackendEstimatorV2._calc_expval_mapNs"NP vx00 O OKE4}-K}-K}-K! meas_circuitindicespaulispreprocessed_circuits circuit_copyrcregs r-rz/BackendEstimatorV2._create_measurement_circuitsis"/1 = ) 3!11T1BB 3 3r}33CE::BM>L"'*F$*  ;$)++%R$*RRR,  % %f - - -  fd K K K$0$9L ! ! ( ( 6 6 6 6$$r,)r7r8r6r9)rLr )rLr8)rRrSrQrTrLrU)rRrb)rRrbrLrg)rRrbrwrnrLr|)rZrrLr/) rZrrrrr/rwrnrLr)rrrrrrrLr0)rrrlrrLr)rrrrrrrLr0)r&r'r(r)rKpropertyr6r7rar\r]rtr~rrrrr+r,r-r5r5KsW22p $ %1%1%1%1%1%1NXX NR      AAAA .GGGG># # # # J>6F%F%F%F%F%F%r,r5rrnrr ctj|j|j|jz}tj|sdg}t t|dtt|d|}t|D]_\}}|j|r7|j|r| || || ||`||fS)Nrq__c_)rarangerrranyrr rrrrdsdghmeasure)rr qubit_indicesrclbitrfs r-rrs Ie.//%'0ABM 6-  ! C((*;C rs[CB"""""" ######$$$$$$!!!!!!MMMMMMMMMM))))))1111111100000000 <<<<<<<<AAAAAAZZZZZZZZZZ!!!!!!MMMMMMMMMMMM444444222222''''''  &  H H H H H H H  Hd%d%d%d%d%d%d%d%N ''''''r,