cirq-core 1.7.0.dev20250824200454__py3-none-any.whl → 1.7.0.dev20250825201352__py3-none-any.whl

cirq/_version.py

@@ -28,4 +28,4 @@ if sys.version_info < (3, 11, 0):  # pragma: no cover
         'of Cirq (e.g. "python -m pip install cirq==1.5.0")'
     )
 
-__version__ = "1.7.0.dev20250824200454"
+__version__ = "1.7.0.dev20250825201352"

cirq/_version_test.py

@@ -3,4 +3,4 @@ import cirq
 
 
 def test_version() -> None:
-    assert cirq.__version__ == "1.7.0.dev20250824200454"
+    assert cirq.__version__ == "1.7.0.dev20250825201352"

cirq/experiments/qubit_characterizations.py

@@ -18,7 +18,7 @@ import dataclasses
 import functools
 import itertools
 import uuid
-from typing import Any, cast, Iterator, Mapping, Sequence, TYPE_CHECKING
+from typing import Any, cast, Iterator, Mapping, Optional, Sequence, TYPE_CHECKING
 
 import attrs
 import numpy as np
@@ -73,7 +73,12 @@ class Cliffords:
 class RandomizedBenchMarkResult:
     """Results from a randomized benchmarking experiment."""
 
-    def __init__(self, num_cliffords: Sequence[int], ground_state_probabilities: Sequence[float]):
+    def __init__(
+        self,
+        num_cliffords: Sequence[int],
+        ground_state_probabilities: Sequence[float],
+        ground_state_probabilities_std: Optional[Sequence[float]] | None = None,
+    ):
         """Inits RandomizedBenchMarkResult.
 
         Args:
@@ -81,9 +86,20 @@ class RandomizedBenchMarkResult:
                 study.
             ground_state_probabilities: The corresponding average ground state
                 probabilities.
+            ground_state_probabilities_std: The standard deviation of the probabilities.
         """
         self._num_cfds_seq = num_cliffords
         self._gnd_state_probs = ground_state_probabilities
+        if ground_state_probabilities_std is None or np.all(
+            np.isclose(ground_state_probabilities_std, 0)
+        ):
+            self._gnd_state_probs_std = None
+        else:
+            self._gnd_state_probs_std = np.array(ground_state_probabilities_std)
+            zeros = np.isclose(self._gnd_state_probs_std, 0)
+            self._gnd_state_probs_std[zeros] = self._gnd_state_probs_std[
+                np.logical_not(zeros)
+            ].min()
 
     @property
     def data(self) -> Sequence[tuple[int, float]]:
@@ -142,6 +158,7 @@ class RandomizedBenchMarkResult:
             f=exp_fit,
             xdata=self._num_cfds_seq,
             ydata=self._gnd_state_probs,
+            sigma=self._gnd_state_probs_std,
             p0=[0.5, 0.5, 1.0 - 1e-3],
             bounds=([0, -1, 0], [1, 1, 1]),
         )
@@ -534,6 +551,7 @@ def parallel_single_qubit_rb(
     # run circuits
     results = sampler.run_batch(circuits_all, repetitions=parameters.repetitions)
     gnd_probs: dict = {q: [] for q in qubits}
+    gnd_probs_std: dict = {q: [] for q in qubits}
     idx = 0
     for num_cliffords in parameters.num_clifford_range:
         excited_probs: dict[cirq.Qid, list[float]] = {q: [] for q in qubits}
@@ -544,9 +562,17 @@ def parallel_single_qubit_rb(
             idx += 1
         for qubit in qubits:
             gnd_probs[qubit].append(1.0 - np.mean(excited_probs[qubit]))
+            gnd_probs_std[qubit].append(
+                np.std(excited_probs[qubit]) / np.sqrt(parameters.repetitions)
+            )
 
     return ParallelRandomizedBenchmarkingResult(
-        {q: RandomizedBenchMarkResult(parameters.num_clifford_range, gnd_probs[q]) for q in qubits}
+        {
+            q: RandomizedBenchMarkResult(
+                parameters.num_clifford_range, gnd_probs[q], gnd_probs_std[q]
+            )
+            for q in qubits
+        }
     )
 
 
@@ -595,6 +621,7 @@ def two_qubit_randomized_benchmarking(
     cliffords = _single_qubit_cliffords()
     cfd_matrices = _two_qubit_clifford_matrices(first_qubit, second_qubit, cliffords)
     gnd_probs = []
+    gnd_probs_std = []
     for num_cfds in num_clifford_range:
         gnd_probs_l = []
         for _ in range(num_circuits):
@@ -606,8 +633,9 @@ def two_qubit_randomized_benchmarking(
             gnds = [(not r[0] and not r[1]) for r in results.measurements['z']]
             gnd_probs_l.append(np.mean(gnds))
         gnd_probs.append(float(np.mean(gnd_probs_l)))
+        gnd_probs_std.append(float(np.std(gnd_probs_l) / np.sqrt(repetitions)))
 
-    return RandomizedBenchMarkResult(num_clifford_range, gnd_probs)
+    return RandomizedBenchMarkResult(num_clifford_range, gnd_probs, gnd_probs_std)
 
 
 def single_qubit_state_tomography(

cirq/experiments/qubit_characterizations_test.py

@@ -140,6 +140,22 @@ def test_parallel_single_qubit_parallel_single_qubit_randomized_benchmarking() -
     _ = results.plot_integrated_histogram()
 
 
+@mock.patch.dict(os.environ, clear='CIRQ_TESTING')
+def test_parallel_single_qubit_randomized_benchmarking_with_noise() -> None:
+    simulator = sim.Simulator(noise=cirq.depolarize(1e-3), seed=0)
+    qubits = (GridQubit(0, 0), GridQubit(0, 1))
+    num_cfds = range(5, 7, 1)
+    results = parallel_single_qubit_randomized_benchmarking(
+        simulator, num_clifford_range=num_cfds, repetitions=10, qubits=qubits
+    )
+    for qubit in qubits:
+        g_pops = np.asarray(results.results_dictionary[qubit].data)[:, 1]
+        assert np.isclose(np.mean(g_pops), 0.99, atol=1e-2)
+        _ = results.plot_single_qubit(qubit)
+    pauli_errors = results.pauli_error()
+    assert len(pauli_errors) == len(qubits)
+
+
 def test_two_qubit_randomized_benchmarking() -> None:
     # Check that the ground state population at the end of the Clifford
     # sequences is always unity.

cirq/protocols/unitary_protocol.py

@@ -19,7 +19,6 @@ from typing import Any, Protocol, TypeVar
 
 import numpy as np
 
-from cirq import linalg
 from cirq._doc import doc_private
 from cirq.protocols import qid_shape_protocol
 from cirq.protocols.apply_unitary_protocol import apply_unitaries, ApplyUnitaryArgs
@@ -111,11 +110,8 @@ def unitary(
     Raises:
         TypeError: `val` doesn't have a unitary effect and no default value was
             specified.
-        ValueError: `val` is a numpy array that is not unitary.
     """
     if isinstance(val, np.ndarray):
-        if not linalg.is_unitary(val):
-            raise ValueError("The provided numpy array is not unitary.")
         return val
 
     strats = [

cirq/protocols/unitary_protocol_test.py

@@ -163,12 +163,9 @@ def test_unitary():
 
     # Test that numpy arrays are handled directly
     test_matrix = np.array([[1, 0], [0, 1]])
-    assert cirq.unitary(test_matrix, NotImplemented) is test_matrix
-
-    # Test that non-unitary numpy arrays raise ValueError
+    assert cirq.unitary(test_matrix) is test_matrix
     non_unitary_matrix = np.array([[1, 1], [0, 1]])
-    with pytest.raises(ValueError, match="The provided numpy array is not unitary"):
-        _ = cirq.unitary(non_unitary_matrix)
+    assert cirq.unitary(non_unitary_matrix) is non_unitary_matrix
 
     assert cirq.unitary(NoMethod(), None) is None
     assert cirq.unitary(ReturnsNotImplemented(), None) is None

{cirq_core-1.7.0.dev20250824200454.dist-info → cirq_core-1.7.0.dev20250825201352.dist-info}/METADATA

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: cirq-core
-Version: 1.7.0.dev20250824200454
+Version: 1.7.0.dev20250825201352
 Summary: A framework for creating, editing, and invoking Noisy Intermediate Scale Quantum (NISQ) circuits.
 Home-page: http://github.com/quantumlib/cirq
 Author: The Cirq Developers

{cirq_core-1.7.0.dev20250824200454.dist-info → cirq_core-1.7.0.dev20250825201352.dist-info}/RECORD

@@ -4,8 +4,8 @@ cirq/_compat_test.py,sha256=emXpdD5ZvwLRlFAoQB8YatmZyU3b4e9jg6FppMTUhkU,33900
 cirq/_doc.py,sha256=BrnoABo1hk5RgB3Cgww4zLHUfiyFny0F1V-tOMCbdaU,2909
 cirq/_import.py,sha256=ixBu4EyGl46Ram2cP3p5eZVEFDW5L2DS-VyTjz4N9iw,8429
 cirq/_import_test.py,sha256=oF4izzOVZLc7NZ0aZHFcGv-r01eiFFt_JORx_x7_D4s,1089
-cirq/_version.py,sha256=bHMW6O934V4_oHnqf9qAPMjgAiM-Uuvt_oIxTGv0TzY,1206
-cirq/_version_test.py,sha256=K2E-4x0afwrLJNyD4mdh-YByGwtPRkI8xBoHUnSm8s8,155
+cirq/_version.py,sha256=c8XuyYustWARgawSw_lq0mSfvX3xI9eVP8Zbi1FUbn4,1206
+cirq/_version_test.py,sha256=pLSg-ptHB9RvGeE-NkfHemLecVPZiskvVzvTppAqaxY,155
 cirq/conftest.py,sha256=wSDKNdIQRDfLnXvOCWD3erheOw8JHRhdfQ53EyTUIXg,1239
 cirq/json_resolver_cache.py,sha256=A5DIgFAY1hUNt9vai_C3-gGBv24116CJMzQxMcXOax4,13726
 cirq/py.typed,sha256=VFSlmh_lNwnaXzwY-ZuW-C2Ws5PkuDoVgBdNCs0jXJE,63
@@ -199,8 +199,8 @@ cirq/experiments/n_qubit_tomography.py,sha256=16u0Tv14SyUM9WCk-ZxbBit9cl93MbZodG
 cirq/experiments/n_qubit_tomography_test.py,sha256=8wIgs0O8DtlCGOyC0MZA_d3tLNoURX1ARcqnnp1360g,4439
 cirq/experiments/purity_estimation.py,sha256=20MRvxDyxlcnWkwI94_PzMexOytS3NgnaaVbtFwLv2c,2516
 cirq/experiments/purity_estimation_test.py,sha256=ZDU6k-pja5W5kQ-5cy4id1IWUQ86WMiQYNOiXyAI0YU,967
-cirq/experiments/qubit_characterizations.py,sha256=RO-gxOI_7n4HhGOmumgLGacV4Ve8bAY_s7XyxkOTUjo,40887
-cirq/experiments/qubit_characterizations_test.py,sha256=gcMeYJLMWvboZGJqOicNAmGV736Ty1I5D5Ihd4EHOf0,11066
+cirq/experiments/qubit_characterizations.py,sha256=FaCyAiBbiuzQvVT8YNjxZnuNvyrGMlxraQTjnoPUlKA,41990
+cirq/experiments/qubit_characterizations_test.py,sha256=z5WFvZmFtvm3j2W_aF6A-dzPKvoYcW8KdD7zccW0aKk,11775
 cirq/experiments/random_quantum_circuit_generation.py,sha256=90iD2KulPXKy9mu5MiHnir2ozCJ5bgGBzxNv9odw5FY,28016
 cirq/experiments/random_quantum_circuit_generation_test.py,sha256=R4ShY12vAfLd4fTT9mX0VQM9VQol58ZdPsqBkZvFds4,16737
 cirq/experiments/readout_confusion_matrix.py,sha256=qK6qDoNRuIHGTV7PS8xGkS7tCXGGyoDfUJJj6575-cE,20664
@@ -446,8 +446,8 @@ cirq/protocols/resolve_parameters.py,sha256=N1NJJ672v3ISQRusO-t6nwm1hOmOqMVBeXtq
 cirq/protocols/resolve_parameters_test.py,sha256=2R2T2p4NkbD4IV2_4i8WkvSHu3OqjXo-Bf856Rwb-3w,4933
 cirq/protocols/trace_distance_bound.py,sha256=xF_qfkFV_T7O3-5lBITupq6ulBYuzFXDHaYDv7it96E,4150
 cirq/protocols/trace_distance_bound_test.py,sha256=0bI9uYttJj5eayM05kShPh9qkxeKG1egcZ9fXJPZWNU,1980
-cirq/protocols/unitary_protocol.py,sha256=TG-1JLsWWRhsVeI7fe0lhB8VdQx6COeWTzl2F-oS9HM,8433
-cirq/protocols/unitary_protocol_test.py,sha256=ZiU-74fCQ5CJ_KvGOBDtIu7A8K5f9lh1FZerx291Gek,10619
+cirq/protocols/unitary_protocol.py,sha256=KsXGUA7bkHOkL9GqJXpyFmhA0HVQeLdZfypZJC7yU6g,8233
+cirq/protocols/unitary_protocol_test.py,sha256=oj31ZmeGylf4wh2jv_EwsXmJ7YYgsu3YM43A-yXrPps,10480
 cirq/protocols/json_test_data/AmplitudeDampingChannel.json,sha256=x3szAuG8j_1uAK5ghFapaB410g0twQ83aQNsvItXVdo,60
 cirq/protocols/json_test_data/AmplitudeDampingChannel.repr,sha256=n_tJNGHkWlxYunXGMFtFO6-RuIv0y8Ki0YqE8w3hOl0,30
 cirq/protocols/json_test_data/AnyIntegerPowerGateFamily.json,sha256=Qf8FTwvPV7en7WcPbhP1kvjLUUPYgbICoPFqT6w86hw,68
@@ -1234,8 +1234,8 @@ cirq/work/sampler.py,sha256=rxbMWvrhu3gfNSBjZKozw28lLKVvBAS_1EGyPdYe8Xg,19041
 cirq/work/sampler_test.py,sha256=SsMrRvLDYELyOAWLKISjkdEfrBwLYWRsT6D8WrsLM3Q,13533
 cirq/work/zeros_sampler.py,sha256=Fs2JWwq0n9zv7_G5Rm-9vPeHUag7uctcMOHg0JTkZpc,2371
 cirq/work/zeros_sampler_test.py,sha256=lQLgQDGBLtfImryys2HzQ2jOSGxHgc7-koVBUhv8qYk,3345
-cirq_core-1.7.0.dev20250824200454.dist-info/licenses/LICENSE,sha256=tAkwu8-AdEyGxGoSvJ2gVmQdcicWw3j1ZZueVV74M-E,11357
-cirq_core-1.7.0.dev20250824200454.dist-info/METADATA,sha256=jdcllNVcO0GMJmNM9l9lBQrjyLpMdb2LaCTMxplz2lU,4819
-cirq_core-1.7.0.dev20250824200454.dist-info/WHEEL,sha256=_zCd3N1l69ArxyTb8rzEoP9TpbYXkqRFSNOD5OuxnTs,91
-cirq_core-1.7.0.dev20250824200454.dist-info/top_level.txt,sha256=Sz9iOxHU0IEMLSFGwiwOCaN2e9K-jFbBbtpPN1hB73g,5
-cirq_core-1.7.0.dev20250824200454.dist-info/RECORD,,
+cirq_core-1.7.0.dev20250825201352.dist-info/licenses/LICENSE,sha256=tAkwu8-AdEyGxGoSvJ2gVmQdcicWw3j1ZZueVV74M-E,11357
+cirq_core-1.7.0.dev20250825201352.dist-info/METADATA,sha256=AZ2Qp-ke_tazWFur88_hu2iaYqCHuhxPkxYZuwIrW4A,4819
+cirq_core-1.7.0.dev20250825201352.dist-info/WHEEL,sha256=_zCd3N1l69ArxyTb8rzEoP9TpbYXkqRFSNOD5OuxnTs,91
+cirq_core-1.7.0.dev20250825201352.dist-info/top_level.txt,sha256=Sz9iOxHU0IEMLSFGwiwOCaN2e9K-jFbBbtpPN1hB73g,5
+cirq_core-1.7.0.dev20250825201352.dist-info/RECORD,,