iqm-benchmarks 2.3__py3-none-any.whl → 2.5__py3-none-any.whl

iqm/benchmarks/randomized_benchmarking/mirror_rb/mirror_rb.py

@@ -10,8 +10,8 @@ import warnings
 
 import numpy as np
 from qiskit import transpile
-from qiskit.quantum_info import random_clifford, random_pauli
-from qiskit_aer import Aer
+from qiskit.quantum_info import Clifford, random_clifford, random_pauli
+from qiskit_aer import Aer, AerSimulator
 from scipy.spatial.distance import hamming
 import xarray as xr
 
@@ -262,9 +262,13 @@ def generate_pauli_dressed_mrb_circuits(
     pauli_dressed_circuits_untranspiled: List[QuantumCircuit] = []
     pauli_dressed_circuits_transpiled: List[QuantumCircuit] = []
 
+    sim_method = "stabilizer"
+    simulator = AerSimulator(method=sim_method)
+
     for _ in range(pauli_samples_per_circ):
         # Initialize the quantum circuit object
         circ = QuantumCircuit(num_qubits)
+        circ_untransp = QuantumCircuit(num_qubits)
         # Sample all the random Paulis
         paulis = [random_pauli(num_qubits) for _ in range(depth + 1)]
 
@@ -282,6 +286,7 @@ def generate_pauli_dressed_mrb_circuits(
             )
             circ.barrier()
             circ.compose(cycle_layers[k], inplace=True)
+            circ_untransp.compose(cycle_layers[k], inplace=True)
             circ.barrier()
 
         # Apply middle Pauli
@@ -307,9 +312,6 @@ def generate_pauli_dressed_mrb_circuits(
         for i in range(num_qubits):
             circ.compose(clifford_layer[i].to_instruction().inverse(), qubits=[i], inplace=True)
 
-        # Add measurements
-        circ.measure_all()
-
         # Transpile to backend - no optimize SQG should be used!
         if isinstance(backend_arg, str):
             retrieved_backend = get_iqm_backend(backend_arg)
@@ -317,6 +319,13 @@ def generate_pauli_dressed_mrb_circuits(
             assert isinstance(backend_arg, IQMBackendBase)
             retrieved_backend = backend_arg
 
+        circ_untransp = circ.copy()
+        # Add measurements to untranspiled - after!
+        circ_untranspiled = transpile(Clifford(circ_untransp).to_circuit(), simulator)
+        circ_untranspiled.measure_all()
+
+        # Add measurements to transpiled - before!
+        circ.measure_all()
         circ_transpiled = transpile(
             circ,
             backend=retrieved_backend,
@@ -325,7 +334,7 @@ def generate_pauli_dressed_mrb_circuits(
             routing_method=routing_method,
         )
 
-        pauli_dressed_circuits_untranspiled.append(circ)
+        pauli_dressed_circuits_untranspiled.append(circ_untranspiled)
         pauli_dressed_circuits_transpiled.append(circ_transpiled)
 
     # Store the circuit
@@ -418,14 +427,14 @@ def mrb_analysis(run: BenchmarkRunResult) -> BenchmarkAnalysisResult:
     observations = {}
     dataset = run.dataset.copy(deep=True)
 
-    shots = dataset.attrs["shots"]
+    # shots = dataset.attrs["shots"]
     num_circuit_samples = dataset.attrs["num_circuit_samples"]
     num_pauli_samples = dataset.attrs["num_pauli_samples"]
 
     density_2q_gates = dataset.attrs["density_2q_gates"]
     two_qubit_gate_ensemble = dataset.attrs["two_qubit_gate_ensemble"]
 
-    max_gates_per_batch = dataset.attrs["max_gates_per_batch"]
+    # max_gates_per_batch = dataset.attrs["max_gates_per_batch"]
 
     # Analyze the results for each qubit layout of the experiment dataset
     qubits_array = dataset.attrs["qubits_array"]
@@ -445,8 +454,9 @@ def mrb_analysis(run: BenchmarkRunResult) -> BenchmarkAnalysisResult:
     else:
         assigned_mrb_depths = {str(qubits_array[i]): [2 * m for m in depths_array[i]] for i in range(len(depths_array))}
 
-    transpiled_circuits = run.circuits["transpiled_circuits"]
-    simulator = Aer.get_backend("qasm_simulator")
+    mrb_sim_circuits = run.circuits["untranspiled_circuits"]
+    sim_method = "stabilizer"
+    simulator = AerSimulator(method=sim_method)  # Aer.get_backend("stabilizer")
 
     all_noisy_counts: Dict[str, Dict[int, List[Dict[str, int]]]] = {}
     all_noiseless_counts: Dict[str, Dict[int, List[Dict[str, int]]]] = {}
@@ -465,22 +475,12 @@ def mrb_analysis(run: BenchmarkRunResult) -> BenchmarkAnalysisResult:
             all_noisy_counts[str(qubits)][depth] = xrvariable_to_counts(
                 dataset, identifier, num_circuit_samples * num_pauli_samples
             )
-
             qcvv_logger.info(f"Depth {depth}")
-            # Execute the quantum circuits on the simulated, ideal backend
-            # pylint: disable=unbalanced-tuple-unpacking
-            all_noiseless_jobs, _ = submit_execute(
-                {tuple(qubits): transpiled_circuits[f"{str(qubits)}_depth_{str(depth)}"].circuits},
-                simulator,
-                shots,
-                calset_id=None,
-                max_gates_per_batch=max_gates_per_batch,
-            )
 
-            # Retrieve counts
-            all_noiseless_counts[str(qubits)][depth], time_retrieve_noiseless[str(qubits)][depth] = retrieve_all_counts(
-                all_noiseless_jobs
-            )
+            mrb_circs = mrb_sim_circuits[f"{str(qubits)}_depth_{str(depth)}"].circuits
+
+            qcvv_logger.info("Getting simulation counts")
+            all_noiseless_counts[str(qubits)][depth] = simulator.run(mrb_circs).result().get_counts()
 
             # Compute polarizations for the current depth
             polarizations[depth] = compute_polarizations(
@@ -695,7 +695,7 @@ class MirrorRandomizedBenchmarking(Benchmark):
             mrb_untranspiled_circuits_lists: Dict[int, List[QuantumCircuit]] = {}
             time_circuit_generation[str(qubits)] = 0
             for depth in assigned_mrb_depths[str(qubits)]:
-                qcvv_logger.info(f"Depth {depth}")
+                qcvv_logger.info(f"Depth {depth} - Generating all circuits")
                 mrb_circuits[depth], elapsed_time = generate_fixed_depth_mrb_circuits(
                     qubits,
                     self.num_circuit_samples,

{iqm_benchmarks-2.3.dist-info → iqm_benchmarks-2.5.dist-info}/METADATA

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: iqm-benchmarks
-Version: 2.3
+Version: 2.5
 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>
 Project-URL: Homepage, https://github.com/iqm-finland/iqm-benchmarks

{iqm_benchmarks-2.3.dist-info → iqm_benchmarks-2.5.dist-info}/RECORD

@@ -1,7 +1,6 @@
 iqm/benchmarks/__init__.py,sha256=sAS6OKctANSzgVUEUWugQs9mcKwG9KOTBEsbSYGxtPc,2201
 iqm/benchmarks/benchmark.py,sha256=SGhBcSxLPUu-cVXAjG4Db2TRobFCRBYoE1NtTDK1lJg,4432
 iqm/benchmarks/benchmark_definition.py,sha256=AZkvANrf0_0glbq_P_uo_YqbBU9IZa2gJlMVz6qT6VU,10500
-iqm/benchmarks/benchmark_experiment.py,sha256=ynbsQS92v6soIgkDEdm3ygNO7f82-dhg7qB1l8jvSFM,6614
 iqm/benchmarks/circuit_containers.py,sha256=anEtZEsodYqOX-34oZRmuKGeEpp_VfgG5045Mz4-4hI,7562
 iqm/benchmarks/logging_config.py,sha256=U7olP5Kr75AcLJqNODf9VBhJLVqIvA4AYR6J39D5rww,1052
 iqm/benchmarks/readout_mitigation.py,sha256=7FlbSH-RJTtQuRYLChwkQV_vBv0ZfMQTH519cAbyxQ4,12252
@@ -12,7 +11,7 @@ iqm/benchmarks/compressive_gst/gst_analysis.py,sha256=qaF9zDIadPGe9I0l_SC3EfnODZ
 iqm/benchmarks/entanglement/__init__.py,sha256=9T7prOwqMmFWdb4t6ETAHZXKK5o6FvU2DvVb6WhNi-U,682
 iqm/benchmarks/entanglement/ghz.py,sha256=qg6mw1RydosKdc2KoPi6M3q0Ek7eMBPIkWOJC2ZDx8M,39446
 iqm/benchmarks/optimization/__init__.py,sha256=_ajW_OibYLCtzU5AUv5c2zuuVYn8ZNeZUcUUSIGt51M,747
-iqm/benchmarks/optimization/qscore.py,sha256=pTDRrNI2kkujqb5IoaPTiq9MS8OZd9M8Gv0G2Jn3Ka0,28078
+iqm/benchmarks/optimization/qscore.py,sha256=6I13YbFvFq1RcX4mYTw6S4ALb4Ix9t7gjlRzZbaTARM,34038
 iqm/benchmarks/quantum_volume/__init__.py,sha256=i-Q4SpDWELBw7frXnxm1j4wJRcxbIyrS5uEK_v06YHo,951
 iqm/benchmarks/quantum_volume/clops.py,sha256=j6BPEj1rKBAHigox7nrvaTLyb4iCrHadBl2d1yiETDA,30956
 iqm/benchmarks/quantum_volume/quantum_volume.py,sha256=PPg1kZA5Lx9b24iFVgP_VcgymI6CsD5LfB5Ppfk8TBM,36814
@@ -26,7 +25,7 @@ iqm/benchmarks/randomized_benchmarking/clifford_rb/clifford_rb.py,sha256=q8bHcrV
 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/mirror_rb/__init__.py,sha256=ZekEqI_89nXzGO1vjM-b5Uwwicy59M4fYHXfA-f0MIg,674
-iqm/benchmarks/randomized_benchmarking/mirror_rb/mirror_rb.py,sha256=kYlcH0uFdFZmZhkPPEliUTftfDTpefjFdlTgDDo_NXM,33819
+iqm/benchmarks/randomized_benchmarking/mirror_rb/mirror_rb.py,sha256=AmPAp80wWc_6ao4OMZBUNGX0_WlSnz9utN4436ZcjJ0,33990
 mGST/LICENSE,sha256=TtHNq55cUcbglb7uhVudeBLUh_qPdUoAEvU0BBwFz-k,1098
 mGST/README.md,sha256=v_5kw253csHF4-RfE-44KqFmBXIsSMRmOtN0AUPrRxE,5050
 mGST/additional_fns.py,sha256=_SEJ10FRNM7_CroysT8hCLZTfpm6ZhEIDCY5zPTnhjo,31390
@@ -37,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=-XBy3OCJIMOsA8cApwKjhVKBwnjSAoxm-voHNbRWytM,25803
-iqm_benchmarks-2.3.dist-info/LICENSE,sha256=2Ncb40-hqkTil78RPv3-YiJfKaJ8te9USJgliKqIdSY,11558
-iqm_benchmarks-2.3.dist-info/METADATA,sha256=9X5g_7kRRtPKs581ILEaU7C00SDdWERCXjoa_OCi4Jw,9103
-iqm_benchmarks-2.3.dist-info/WHEEL,sha256=PZUExdf71Ui_so67QXpySuHtCi3-J3wvF4ORK6k_S8U,91
-iqm_benchmarks-2.3.dist-info/top_level.txt,sha256=3G23Z-1LGf-IOzTCUl6QwWqiQ3USz25Zt90Ihq192to,9
-iqm_benchmarks-2.3.dist-info/RECORD,,
+iqm_benchmarks-2.5.dist-info/LICENSE,sha256=2Ncb40-hqkTil78RPv3-YiJfKaJ8te9USJgliKqIdSY,11558
+iqm_benchmarks-2.5.dist-info/METADATA,sha256=c3KkSGEMw7ux1jcjbN4k-JU-hNfQSLcg8TOfQZb2npo,9103
+iqm_benchmarks-2.5.dist-info/WHEEL,sha256=A3WOREP4zgxI0fKrHUG8DC8013e3dK3n7a6HDbcEIwE,91
+iqm_benchmarks-2.5.dist-info/top_level.txt,sha256=3G23Z-1LGf-IOzTCUl6QwWqiQ3USz25Zt90Ihq192to,9
+iqm_benchmarks-2.5.dist-info/RECORD,,

{iqm_benchmarks-2.3.dist-info → iqm_benchmarks-2.5.dist-info}/WHEEL

@@ -1,5 +1,5 @@
 Wheel-Version: 1.0
-Generator: setuptools (75.6.0)
+Generator: setuptools (75.7.0)
 Root-Is-Purelib: true
 Tag: py3-none-any
 

iqm/benchmarks/benchmark_experiment.py

@@ -1,163 +0,0 @@
-# Copyright 2024 IQM Benchmarks developers
-#
-# Licensed under the Apache License, Version 2.0 (the "License");
-# you may not use this file except in compliance with the License.
-# You may obtain a copy of the License at
-#
-#     http://www.apache.org/licenses/LICENSE-2.0
-#
-# Unless required by applicable law or agreed to in writing, software
-# distributed under the License is distributed on an "AS IS" BASIS,
-# WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
-# See the License for the specific language governing permissions and
-# limitations under the License.
-
-"""
-Generic Benchmark Experiment class
-"""
-
-from copy import deepcopy
-from json import dump
-from pathlib import Path
-import pickle
-from time import strftime
-from typing import List, Optional, OrderedDict, Union
-
-from iqm.benchmarks.benchmark import BenchmarkBase, BenchmarkConfigurationBase
-from iqm.benchmarks.logging_config import qcvv_logger
-from iqm.benchmarks.utils import get_iqm_backend
-from iqm.qiskit_iqm.iqm_backend import IQMBackendBase
-
-
-class BenchmarkExperiment:
-    """
-    A Benchmark Experiment wraps the execution of one or more benchmarks, checks their requirements are met, stores the
-    execution results and plots relevant figures.
-    """
-
-    def __init__(
-        self,
-        backend: Union[str, IQMBackendBase],
-        benchmark_configurations: List[BenchmarkConfigurationBase],
-        device_id: Optional[str] = None,
-    ):
-        """Construct the BenchmarkExperiment class.
-
-        Args:
-            backend (str | IQMBackendBase): the backend to execute the benchmarks on
-            benchmark_configurations (List[BenchmarkConfigurationBase]): the configuration(s) of the benchmark(s)
-            device_id (Optional[str], optional): the identifier of the device. Defaults to None.
-
-        Raises:
-            ValueError: backend not supported. Try 'garnet' or 'iqmfakeadonis'
-        """
-        self.timestamp = strftime("%Y%m%d-%H%M%S")
-
-        if isinstance(backend, str):
-            self.backend = get_iqm_backend(backend)
-        else:
-            assert isinstance(backend, IQMBackendBase)
-            self.backend = backend
-
-        self.device_id = device_id if device_id is not None else self.backend.name
-
-        benchmarks: OrderedDict[str, BenchmarkBase] = OrderedDict(
-            (config.benchmark.name(), config.benchmark(self.backend, config)) for config in benchmark_configurations
-        )
-
-        for benchmark in benchmarks.values():
-            benchmarks_copy = deepcopy(benchmarks)
-            benchmarks_copy = benchmark.check_requirements(benchmarks_copy)
-        self.benchmarks = benchmarks_copy
-
-    def run_experiment(self) -> None:
-        """Run the Benchmark experiment, and store the configuration, raw data, results and figures."""
-
-        for name, benchmark in self.benchmarks.items():
-            qcvv_logger.info("\nNow executing " + name)
-            # Create the directory for results
-            results_dir = f"Outputs/{self.device_id}/{self.timestamp}/{name}/"
-            Path(results_dir).mkdir(parents=True, exist_ok=True)
-
-            # Execute the current benchmark
-            benchmark.generate_requirements(self.benchmarks)
-            benchmark.execute_full_benchmark()
-
-            # Create configuration JSON file
-            with open(
-                f"{results_dir}{self.device_id}_{self.timestamp}_{name}_configuration.json",
-                "w",
-                encoding="utf-8",
-            ) as f_configuration:
-                dump(benchmark.serializable_configuration.model_dump(), f_configuration)
-
-            # Create untranspiled circuit files
-            if benchmark.untranspiled_circuits:
-                for key_qubits in benchmark.untranspiled_circuits.keys():
-                    for key_depth in benchmark.untranspiled_circuits[key_qubits].keys():
-                        with open(
-                            f"{results_dir}{self.device_id}_{self.timestamp}_{name}_qubits_{key_qubits}_depth_{key_depth}_untranspiled.pkl",
-                            "wb",
-                        ) as f_circuits:
-                            pickle.dump(
-                                benchmark.untranspiled_circuits[key_qubits][key_depth],
-                                f_circuits,
-                            )
-
-            # Create transpiled circuit files
-            for key_qubits in benchmark.transpiled_circuits.keys():
-                for key_depth in benchmark.transpiled_circuits[key_qubits].keys():
-                    with open(
-                        f"{results_dir}{self.device_id}_{self.timestamp}_{name}_qubits_{key_qubits}_depth_{key_depth}_transpiled.pkl",
-                        "wb",
-                    ) as f_circuits:
-                        pickle.dump(
-                            benchmark.transpiled_circuits[key_qubits][key_depth],
-                            f_circuits,
-                        )
-
-            # Create raw result pickle files
-            with open(
-                f"{results_dir}{self.device_id}_{self.timestamp}_{name}_raw_results.pkl",
-                "wb",
-            ) as f_raw_results:
-                pickle.dump(
-                    benchmark.raw_results,
-                    f_raw_results,
-                )
-
-            # Create raw data JSON file
-            with open(
-                f"{results_dir}{self.device_id}_{self.timestamp}_{name}_raw_data.json",
-                "w",
-                encoding="utf-8",
-            ) as f_raw_data:
-                dump(benchmark.raw_data, f_raw_data)
-
-            # Create job metadata JSON file
-            with open(
-                f"{results_dir}{self.device_id}_{self.timestamp}_{name}_job_metadata.json",
-                "w",
-                encoding="utf-8",
-            ) as f_job_metadata:
-                dump(benchmark.job_meta, f_job_metadata)
-
-            # Create results JSON file
-            with open(
-                f"{results_dir}{self.device_id}_{self.timestamp}_{name}_results.json",
-                "w",
-                encoding="utf-8",
-            ) as f_results:
-                dump(benchmark.results, f_results)
-
-            # Create figures
-            Path(f"{results_dir}figures/").mkdir(parents=True, exist_ok=True)
-            for fig_name, fig in benchmark.figures.items():
-                fig.savefig(
-                    f"{results_dir}figures/{self.device_id}_{self.timestamp}_{name}_{fig_name}",
-                    dpi=250,
-                    bbox_inches="tight",
-                )
-
-            # Save benchmark
-            self.benchmarks[benchmark.name()] = benchmark