cirq-core 1.4.0.dev20240430181951__py3-none-any.whl → 1.4.0.dev20240506164407__py3-none-any.whl

cirq/_version.py

@@ -1 +1 @@
-__version__ = "1.4.0.dev20240430181951"
+__version__ = "1.4.0.dev20240506164407"

cirq/ops/clifford_gate.py

@@ -400,7 +400,11 @@ class CliffordGate(raw_types.Gate, CommonCliffordGates):
         # By definition, Clifford Gate should always return True.
         return True
 
-    def __pow__(self, exponent) -> 'CliffordGate':
+    def __pow__(self, exponent: float) -> 'CliffordGate':
+        if exponent != int(exponent):
+            return NotImplemented
+        exponent = int(exponent)
+
         if exponent == -1:
             return CliffordGate.from_clifford_tableau(self.clifford_tableau.inverse())
         if exponent == 0:
@@ -409,18 +413,24 @@ class CliffordGate(raw_types.Gate, CommonCliffordGates):
             )
         if exponent == 1:
             return self
-        if exponent > 0 and int(exponent) == exponent:
-            base_tableau = self.clifford_tableau.copy()
-            for _ in range(int(exponent) - 1):
-                base_tableau = base_tableau.then(self.clifford_tableau)
-            return CliffordGate.from_clifford_tableau(base_tableau)
-        if exponent < 0 and int(exponent) == exponent:
-            base_tableau = self.clifford_tableau.copy()
-            for _ in range(int(-exponent) - 1):
-                base_tableau = base_tableau.then(self.clifford_tableau)
-            return CliffordGate.from_clifford_tableau(base_tableau.inverse())
 
-        return NotImplemented
+        base_tableau = self.clifford_tableau.copy()
+        if exponent < 0:
+            base_tableau = base_tableau.inverse()
+            exponent = abs(exponent)
+
+        # https://cp-algorithms.com/algebra/binary-exp.html
+        aux = qis.CliffordTableau(
+            num_qubits=self.clifford_tableau.n
+        )  # this tableau collects the odd terms
+        while exponent > 1:
+            if exponent & 1:
+                aux = aux.then(base_tableau)
+            base_tableau = base_tableau.then(base_tableau)
+            exponent >>= 1
+
+        base_tableau = base_tableau.then(aux)
+        return CliffordGate.from_clifford_tableau(base_tableau)
 
     def __repr__(self) -> str:
         return f"Clifford Gate with Tableau:\n {self.clifford_tableau._str_full_()}"

cirq/qis/clifford_tableau.py

@@ -129,7 +129,9 @@ class StabilizerState(
 
 class CliffordTableau(StabilizerState):
     """Tableau representation of a stabilizer state
-    (based on Aaronson and Gottesman 2006).
+
+    References:
+        - [Aaronson and Gottesman](https://arxiv.org/abs/quant-ph/0406196)
 
     The tableau stores the stabilizer generators of
     the state using three binary arrays: xs, zs, and rs.

cirq/transformers/gauge_compiling/gauge_compiling.py

@@ -71,6 +71,7 @@ class ConstantGauge(Gauge):
     post_q1: Tuple[ops.Gate, ...] = field(
         default=(), converter=lambda g: (g,) if isinstance(g, ops.Gate) else tuple(g)
     )
+    swap_qubits: bool = False
 
     def sample(self, gate: ops.Gate, prng: np.random.Generator) -> "ConstantGauge":
         return self
@@ -85,6 +86,41 @@ class ConstantGauge(Gauge):
         """A tuple (ops to apply to q0, ops to apply to q1)."""
         return self.post_q0, self.post_q1
 
+    def on(self, q0: ops.Qid, q1: ops.Qid) -> ops.Operation:
+        """Returns the operation that replaces the two qubit gate."""
+        if self.swap_qubits:
+            return self.two_qubit_gate(q1, q0)
+        return self.two_qubit_gate(q0, q1)
+
+
+@frozen
+class SameGateGauge(Gauge):
+    """Same as ConstantGauge but the new two-qubit gate equals the old gate."""
+
+    pre_q0: Tuple[ops.Gate, ...] = field(
+        default=(), converter=lambda g: (g,) if isinstance(g, ops.Gate) else tuple(g)
+    )
+    pre_q1: Tuple[ops.Gate, ...] = field(
+        default=(), converter=lambda g: (g,) if isinstance(g, ops.Gate) else tuple(g)
+    )
+    post_q0: Tuple[ops.Gate, ...] = field(
+        default=(), converter=lambda g: (g,) if isinstance(g, ops.Gate) else tuple(g)
+    )
+    post_q1: Tuple[ops.Gate, ...] = field(
+        default=(), converter=lambda g: (g,) if isinstance(g, ops.Gate) else tuple(g)
+    )
+    swap_qubits: bool = False
+
+    def sample(self, gate: ops.Gate, prng: np.random.Generator) -> ConstantGauge:
+        return ConstantGauge(
+            two_qubit_gate=gate,
+            pre_q0=self.pre_q0,
+            pre_q1=self.pre_q1,
+            post_q0=self.post_q0,
+            post_q1=self.post_q1,
+            swap_qubits=self.swap_qubits,
+        )
+
 
 def _select(choices: Sequence[Gauge], probabilites: np.ndarray, prng: np.random.Generator) -> Gauge:
     return choices[prng.choice(len(choices), p=probabilites)]
@@ -154,7 +190,7 @@ class GaugeTransformer:
                     gauge = self.gauge_selector(rng).sample(op.gate, rng)
                     q0, q1 = op.qubits
                     left.extend([g(q) for g in gs] for q, gs in zip(op.qubits, gauge.pre))
-                    center.append(gauge.two_qubit_gate(q0, q1))
+                    center.append(gauge.on(q0, q1))
                     right.extend([g(q) for g in gs] for q, gs in zip(op.qubits, gauge.post))
                 else:
                     center.append(op)

cirq/transformers/gauge_compiling/spin_inversion_gauge.py

@@ -17,17 +17,17 @@
 from cirq.transformers.gauge_compiling.gauge_compiling import (
     GaugeTransformer,
     GaugeSelector,
-    ConstantGauge,
+    SameGateGauge,
 )
 from cirq import ops
 
 SpinInversionGaugeSelector = GaugeSelector(
     gauges=[
-        ConstantGauge(two_qubit_gate=ops.ZZ, pre_q0=ops.X, post_q0=ops.X),
-        ConstantGauge(two_qubit_gate=ops.ZZ, pre_q1=ops.X, post_q1=ops.X),
+        SameGateGauge(pre_q0=ops.X, post_q0=ops.X, pre_q1=ops.X, post_q1=ops.X),
+        SameGateGauge(),
     ]
 )
 
 SpinInversionGaugeTransformer = GaugeTransformer(
-    target=ops.ZZ, gauge_selector=SpinInversionGaugeSelector
+    target=ops.GateFamily(ops.ZZPowGate), gauge_selector=SpinInversionGaugeSelector
 )

cirq/transformers/gauge_compiling/spin_inversion_gauge_test.py

@@ -12,12 +12,26 @@
 # See the License for the specific language governing permissions and
 # limitations under the License.
 
-
 import cirq
 from cirq.transformers.gauge_compiling import SpinInversionGaugeTransformer
 from cirq.transformers.gauge_compiling.gauge_compiling_test_utils import GaugeTester
 
 
-class TestSpinInversionGauge(GaugeTester):
+class TestSpinInversionGauge_0(GaugeTester):
     two_qubit_gate = cirq.ZZ
     gauge_transformer = SpinInversionGaugeTransformer
+
+
+class TestSpinInversionGauge_1(GaugeTester):
+    two_qubit_gate = cirq.ZZ**0.1
+    gauge_transformer = SpinInversionGaugeTransformer
+
+
+class TestSpinInversionGauge_2(GaugeTester):
+    two_qubit_gate = cirq.ZZ**-1
+    gauge_transformer = SpinInversionGaugeTransformer
+
+
+class TestSpinInversionGauge_3(GaugeTester):
+    two_qubit_gate = cirq.ZZ**0.3
+    gauge_transformer = SpinInversionGaugeTransformer

cirq/transformers/gauge_compiling/sqrt_cz_gauge.py

@@ -12,18 +12,53 @@
 # See the License for the specific language governing permissions and
 # limitations under the License.
 
-"""A Gauge transformer for CZ**0.5 gate."""
+"""A Gauge transformer for CZ**0.5 and CZ**-0.5 gates."""
+
+
+from typing import TYPE_CHECKING
+import numpy as np
 
 from cirq.transformers.gauge_compiling.gauge_compiling import (
     GaugeTransformer,
     GaugeSelector,
     ConstantGauge,
+    Gauge,
 )
-from cirq.ops.common_gates import CZ
-from cirq import ops
+from cirq.ops import CZ, S, X, Gateset
+
+if TYPE_CHECKING:
+    import cirq
+
+_SQRT_CZ = CZ**0.5
+_ADJ_S = S**-1
 
-SqrtCZGaugeSelector = GaugeSelector(
-    gauges=[ConstantGauge(pre_q0=ops.X, post_q0=ops.X, post_q1=ops.Z**0.5, two_qubit_gate=CZ**-0.5)]
-)
 
-SqrtCZGaugeTransformer = GaugeTransformer(target=CZ**0.5, gauge_selector=SqrtCZGaugeSelector)
+class SqrtCZGauge(Gauge):
+
+    def weight(self) -> float:
+        return 3.0
+
+    def sample(self, gate: 'cirq.Gate', prng: np.random.Generator) -> ConstantGauge:
+        if prng.choice([True, False]):
+            return ConstantGauge(two_qubit_gate=gate)
+        swap_qubits = prng.choice([True, False])
+        if swap_qubits:
+            return ConstantGauge(
+                pre_q1=X,
+                post_q1=X,
+                post_q0=S if gate == _SQRT_CZ else _ADJ_S,
+                two_qubit_gate=gate**-1,
+                swap_qubits=True,
+            )
+        else:
+            return ConstantGauge(
+                pre_q0=X,
+                post_q0=X,
+                post_q1=S if gate == _SQRT_CZ else _ADJ_S,
+                two_qubit_gate=gate**-1,
+            )
+
+
+SqrtCZGaugeTransformer = GaugeTransformer(
+    target=Gateset(_SQRT_CZ, _SQRT_CZ**-1), gauge_selector=GaugeSelector(gauges=[SqrtCZGauge()])
+)

cirq/transformers/gauge_compiling/sqrt_cz_gauge_test.py

@@ -12,7 +12,6 @@
 # See the License for the specific language governing permissions and
 # limitations under the License.
 
-
 import cirq
 from cirq.transformers.gauge_compiling import SqrtCZGaugeTransformer
 from cirq.transformers.gauge_compiling.gauge_compiling_test_utils import GaugeTester
@@ -21,3 +20,8 @@ from cirq.transformers.gauge_compiling.gauge_compiling_test_utils import GaugeTe
 class TestSqrtCZGauge(GaugeTester):
     two_qubit_gate = cirq.CZ**0.5
     gauge_transformer = SqrtCZGaugeTransformer
+
+
+class TestAdjointSqrtCZGauge(GaugeTester):
+    two_qubit_gate = cirq.CZ**-0.5
+    gauge_transformer = SqrtCZGaugeTransformer

{cirq_core-1.4.0.dev20240430181951.dist-info → cirq_core-1.4.0.dev20240506164407.dist-info}/METADATA

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: cirq-core
-Version: 1.4.0.dev20240430181951
+Version: 1.4.0.dev20240506164407
 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.4.0.dev20240430181951.dist-info → cirq_core-1.4.0.dev20240506164407.dist-info}/RECORD

@@ -4,7 +4,7 @@ cirq/_compat_test.py,sha256=Qq3ZcfgD-Nb81cEppQdJqhAyrVqXKtfXZYGXT0p-Wh0,34718
 cirq/_doc.py,sha256=yDyWUD_2JDS0gShfGRb-rdqRt9-WeL7DhkqX7np0Nko,2879
 cirq/_import.py,sha256=p9gMHJscbtDDkfHOaulvd3Aer0pwUF5AXpL89XR8dNw,8402
 cirq/_import_test.py,sha256=6K_v0riZJXOXUphHNkGA8MY-JcmGlezFaGmvrNhm3OQ,1015
-cirq/_version.py,sha256=sl6JoWpxiZDtf2vZLNMuq-vBb6mXIBq5SXaP3tQKItw,40
+cirq/_version.py,sha256=ssE6yg0quJRx6t90raQtx0xyl-ZWXvP37SExL_xDZBA,40
 cirq/_version_test.py,sha256=yYS6xm5-nuBRQJa9R3n41WdvFtVyY7Lb5Q8bea3VgBI,133
 cirq/conftest.py,sha256=X7yLFL8GLhg2CjPw0hp5e_dGASfvHx1-QT03aUbhKJw,1168
 cirq/json_resolver_cache.py,sha256=9g_JQMmfBzTuV-3s2flUbXIgcLjs4K7LjAFFgngdG1U,13204
@@ -269,7 +269,7 @@ cirq/ops/boolean_hamiltonian.py,sha256=hjmzBvWiZnb18JfdhId-dVxt2VmbyAbqMSFAPWV7s
 cirq/ops/boolean_hamiltonian_test.py,sha256=1ey5yfYZPKZDsfM3jpCPAOpbPs_y8i4K_WvDK2d5_4Y,8518
 cirq/ops/classically_controlled_operation.py,sha256=bMqKutwzqbvN2r7mOVOI12HTPDGSJLkhQQgfAcTtbDU,9166
 cirq/ops/classically_controlled_operation_test.py,sha256=MNU0Adff4kTosqsrQ3PUT3ARcZee_PkchT6u0xDl8Qg,48039
-cirq/ops/clifford_gate.py,sha256=xTAGYR24u__ctlY2fNsWYTQcJof4euRYMGN8cdERV2E,39235
+cirq/ops/clifford_gate.py,sha256=ARbHGptZxQo6kX9QTrh_kZ89N7lW6V3GQSE3o-1kixU,39366
 cirq/ops/clifford_gate_test.py,sha256=NF_if1X8LCMA9hy0vBO7lxvVPdumlvMMnI2XRQ-RLpk,37472
 cirq/ops/common_channels.py,sha256=mVGEta2wnaWOhRCXQO6gUj9Zll3rtXqpY-PRYcJ8m1U,38282
 cirq/ops/common_channels_test.py,sha256=EL_PxbqD3KPC8ioihiukhmW8bUdclqqigKoFyUQpmIM,29690
@@ -873,7 +873,7 @@ cirq/protocols/json_test_data/sympy.pi.repr,sha256=ZQS0my0esr3dWTZ3mWlqgR63uorPC
 cirq/qis/__init__.py,sha256=C_gnZt9CNWEE0I1IAZE0liNVIv-aphUohIHJih-szDg,1836
 cirq/qis/channels.py,sha256=AxKgLUbWLrb1pz9xLtSpYm_stjN-IWOwLFYC9YZSons,12798
 cirq/qis/channels_test.py,sha256=iqkSyfvx_c2ZpJKAVEsFyZRmKpzNrJsyFjHbaJUMqcQ,14454
-cirq/qis/clifford_tableau.py,sha256=Ik1ogTElaB3nJCPCGj1yqxMdpu404AbZjkVOVYZOins,25305
+cirq/qis/clifford_tableau.py,sha256=X25eCXWEUWzoB-x6cna4LdLMFVNxUiIEUxKHJkkc2aE,25353
 cirq/qis/clifford_tableau_test.py,sha256=Z-5YSysJw2H3AfgZxhptvUMr5GUjUFgPcAUdcust9Qk,18441
 cirq/qis/measures.py,sha256=-e2mStl9882z4mbLwWtbdmmdJilmf6Ex2b_jSXhMX3c,12194
 cirq/qis/measures_test.py,sha256=3LTTGdvuFfZWmFyWFeUHpg1yGc9z0VE8j7Kd7J7y8i4,10092
@@ -1074,16 +1074,16 @@ cirq/transformers/analytical_decompositions/two_qubit_to_sqrt_iswap_test.py,sha2
 cirq/transformers/gauge_compiling/__init__.py,sha256=dmgVx7_yfsO9TarOCMBfVuZjza_6_MmUyAfM_eo698E,1146
 cirq/transformers/gauge_compiling/cz_gauge.py,sha256=TNZviXFu4j-lCF87QMGYVdb2RC_ePHLdI6FRCqh9go4,2550
 cirq/transformers/gauge_compiling/cz_gauge_test.py,sha256=_2RMzwuMoqytgsVZEET2m6QsGoZ2_uWgBfSGGvhZEWo,854
-cirq/transformers/gauge_compiling/gauge_compiling.py,sha256=GCYpaZIUxkcz16s8DI64dtpHhSGiIy6JIplvM6GUklc,6662
+cirq/transformers/gauge_compiling/gauge_compiling.py,sha256=7jUg2MHVX2V8SzH7r27e_BUUA_p5l15-Qv8lLlhGasQ,7948
 cirq/transformers/gauge_compiling/gauge_compiling_test.py,sha256=xVeoaaMuenlp_pxqmZO89kJEeBV0iapK7-JDTUuC7Lg,1565
 cirq/transformers/gauge_compiling/gauge_compiling_test_utils.py,sha256=6mpF2ftcDtWQEsZNr-m16GkDOn5WBWV08JB4ZS5rf8M,3408
 cirq/transformers/gauge_compiling/gauge_compiling_test_utils_test.py,sha256=HQw1xeko9OAD0PCuKcehUiy8mbxeRza4nE-6E8MGVME,1584
 cirq/transformers/gauge_compiling/iswap_gauge.py,sha256=UGJ_061h65Rfgb9LWREjxC8OSt01ZqP9TGnacL8VAuk,3500
 cirq/transformers/gauge_compiling/iswap_gauge_test.py,sha256=HEIIwKlX5ixau1e_etSUj5NvYOVTT-Gc3kuHcyKAeJ4,866
-cirq/transformers/gauge_compiling/spin_inversion_gauge.py,sha256=BMt1g2h23J9jb6LTlZsNczMKIY8J0duHhZNu0jRB48M,1107
-cirq/transformers/gauge_compiling/spin_inversion_gauge_test.py,sha256=PtbHNvh6Oz-n6Tn0fdvpgHGgaPuJQPlJ8UfHd59-dks,887
-cirq/transformers/gauge_compiling/sqrt_cz_gauge.py,sha256=l0Dkz2a7_ljKzzP-Rr9f_yY4pqJ8TGY_R-FVEkBb77I,1050
-cirq/transformers/gauge_compiling/sqrt_cz_gauge_test.py,sha256=wB_WezbTBbkwnTs7hz_tXmU4syuMulGirDDDjzqvNt0,871
+cirq/transformers/gauge_compiling/spin_inversion_gauge.py,sha256=gfjSlQdo13GfBPlrkQoHPWWzouiV7yYr7JAaB85NSGY,1086
+cirq/transformers/gauge_compiling/spin_inversion_gauge_test.py,sha256=FJwl11VPNXkf_v0WjCg9VxNKAplZ_5bZR7guOfnK-aI,1292
+cirq/transformers/gauge_compiling/sqrt_cz_gauge.py,sha256=32OGTcYT3tBFEQ1GQlrssc1wtwCcSvk4ZC0I1XD1QXg,1869
+cirq/transformers/gauge_compiling/sqrt_cz_gauge_test.py,sha256=RwjadOOJfa2Qf7iryTIMJLPzeDMNqFkP6Tewjq68gJI,997
 cirq/transformers/gauge_compiling/sqrt_iswap_gauge.py,sha256=dqQa-UWq31bE_jF3KMhU76sND5GuqTpEp9-wVuXdWVM,3126
 cirq/transformers/gauge_compiling/sqrt_iswap_gauge_test.py,sha256=0CLZoLw-WK3aKEIoaKBrQZ-qvaprOVLad-dVyWFmSiU,882
 cirq/transformers/heuristic_decompositions/__init__.py,sha256=2KBAUk-vsVku-HuahLbDQxp8bjEQMZUUNiBz_cvUMUw,831
@@ -1169,8 +1169,8 @@ cirq/work/sampler.py,sha256=JEAeQQRF3bqlO9AkOf4XbrTATDI5f5JgyM_FAUCNxao,19751
 cirq/work/sampler_test.py,sha256=B2ZsuqGT854gQtBIAh8k0LiG9Vj5wSzcGvkxOUoTcW4,13217
 cirq/work/zeros_sampler.py,sha256=x1C7cup66a43n-3tm8QjhiqJa07qcJW10FxNp9jJ59Q,2356
 cirq/work/zeros_sampler_test.py,sha256=JIkpBBFPJe5Ba4142vzogyWyboG1Q1ZAm0UVGgOoZn8,3279
-cirq_core-1.4.0.dev20240430181951.dist-info/LICENSE,sha256=tAkwu8-AdEyGxGoSvJ2gVmQdcicWw3j1ZZueVV74M-E,11357
-cirq_core-1.4.0.dev20240430181951.dist-info/METADATA,sha256=ridrMaAz813yPBhStGmnK4X6S_TW3CoO7VnVs8HH4nI,2098
-cirq_core-1.4.0.dev20240430181951.dist-info/WHEEL,sha256=GJ7t_kWBFywbagK5eo9IoUwLW6oyOeTKmQ-9iHFVNxQ,92
-cirq_core-1.4.0.dev20240430181951.dist-info/top_level.txt,sha256=Sz9iOxHU0IEMLSFGwiwOCaN2e9K-jFbBbtpPN1hB73g,5
-cirq_core-1.4.0.dev20240430181951.dist-info/RECORD,,
+cirq_core-1.4.0.dev20240506164407.dist-info/LICENSE,sha256=tAkwu8-AdEyGxGoSvJ2gVmQdcicWw3j1ZZueVV74M-E,11357
+cirq_core-1.4.0.dev20240506164407.dist-info/METADATA,sha256=6ilalyHUj2Tq2g4Fwy6ptYf4dX8Hq7prrwYw5AYJT28,2098
+cirq_core-1.4.0.dev20240506164407.dist-info/WHEEL,sha256=GJ7t_kWBFywbagK5eo9IoUwLW6oyOeTKmQ-9iHFVNxQ,92
+cirq_core-1.4.0.dev20240506164407.dist-info/top_level.txt,sha256=Sz9iOxHU0IEMLSFGwiwOCaN2e9K-jFbBbtpPN1hB73g,5
+cirq_core-1.4.0.dev20240506164407.dist-info/RECORD,,