pyqrack-complex128 2.26.0__tar.gz → 2.26.2__tar.gz

{pyqrack_complex128-2.26.0/pyqrack_complex128.egg-info → pyqrack_complex128-2.26.2}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: pyqrack-complex128
-Version: 2.26.0
+Version: 2.26.2
 Summary: pyqrack - Pure Python vm6502q/qrack Wrapper
 Home-page: https://github.com/vm6502q/pyqrack
 Author: Daniel Strano

{pyqrack_complex128-2.26.0 → pyqrack_complex128-2.26.2}/pyqrack/qrack_ace_backend.py

@@ -1384,6 +1384,15 @@ class QrackAceBackend:
     def _apply_coupling(self, pauli, anti, qb1, hq1, qb2, hq2, lq1_lr, lq1=None, lq2=None):
         shadow_targets = []
 
+        if len(hq1) == len(hq2) and all(hq1[i][0] == hq2[i][0] for i in range(len(hq1))):
+            for i in range(len(hq1)):
+                sim_id, idx1 = hq1[i]
+                gate_fn, shadow_fn = self._get_gate(pauli, anti, sim_id)
+                _, idx2 = hq2[i]
+                self.sim[sim_id].swap(idx1, idx2)
+                gate_fn([idx1], idx2)
+            return
+
         # Logical-level error-detection gadget: an earlier version of this
         # also captured lq1's own logical value into a dedicated ancilla
         # (via a nested cx()), on the theory that the shadow-coupling
@@ -1410,7 +1419,7 @@ class QrackAceBackend:
         # a lq2-invariant check for that gate family, same reasoning as
         # the single-replica CZ dual-check elsewhere.
         anc2 = None
-        if self.is_error_detection and not self._in_gadget_capture and (len(hq1) > 1):
+        if self.is_error_detection and (not self._in_gadget_capture) and (len(hq2) > 1):
             anc_sim = hq1[0][0]
             t_sims = {y[0] for y in hq2 if y[0] >= 0}
             for x in hq1:
@@ -1810,42 +1819,6 @@ class QrackAceBackend:
         hq1 = self._unpack(lq1)
         hq2 = self._unpack(lq2)
 
-        sims1 = {r[0] for r in hq1}
-        sims2 = {r[0] for r in hq2}
-
-        if sims1.isdisjoint(sims2):
-            self._correct(lq1)
-            self._correct(lq2)
-
-            # Once fully resolved, the swap itself is nothing but a
-            # reference exchange: self._qubits[lq] is a reference to the
-            # replica-list object, so swapping which index points to
-            # which object IS the swap -- no separate content to copy.
-            self._qubits[lq1], self._qubits[lq2] = self._qubits[lq2], self._qubits[lq1]
-
-            lhv1 = self._lhv.pop(lq1, None)
-            lhv2 = self._lhv.pop(lq2, None)
-            if lhv2 is not None:
-                self._lhv[lq1] = lhv2
-            if lhv1 is not None:
-                self._lhv[lq2] = lhv1
-
-            # Boundary repetition code needs no bookkeeping here, unlike
-            # self._qubits/self._lhv just above: it's on-demand now, not
-            # a permanent, per-logical-index dict of pre-allocated
-            # partners -- eligibility is just len(self._qubits[lq]) > 1,
-            # computed fresh on every coupling event, so the reference
-            # swap of self._qubits above already, automatically carries
-            # the right eligibility with it. Nothing else to move.
-            return
-
-        # Boundary-to-boundary handling
-        if (len(hq1) > 1) and (len(hq2) > 1):
-            self.cx(lq1, lq2)
-            self.cx(lq2, lq1)
-            self.cx(lq1, lq2)
-            return
-
         # Fast/exact path: every replica of lq1 lines up, position-for-
         # position, with the corresponding replica of lq2 on the SAME
         # underlying simulator. A native swap() there is pure index
@@ -1872,8 +1845,6 @@ class QrackAceBackend:
         # two sides can land on different simulators (the partial-match
         # branches below) needs separate, careful treatment -- not
         # reused as-is from _cpauli's helper.
-        self._correct(lq1)
-        self._correct(lq2)
 
         if len(hq1) == len(hq2) and all(hq1[i][0] == hq2[i][0] for i in range(len(hq1))):
             for i in range(len(hq1)):
@@ -1887,6 +1858,39 @@ class QrackAceBackend:
             # with or without this call).
             return
 
+        sims1 = {r[0] for r in hq1}
+        sims2 = {r[0] for r in hq2}
+
+        if sims1.isdisjoint(sims2):
+            # Once fully resolved, the swap itself is nothing but a
+            # reference exchange: self._qubits[lq] is a reference to the
+            # replica-list object, so swapping which index points to
+            # which object IS the swap -- no separate content to copy.
+            self._qubits[lq1], self._qubits[lq2] = self._qubits[lq2], self._qubits[lq1]
+
+            lhv1 = self._lhv.pop(lq1, None)
+            lhv2 = self._lhv.pop(lq2, None)
+            if lhv2 is not None:
+                self._lhv[lq1] = lhv2
+            if lhv1 is not None:
+                self._lhv[lq2] = lhv1
+
+            # Boundary repetition code needs no bookkeeping here, unlike
+            # self._qubits/self._lhv just above: it's on-demand now, not
+            # a permanent, per-logical-index dict of pre-allocated
+            # partners -- eligibility is just len(self._qubits[lq]) > 1,
+            # computed fresh on every coupling event, so the reference
+            # swap of self._qubits above already, automatically carries
+            # the right eligibility with it. Nothing else to move.
+            return
+
+        # Boundary-to-boundary handling
+        if (len(hq1) > 1) and (len(hq2) > 1):
+            self.cx(lq1, lq2)
+            self.cx(lq2, lq1)
+            self.cx(lq1, lq2)
+            return
+
         anc1, anc2 = None, None
         if self.is_error_detection and not self._in_gadget_capture and ((len(hq1) > 1) or (len(hq2) > 1)):
             t_sims = {y[0] for y in hq2 if y[0] >= 0}
@@ -1902,6 +1906,9 @@ class QrackAceBackend:
             self.cx(lq2, anc2)
             self._in_gadget_capture = False
 
+        self._correct(lq1)
+        self._correct(lq2)
+
         # Partial-match cases: one side is a simple (single-replica) qubit,
         # the other has multiple replicas. For the multi-replica side's
         # non-matching replicas, the correct decomposition (verified
@@ -2001,10 +2008,6 @@ class QrackAceBackend:
             self.mcx([c1, c2], t)
             return
 
-        if not self._in_gadget_capture:
-            self._correct(c1)
-            self._correct(c2)
-
         # Control-protection ancillas (anc1/anc1b/anc2/anc2b, checking
         # whether c1/c2's own values held still) removed entirely, same
         # basis as _apply_coupling: instrumented and stress-tested
@@ -2028,7 +2031,7 @@ class QrackAceBackend:
         # bulk on this path (the len(hqt) > 1 case returns earlier,
         # above) and nothing else in this method uses that array anymore.
         anc_and = None
-        if self.is_error_detection and not self._in_gadget_capture and ((len(hq1) > 1) or (len(hq2) > 1) or (len(hqt) > 1)):
+        if self.is_error_detection and not self._in_gadget_capture and (len(hqt) > 1):
             found = False
             anc_sim = hq1[0][0]
             t_sims = {y[0] for y in hqt if y[0] >= 0}
@@ -2070,6 +2073,9 @@ class QrackAceBackend:
                 self.ccx(c1, c2, anc_and)
             self._in_gadget_capture = False
 
+        self._correct(c1)
+        self._correct(c2)
+
         qb1, _ = QrackAceBackend._get_qb_lhv_indices(hq1)
         qb2, _ = QrackAceBackend._get_qb_lhv_indices(hq2)
         qbt, _ = QrackAceBackend._get_qb_lhv_indices(hqt)
@@ -2136,6 +2142,8 @@ class QrackAceBackend:
             for (b1, b2, bt) in witnesses:
                 self.sim[b1[0]].mcx([b1[1], b2[1]], bt[1])
 
+        self._correct(t)
+
         self._in_gadget_capture = True
 
         if anc_and is not None:
@@ -2153,9 +2161,6 @@ class QrackAceBackend:
 
         self._in_gadget_capture = False
 
-        if not self._in_gadget_capture:
-            self._correct(t)
-
     def ccz(self, c1, c2, t):
         self.h(t)
         self.ccx(c1, c2, t)

{pyqrack_complex128-2.26.0 → pyqrack_complex128-2.26.2/pyqrack_complex128.egg-info}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: pyqrack-complex128
-Version: 2.26.0
+Version: 2.26.2
 Summary: pyqrack - Pure Python vm6502q/qrack Wrapper
 Home-page: https://github.com/vm6502q/pyqrack
 Author: Daniel Strano

{pyqrack_complex128-2.26.0 → pyqrack_complex128-2.26.2}/setup.py

@@ -7,7 +7,7 @@ from setuptools import setup
 from setuptools.command.build_py import build_py
 
 
-VERSION = "2.26.0"
+VERSION = "2.26.2"
 
 # Read long description from README.
 README_PATH = os.path.join(os.path.abspath(os.path.dirname(__file__)), 'README.md')