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

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

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

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

@@ -39,7 +39,8 @@ class QrackAceBackend:
     """
 
     def __init__(self, qubit_count=-1, alternating_codes=True, toClone=None):
-        self.sim = toClone.sim.clone() if toClone else QrackSimulator(3 * qubit_count)
+        self.sim = toClone.sim.clone() if toClone else QrackSimulator(3 * qubit_count + 1)
+        self._ancilla = 3 * qubit_count
         self._factor_width(qubit_count)
         self.alternating_codes = alternating_codes
 
@@ -118,65 +119,109 @@ class QrackAceBackend:
             self._cx_shadow(hq[1], hq[2])
             self.sim.mcx([hq[0]], hq[1])
 
+    def _correct(self, lq):
+        hq = self._unpack(lq)
+        shots = 1024
+        samples = self.sim.measure_shots(hq, shots)
+        syndrome = [0, 0, 0]
+        for sample in samples:
+            match sample:
+                case 1:
+                    syndrome[0] += 1
+                case 2:
+                    syndrome[1] += 1
+                case 4:
+                    syndrome[2] += 1
+                case 6:
+                    syndrome[0] += 1
+                case 5:
+                    syndrome[1] += 1
+                case 3:
+                    syndrome[2] += 1
+        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
+            if error_bit == single_bit:
+                self.sim.x(hq[error_bit])
+            else:
+                self.sim.mcx([hq[1]], self._ancilla)
+                self.sim.mcx([hq[2]], self._ancilla)
+                self.sim.force_m(self._ancilla, True)
+                self.sim.x(hq[error_bit])
+
+
     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):
         gate = None
@@ -198,6 +243,8 @@ class QrackAceBackend:
         lq2_col = lq2 // self.row_length
         lq2_row = lq2 % self.row_length
 
+        hq1 = None
+        hq2 = None
         if (lq2_col == lq1_col) and (((lq1_row + 1) % self.row_length) == lq2_row):
             hq1 = self._unpack(lq1, True)
             hq2 = self._unpack(lq2, False)
@@ -224,6 +271,10 @@ 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.11 → pyqrack_cuda-1.44.13/pyqrack_cuda.egg-info}/PKG-INFO

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

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

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