pyqrack-complex128 2.26.1__tar.gz → 2.26.4__tar.gz

{pyqrack_complex128-2.26.1/pyqrack_complex128.egg-info → pyqrack_complex128-2.26.4}/PKG-INFO

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

{pyqrack_complex128-2.26.1 → pyqrack_complex128-2.26.4}/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:
@@ -1533,7 +1542,13 @@ class QrackAceBackend:
         anc1 = None
         if self.is_boundary_repetition_code and (not self._in_gadget_capture) and lq1_lr and (not lq2_lr):
             # Encode:
-            anc1 = self._rep_code_ancilla_lq[hq1[0][0]]
+            anc_sim = hq1[0][0]
+            t_sims = {y[0] for y in hq2 if y[0] >= 0}
+            for x in hq1:
+                if x[0] in t_sims:
+                    anc_sim = x[0]
+                    break
+            anc1 = self._detect_ancilla2_lq[anc_sim]
             self._in_gadget_capture = True
             self.cx(lq2, anc1)
             self._in_gadget_capture = False
@@ -1810,9 +1825,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}
-
         # 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
@@ -1852,6 +1864,9 @@ 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
@@ -2022,7 +2037,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}
@@ -2864,7 +2879,7 @@ class QrackAceBackend:
         return tgt
 
     # Designed by Dan, Elara (ChatGPT), and (Anthropic) Claude:
-    def create_noise_model(self, x=0.5, y=0.5):
+    def create_noise_model(self, x=0.25, y=0.25):
         if not _IS_QISKIT_AER_AVAILABLE:
             raise RuntimeError(
                 "Before trying to run_qiskit_circuit() with QrackAceBackend, you must install Qiskit Aer!"
@@ -2872,17 +2887,10 @@ class QrackAceBackend:
         noise_model = NoiseModel()
 
         def _uncommon_sim_fraction(lq1, lq2):
-            sims1 = [qb[0] for qb in self._qubits[lq1]]
-            sims2 = [qb[0] for qb in self._qubits[lq2]]
-            n = 0
-            for s in sims1:
-                if s not in sims2:
-                    n += 1
-            for s in sims2:
-                if s not in sims1:
-                    n += 1
-
-            return n / (len(sims1) + len(sims2))
+            sims1 = {qb[0] for qb in self._qubits[lq1]}
+            sims2 = {qb[0] for qb in self._qubits[lq2]}
+
+            return len(sims1 ^ sims2) / len(sims1.union(sims2))
 
         coupling_map = self.get_logical_coupling_map()
 
@@ -2938,6 +2946,7 @@ class QrackAceBackend:
 
             sp = (1 - self._sdrp / 2) ** c
             p = (1 - x) ** u
+            d = 1 - ((1 - y) ** u)
 
             # cx/cy/cz: single-qubit, TARGET-only error, with the gate's
             # own generator as the Pauli type -- derived directly from
@@ -2962,7 +2971,7 @@ class QrackAceBackend:
                 # anc2 verifying b) are structurally live exactly when a
                 # itself is boundary -- see the comment above this loop.
                 # Damp accordingly.
-                p_shadow_ccz *= y
+                p_shadow_ccz *= d
             noise_model.add_quantum_error(
                 pauli_error([("ZI", p_shadow_ccz), ("II", 1 - p_shadow_ccz)]), "cz", [a, b]
             )
@@ -3048,7 +3057,7 @@ class QrackAceBackend:
             if (is_a_simple != is_b_simple) and has_match:
                 p_net_swap = 2 * p * (1 - p)
                 if self.is_error_detection:
-                    p_net_swap *= y
+                    p_net_swap *= d
                 if is_a_simple:
                     # a is bulk, b is boundary -> error lands on b
                     noise_model.add_quantum_error(
@@ -3099,7 +3108,7 @@ class QrackAceBackend:
                     # from a different gate (the now-protected native
                     # swap, not
                     # this cz component).
-                    p_cz *= y
+                    p_cz *= d
                 p_net_iswap = p_net_swap * (1 - p_cz) + p_cz * (1 - p_net_swap)
                 if is_a_simple:
                     terms = [("II", 1 - p_net_iswap), ("ZI", p_net_iswap)]

{pyqrack_complex128-2.26.1 → pyqrack_complex128-2.26.4/pyqrack_complex128.egg-info}/PKG-INFO

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

{pyqrack_complex128-2.26.1 → pyqrack_complex128-2.26.4}/pyqrack_complex128.egg-info/SOURCES.txt

@@ -17,9 +17,6 @@ pyqrack/qrack_stabilizer.py
 pyqrack/quimb_circuit_type.py
 pyqrack/qrack_system/__init__.py
 pyqrack/qrack_system/qrack_system.py
-pyqrack/qrack_system/qrack_cl_precompile/qrack_cl_precompile
-pyqrack/qrack_system/qrack_lib/libqrack_pinvoke.so
-pyqrack/qrack_system/qrack_lib/libqrack_pinvoke.so.10.13.1
 pyqrack/stats/__init__.py
 pyqrack/stats/load_quantized_data.py
 pyqrack/stats/quantize_by_range.py

{pyqrack_complex128-2.26.1 → pyqrack_complex128-2.26.4}/setup.py

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