pyqrack-cuda 1.44.29__tar.gz → 1.44.31__tar.gz

{pyqrack_cuda-1.44.29/pyqrack_cuda.egg-info → pyqrack_cuda-1.44.31}/PKG-INFO

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

{pyqrack_cuda-1.44.29 → pyqrack_cuda-1.44.31}/pyqrack/qrack_ace_backend.py

@@ -163,27 +163,50 @@ class QrackAceBackend:
         single_bit_value = self.sim.prob(hq[single_bit])
         single_bit_polarization = max(single_bit_value, 1 - single_bit_value)
         samples = self.sim.measure_shots([hq[other_bits[0]], hq[other_bits[1]]], shots)
-        syndrome_indices = [other_bits[1], other_bits[0]] if (single_bit_value >= 0.5) else [other_bits[0], other_bits[1]]
+        syndrome_indices = (
+            [other_bits[1], other_bits[0]]
+            if (single_bit_value >= 0.5)
+            else [other_bits[0], other_bits[1]]
+        )
         syndrome = [0, 0, 0]
+        values = []
         for sample in samples:
             match sample:
                 case 0:
-                    syndrome[single_bit] += single_bit_value
+                    value = single_bit_value
+                    syndrome[single_bit] += value
                 case 1:
-                    syndrome[syndrome_indices[0]] += single_bit_polarization
+                    value = single_bit_polarization
+                    syndrome[syndrome_indices[0]] += value
                 case 2:
-                    syndrome[syndrome_indices[1]] += single_bit_polarization
+                    value = single_bit_polarization
+                    syndrome[syndrome_indices[1]] += value
                 case 3:
-                    syndrome[single_bit] += 1 - single_bit_value
+                    value = 1 - single_bit_value
+                    syndrome[single_bit] += value
+            values.append(value)
 
         # Suggestion from Elara (custom OpenAI GPT):
-        # Assume binomial statistics and compute the standard deviation.
-        # Only correct if we're outside a confidence interval.
+        # Compute the standard deviation and only correct if we're outside a confidence interval.
         # (This helps avoid limit-point over-correction.)
-        syndrome_std_dev = (sum(syndrome) - shots / 2) / math.sqrt(shots)
+        syndrome_sum = sum(syndrome)
+        z_score_numer = syndrome_sum - shots / 2
+        z_score = 0
+        if z_score_numer > 0:
+            syndrome_mean = syndrome_sum / shots
+            syndrome_variance = (
+                sum((value - syndrome_mean) ** 2 for value in values) / shots
+            )
+            z_score_denom = math.sqrt(syndrome_variance)
+            z_score = (
+                math.inf
+                if math.isclose(z_score_denom, 0)
+                else (z_score_numer / z_score_denom)
+            )
+
         force_syndrome = True
-        # (~68% confidence interval, but this is the margin for average sampling error.)
-        if syndrome_std_dev > 1:
+        # (From Elara, this is the value that minimizes the sum of Type I and Type II error.)
+        if z_score >= (497 / 999):
             # There is an error.
             error_bit = syndrome.index(max(syndrome))
             if error_bit == single_bit:

{pyqrack_cuda-1.44.29 → pyqrack_cuda-1.44.31/pyqrack_cuda.egg-info}/PKG-INFO

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

{pyqrack_cuda-1.44.29 → pyqrack_cuda-1.44.31}/setup.py

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