pyqrack-cuda 1.44.13__tar.gz → 1.44.15__tar.gz

{pyqrack_cuda-1.44.13/pyqrack_cuda.egg-info → pyqrack_cuda-1.44.15}/PKG-INFO

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

{pyqrack_cuda-1.44.13 → pyqrack_cuda-1.44.15}/pyqrack/qrack_ace_backend.py

@@ -120,6 +120,9 @@ class QrackAceBackend:
             self.sim.mcx([hq[0]], hq[1])
 
     def _correct(self, lq):
+        # We can't use true syndrome-based error correction,
+        # because one of the qubits in the code is separated.
+        # However, we can get pretty close!
         hq = self._unpack(lq)
         shots = 1024
         samples = self.sim.measure_shots(hq, shots)
@@ -138,92 +141,111 @@ class QrackAceBackend:
                     syndrome[1] += 1
                 case 3:
                     syndrome[2] += 1
+
+        row = (hq[0] // 3) // self.row_length
+        even_row = not (row & 1)
+        single_bit = 0
+        other_bits = []
+        if (not self.alternating_codes or even_row):
+            single_bit = 2
+            other_bits = [0, 1]
+        else:
+            single_bit = 0
+            other_bits = [1, 2]
+
         max_syndrome = max(syndrome)
-        if (2 * max_syndrome) > shots:
-            error_bit = syndrome.index(max_syndrome)
-            row = (hq[0] // 3) // self.row_length
-            even_row = not (row & 1)
-            single_bit = 2 if (not self.alternating_codes or even_row) else 0
+        error_bit = syndrome.index(max_syndrome)
+        force_syndrome = True
+        if (2 * max_syndrome) >= shots:
+            # There is an error.
             if error_bit == single_bit:
+                # The stand-alone bit carries the error.
                 self.sim.x(hq[error_bit])
             else:
-                self.sim.mcx([hq[1]], self._ancilla)
-                self.sim.mcx([hq[2]], self._ancilla)
+                # The coherent bits carry the error.
+                force_syndrome = False
+                # Form their syndrome.
+                self.sim.mcx([hq[other_bits[0]]], self._ancilla)
+                self.sim.mcx([hq[other_bits[1]]], self._ancilla)
+                # Force the syndrome pathological
                 self.sim.force_m(self._ancilla, True)
+                # Reset the ancilla.
+                self.sim.x(self._ancilla)
+                # Correct the bit flip.
                 self.sim.x(hq[error_bit])
 
+        # There is no error.
+        if force_syndrome:
+            # Form the syndrome of the coherent bits.
+            self.sim.mcx([hq[other_bits[0]]], self._ancilla)
+            self.sim.mcx([hq[other_bits[1]]], self._ancilla)
+            # Force the syndrome non-pathological.
+            self.sim.force_m(self._ancilla, False)
+
 
     def u(self, th, ph, lm, lq):
         hq = self._unpack(lq)
         self._decode(hq)
         self.sim.u(hq[0], th, ph, lm)
         self._encode(hq)
-        self._correct(lq)
 
     def r(self, p, th, lq):
         hq = self._unpack(lq)
         self._decode(hq)
         self.sim.r(p, th, hq[0])
         self._encode(hq)
-        self._correct(lq)
 
     def s(self, lq):
         hq = self._unpack(lq)
         self._decode(hq)
         self.sim.s(hq[0])
         self._encode(hq)
-        self._correct(lq)
 
     def adjs(self, lq):
         hq = self._unpack(lq)
         self._decode(hq)
         self.sim.adjs(hq[0])
         self._encode(hq)
-        self._correct(lq)
 
     def x(self, lq):
         hq = self._unpack(lq)
         self._decode(hq)
         self.sim.x(hq[0])
         self._encode(hq)
-        self._correct(lq)
 
     def y(self, lq):
         hq = self._unpack(lq)
         self._decode(hq)
         self.sim.y(hq[0])
         self._encode(hq)
-        self._correct(lq)
 
     def z(self, lq):
         hq = self._unpack(lq)
         self._decode(hq)
         self.sim.z(hq[0])
         self._encode(hq)
-        self._correct(lq)
 
     def h(self, lq):
         hq = self._unpack(lq)
         self._decode(hq)
         self.sim.h(hq[0])
         self._encode(hq)
-        self._correct(lq)
 
     def t(self, lq):
         hq = self._unpack(lq)
         self._decode(hq)
         self.sim.t(hq[0])
         self._encode(hq)
-        self._correct(lq)
 
     def adjt(self, lq):
         hq = self._unpack(lq)
         self._decode(hq)
         self.sim.adjt(hq[0])
         self._encode(hq)
-        self._correct(lq)
 
     def _cpauli(self, lq1, lq2, anti, pauli):
+        self._correct(lq1)
+
         gate = None
         shadow = None
         if pauli == Pauli.PauliX:
@@ -271,9 +293,6 @@ class QrackAceBackend:
                 gate([hq1[1]], hq2[1])
             gate([hq1[2]], hq2[2])
 
-        self._correct(lq1)
-        self._correct(lq2)
-
 
     def cx(self, lq1, lq2):
         self._cpauli(lq1, lq2, False, Pauli.PauliX)

{pyqrack_cuda-1.44.13 → pyqrack_cuda-1.44.15/pyqrack_cuda.egg-info}/PKG-INFO

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

{pyqrack_cuda-1.44.13 → pyqrack_cuda-1.44.15}/setup.py

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