iqm-benchmarks 2.9__py3-none-any.whl → 2.11__py3-none-any.whl

iqm/benchmarks/randomized_benchmarking/clifford_rb/clifford_rb.py

@@ -128,6 +128,12 @@ def clifford_rb_analysis(run: BenchmarkRunResult) -> BenchmarkAnalysisResult:
             "avg_gate_fidelity": {"value": fidelity.value, "uncertainty": fidelity.stderr},
         }
 
+        if len(qubits) == 1:
+            fidelity_native = rb_fit_results.params["fidelity_per_native_sqg"]
+            processed_results.update(
+                {"avg_native_gate_fidelity": {"value": fidelity_native.value, "uncertainty": fidelity_native.stderr}}
+            )
+
         dataset.attrs[qubits_idx].update(
             {
                 "decay_rate": {"value": popt["decay_rate"].value, "uncertainty": popt["decay_rate"].stderr},

iqm/benchmarks/randomized_benchmarking/interleaved_rb/interleaved_rb.py

@@ -75,6 +75,11 @@ def interleaved_rb_analysis(run: BenchmarkRunResult) -> BenchmarkAnalysisResult:
 
     interleaved_gate = dataset.attrs["interleaved_gate"]
     interleaved_gate_parameters = dataset.attrs["interleaved_gate_params"]
+    if interleaved_gate_parameters is None:
+        interleaved_gate_string = f"{interleaved_gate}"
+    else:
+        params_string = str(tuple(f"{x:.2f}" for x in interleaved_gate_parameters))
+        interleaved_gate_string = f"{interleaved_gate}{params_string}"
 
     simultaneous_fit = dataset.attrs["simultaneous_fit"]
 
@@ -138,6 +143,7 @@ def interleaved_rb_analysis(run: BenchmarkRunResult) -> BenchmarkAnalysisResult:
             [list_of_fidelities_clifford, list_of_fidelities_interleaved],
             "interleaved",
             simultaneous_fit,
+            interleaved_gate_string,
         )
         rb_fit_results = lmfit_minimizer(fit_parameters, fit_data, sequence_lengths, exponential_rb)
 
@@ -170,10 +176,31 @@ def interleaved_rb_analysis(run: BenchmarkRunResult) -> BenchmarkAnalysisResult:
                 "avg_gate_fidelity": {"value": fidelity.value, "uncertainty": fidelity.stderr},
             }
 
+            if len(qubits) == 1 and rb_type == "clifford":
+                fidelity_native = rb_fit_results.params["fidelity_per_native_sqg"]
+                processed_results[rb_type].update(
+                    {
+                        "avg_gate_fidelity_native": {
+                            "value": fidelity_native.value,
+                            "uncertainty": fidelity_native.stderr,
+                        }
+                    }
+                )
+            elif len(qubits) == 2 and rb_type == "clifford" and interleaved_gate_string == "CZGate":
+                fidelity_native_sqg = rb_fit_results.params["fidelity_per_native_sqg"]
+                processed_results[rb_type].update(
+                    {
+                        "avg_gate_fidelity_native_sqg": {
+                            "value": fidelity_native_sqg.value,
+                            "uncertainty": fidelity_native_sqg.stderr,
+                        }
+                    }
+                )
+
             observations.extend(
                 [
                     BenchmarkObservation(
-                        name=f"{key}_{rb_type}",
+                        name=f"{key}_{rb_type}" if "native" not in key else f"{key}",
                         identifier=BenchmarkObservationIdentifier(qubits),
                         value=values["value"],
                         uncertainty=values["uncertainty"],
@@ -206,12 +233,6 @@ def interleaved_rb_analysis(run: BenchmarkRunResult) -> BenchmarkAnalysisResult:
         obs_dict.update({qubits_idx: processed_results})
 
         # Generate decay plots
-        if interleaved_gate_parameters is None:
-            interleaved_gate_string = f"{interleaved_gate}"
-        else:
-            params_string = str(tuple(f"{x:.2f}" for x in interleaved_gate_parameters))
-            interleaved_gate_string = f"{interleaved_gate}{params_string}"
-
         fig_name, fig = plot_rb_decay(
             "irb",
             [qubits],

iqm/benchmarks/randomized_benchmarking/randomized_benchmarking_common.py

@@ -89,6 +89,7 @@ def fit_decay_lmfit(
     data: List[List[float]] | List[List[List[float]]],
     rb_identifier: str,
     simultaneous_fit_vars: Optional[List[str]] = None,
+    interleaved_gate_str: Optional[str] = None,
 ) -> Tuple[np.ndarray, Parameters]:
     """Perform a fitting routine for 0th-order (Ap^m+B) RB using lmfit
 
@@ -98,6 +99,7 @@ def fit_decay_lmfit(
         data (List[List[float]] | List[List[List[float]]]): the data to be fitted
         rb_identifier (str): the RB identifier, either "stdrb", "irb" or "mrb"
         simultaneous_fit_vars (List[str], optional): the list of variables used to fit simultaneously
+        interleaved_gate_str (Optional[str]): the name of the interleaved gate in IRB
     Returns:
         A tuple of fitting data (list of lists of average fidelities or polarizations) and MRB fit parameters
     """
@@ -125,6 +127,12 @@ def fit_decay_lmfit(
             )
             params.add(f"p_rb", expr=f"1-depolarization_probability_{1}")
             params.add(f"fidelity_per_clifford", expr=f"p_rb + (1 - p_rb) / (2**{n_qubits})")
+            if n_qubits == 1:
+                # The construction of the 1Q and 2Q Clifford gate dictionaries in "generate_2qubit_cliffords.ipynb"
+                # is based on the reference below; thus the expected amount of 1.875 native gates per Clifford
+                # Native gate set being {I, X(pi/2), X(pi), X(-pi/2), Y(pi/2), Y(pi), Y(-pi)} with X(a)=r(a,0), Y(a)=r(a,pi/2)
+                # Ref: Barends et al., Nature 508, 500-503 (2014); Eq.(S3) of arXiv:1402.4848 [quant-ph]
+                params.add("fidelity_per_native_sqg", expr=f"1 - (1 - (p_rb + (1 - p_rb) / (2**{n_qubits})))/1.875")
         else:
             params = create_multi_dataset_params(
                 func, fit_data, initial_guesses=None, constraints=constraints, simultaneously_fit_vars=None
@@ -140,10 +148,24 @@ def fit_decay_lmfit(
             constraints=None,
             simultaneously_fit_vars=simultaneous_fit_vars,
         )
-        params.add(f"p_rb", expr=f"1-depolarization_probability_{1}")
-        params.add(f"fidelity_per_clifford", expr=f"p_rb + (1 - p_rb) / (2**{n_qubits})")
-        params.add(f"p_irb", expr=f"1-depolarization_probability_{2}")
-        params.add(f"interleaved_fidelity", expr=f"p_irb / p_rb + (1 - p_irb / p_rb) / (2**{n_qubits})")
+        params.add(f"p_rb", expr=f"1.0 - depolarization_probability_{1}")
+        params.add(f"fidelity_per_clifford", expr=f"p_rb + (1.0 - p_rb) / (2.0**{n_qubits})")
+        params.add(f"p_irb", expr=f"1.0 - depolarization_probability_{2}")
+        params.add(f"interleaved_fidelity", expr=f"p_irb / p_rb + (1.0 - p_irb / p_rb) / (2.0**{n_qubits})")
+        if n_qubits == 1:
+            # The construction of the 1Q and 2Q Clifford gate dictionaries in "generate_2qubit_cliffords.ipynb"
+            # is based on the reference below; thus the expected amount of 1.875 native gates per Clifford
+            # Native gate set being {I, X(pi/2), X(pi), X(-pi/2), Y(pi/2), Y(pi), Y(-pi)} with X(a)=r(a,0), Y(a)=r(a,pi/2)
+            # Ref: Barends et al., Nature 508, 500-503 (2014); Eq.(S3) of arXiv:1402.4848 [quant-ph]
+            # For IRB it may be used as proxy to how well (or bad) the assumption of uniform fidelity in native sqg gates holds.
+            params.add("fidelity_per_native_sqg", expr=f"1.0 - (1.0 - fidelity_per_clifford)/1.875")
+        elif n_qubits == 2 and interleaved_gate_str == "CZGate":
+            # Here similarly, we may use Eq.(S8 - S11) of arXiv:1402.4848 [quant-ph]
+            params.add(f"fidelity_clifford_and_interleaved", expr=f"p_irb + (1.0 - p_irb) / (2.0 ** {n_qubits})")
+            params.add(
+                "fidelity_per_native_sqg",
+                expr=f"(1.0 / 33.0) * (4.0 * fidelity_clifford_and_interleaved - 10.0 * interleaved_fidelity + 39.0)",
+            )
 
     return fit_data, params
 
@@ -578,6 +600,8 @@ def plot_rb_decay(
     stddevs_from_mean = {}
     fidelity_value = {}
     fidelity_stderr = {}
+    fidelity_native1q_value = {}
+    fidelity_native1q_stderr = {}
     decay_rate = {}
     offset = {}
     amplitude = {}
@@ -618,6 +642,14 @@ def plot_rb_decay(
                 str(q): dataset.attrs[q_idx]["avg_fidelities_stderr"]
                 for q_idx, q in enumerate(qubits_array, qubits_index)
             }
+            fidelity_native1q_value[identifier] = {
+                str(q): observations[q_idx]["avg_native_gate_fidelity"]["value"] if len(q) == 1 else np.nan
+                for q_idx, q in enumerate(qubits_array, qubits_index)
+            }
+            fidelity_native1q_stderr[identifier] = {
+                str(q): observations[q_idx]["avg_native_gate_fidelity"]["uncertainty"] if len(q) == 1 else np.nan
+                for q_idx, q in enumerate(qubits_array, qubits_index)
+            }
         # These are common to both MRB and standard Clifford
         fidelity_value[identifier] = {
             str(q): observations[q_idx]["avg_gate_fidelity"]["value"]
@@ -636,7 +668,7 @@ def plot_rb_decay(
         amplitude[identifier] = {
             str(q): dataset.attrs[q_idx]["fit_amplitude"]["value"] for q_idx, q in enumerate(qubits_array, qubits_index)
         }
-    else:  # id MRB
+    else:  # id IRB
         rb_type_keys = list(observations[0].keys())
         colors = [cmap(i) for i in np.linspace(start=1, stop=0, num=len(rb_type_keys)).tolist()]
         for rb_type in rb_type_keys:
@@ -791,12 +823,12 @@ def plot_rb_decay(
             if identifier == "mrb":
                 plot_label = fr"$\overline{{F}}_{{MRB}} (n={len(qubits)})$ = {100.0 * fidelity_value[key][str(qubits)]:.2f} +/- {100.0 * fidelity_stderr[key][str(qubits)]:.2f} (%)"
             elif key == "interleaved":
-                plot_label = fr"$\overline{{F}}_{{{interleaved_gate}}} ({qubits})$ = {100.0 * fidelity_value[key][str(qubits)]:.2f} +/- {100.0 * fidelity_stderr[key][str(qubits)]:.2f} (%)"
-            else:
-                print(fidelity_value)
-                print(qubits)
-                print(key)
-                plot_label = fr"$\overline{{F}}_{{CRB}}$ = {100.0 * fidelity_value[key][str(qubits)]:.2f} +/- {100.0 * fidelity_stderr[key][str(qubits)]:.2f} (%)"
+                plot_label = fr"$\overline{{F}}_{{{interleaved_gate}}}$ = {100.0 * fidelity_value[key][str(qubits)]:.2f} +/- {100.0 * fidelity_stderr[key][str(qubits)]:.2f} (%)"
+            else:  # if id is "clifford"
+                if len(qubits) == 1 and identifier != "irb":
+                    plot_label = fr"$\overline{{F}}_{{native\_sqg}} \simeq$ {100.0 * fidelity_native1q_value[key][str(qubits)]:.2f} +/- {100.0 * fidelity_native1q_stderr[key][str(qubits)]:.2f} (%)"
+                else:
+                    plot_label = fr"$\overline{{F}}_{{CRB}}$ = {100.0 * fidelity_value[key][str(qubits)]:.2f} +/- {100.0 * fidelity_stderr[key][str(qubits)]:.2f} (%)"
 
             ax.plot(
                 x_linspace,

{iqm_benchmarks-2.9.dist-info → iqm_benchmarks-2.11.dist-info}/METADATA

@@ -1,6 +1,6 @@
 Metadata-Version: 2.2
 Name: iqm-benchmarks
-Version: 2.9
+Version: 2.11
 Summary: A package for implementation of Quantum Characterization, Verification and Validation (QCVV) techniques on IQM's hardware at gate level abstraction
 Author-email: IQM Finland Oy <developers@meetiqm.com>, Aniket Rath <aniket.rath@meetiqm.com>, Jami Rönkkö <jami@meetiqm.com>, Pedro Figueroa Romero <pedro.romero@meetiqm.com>, Vicente Pina Canelles <vicente.pina@meetiqm.com>, Raphael Brieger <raphael.brieger@meetiqm.com>, Stefan Seegerer <stefan.seegerer@meetiqm.com>, Miikka Koistinen <miikka@meetiqm.com>
 Project-URL: Homepage, https://github.com/iqm-finland/iqm-benchmarks

{iqm_benchmarks-2.9.dist-info → iqm_benchmarks-2.11.dist-info}/RECORD

@@ -19,11 +19,11 @@ iqm/benchmarks/randomized_benchmarking/__init__.py,sha256=IkKo-7zUChxZZd3my_csQC
 iqm/benchmarks/randomized_benchmarking/clifford_1q.pkl,sha256=vvSd0pRWxtzyirohO9yf_58mjevkc2-pbuWIEb-4gaw,46928
 iqm/benchmarks/randomized_benchmarking/clifford_2q.pkl,sha256=ZipqU3crPhz2T35qGFgB4GvMyoi_7pnu8NqW5ZP8NXg,90707258
 iqm/benchmarks/randomized_benchmarking/multi_lmfit.py,sha256=Se1ygR4mXn_2_P82Ch31KBnCmY-g_A9NKzE9Ir8nEvw,3247
-iqm/benchmarks/randomized_benchmarking/randomized_benchmarking_common.py,sha256=7vp2BEXSYz3Q7u_42qPKlQtdWPAI-a90tmtlhThWLjg,38809
+iqm/benchmarks/randomized_benchmarking/randomized_benchmarking_common.py,sha256=xMiQyikGEVj3zoFQfVMS50K15MszArc0kiBWvZ7gWL8,41492
 iqm/benchmarks/randomized_benchmarking/clifford_rb/__init__.py,sha256=bTDA156LAl7OLGcMec--1nzDrV1XpPRVq3CquTmucgE,677
-iqm/benchmarks/randomized_benchmarking/clifford_rb/clifford_rb.py,sha256=q8bHcrVr8EjTuLeOdg3RoeoLaZhmieCd09IvH-tSvYs,18038
+iqm/benchmarks/randomized_benchmarking/clifford_rb/clifford_rb.py,sha256=v8GDsEC3JscZVJXc6ZqfJaaSb1LocdFTSeOcwHxp50Y,18317
 iqm/benchmarks/randomized_benchmarking/interleaved_rb/__init__.py,sha256=sq6MgN_hwlpkOj10vyCU4e6eKSX-oLcF2L9na6W2Gt4,681
-iqm/benchmarks/randomized_benchmarking/interleaved_rb/interleaved_rb.py,sha256=PlHn1J8VPaJF5csNH8jxcifz_MdisOEPU54kU-FYoLY,26920
+iqm/benchmarks/randomized_benchmarking/interleaved_rb/interleaved_rb.py,sha256=4hBjIXxqy18WKC8hgpCU6T6Vv2cvJu4AImj7QQDqI1k,27956
 iqm/benchmarks/randomized_benchmarking/mirror_rb/__init__.py,sha256=ZekEqI_89nXzGO1vjM-b5Uwwicy59M4fYHXfA-f0MIg,674
 iqm/benchmarks/randomized_benchmarking/mirror_rb/mirror_rb.py,sha256=pOEnEtuJ0rMnTq5nZ1obsQ6LOIVlgFcgNGc8CULWz6E,33987
 mGST/LICENSE,sha256=TtHNq55cUcbglb7uhVudeBLUh_qPdUoAEvU0BBwFz-k,1098
@@ -36,8 +36,8 @@ mGST/optimization.py,sha256=YHwkzIkYvsZOPjclR-BCQWh24jeqjuXp0BB0WX5Lwow,10559
 mGST/qiskit_interface.py,sha256=L4H-4SdhP_bjSFFvpQoF1E7EyGbIJ_CI_y4a7_YEwmU,10102
 mGST/reporting/figure_gen.py,sha256=6Xd8vwfy09hLY1YbJY6TRevuMsQSU4MsWqemly3ZO0I,12970
 mGST/reporting/reporting.py,sha256=We1cccz9BKbITYcSlZHdmBGdjMWAa1xNZe5tKP-yh_E,26004
-iqm_benchmarks-2.9.dist-info/LICENSE,sha256=2Ncb40-hqkTil78RPv3-YiJfKaJ8te9USJgliKqIdSY,11558
-iqm_benchmarks-2.9.dist-info/METADATA,sha256=9ByGfYfFNY1OMBRjDm-8GQMWUGQSaBJe3a3AmSAFjUU,10170
-iqm_benchmarks-2.9.dist-info/WHEEL,sha256=In9FTNxeP60KnTkGw7wk6mJPYd_dQSjEZmXdBdMCI-8,91
-iqm_benchmarks-2.9.dist-info/top_level.txt,sha256=3G23Z-1LGf-IOzTCUl6QwWqiQ3USz25Zt90Ihq192to,9
-iqm_benchmarks-2.9.dist-info/RECORD,,
+iqm_benchmarks-2.11.dist-info/LICENSE,sha256=2Ncb40-hqkTil78RPv3-YiJfKaJ8te9USJgliKqIdSY,11558
+iqm_benchmarks-2.11.dist-info/METADATA,sha256=vijW3YgqFM3XIVA7uHSImcjeRR5Sr2b9a75a_zetmrQ,10171
+iqm_benchmarks-2.11.dist-info/WHEEL,sha256=In9FTNxeP60KnTkGw7wk6mJPYd_dQSjEZmXdBdMCI-8,91
+iqm_benchmarks-2.11.dist-info/top_level.txt,sha256=3G23Z-1LGf-IOzTCUl6QwWqiQ3USz25Zt90Ihq192to,9
+iqm_benchmarks-2.11.dist-info/RECORD,,