pyqrack-cuda 1.55.0__tar.gz → 1.55.2__tar.gz

{pyqrack_cuda-1.55.0/pyqrack_cuda.egg-info → pyqrack_cuda-1.55.2}/PKG-INFO

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

{pyqrack_cuda-1.55.0 → pyqrack_cuda-1.55.2}/pyqrack/qrack_ace_backend.py

@@ -121,6 +121,7 @@ class LHVQubit:
             self.x()
         return result
 
+
 # Provided by Elara (the custom OpenAI GPT)
 def _cpauli_lhv(prob, targ, axis, anti, theta=math.pi):
     """
@@ -143,6 +144,7 @@ def _cpauli_lhv(prob, targ, axis, anti, theta=math.pi):
     elif axis == Pauli.PauliZ:
         targ.rz(effective_theta)
 
+
 class QrackAceBackend:
     """A back end for elided quantum error correction
 
@@ -186,6 +188,16 @@ class QrackAceBackend:
         self.long_range_columns = long_range_columns
         self.is_transpose = is_transpose
 
+        fppow = 5
+        if "QRACK_FPPOW" in os.environ:
+            fppow = int(os.environ.get("QRACK_FPPOW"))
+        if fppow < 5:
+            self._epsilon = 2**-11
+        elif fppow > 5:
+            self._epsilon = 2**-53
+        else:
+            self._epsilon = 2**-24
+
         self._coupling_map = None
 
         # If there's only one or zero "False" columns,
@@ -218,7 +230,9 @@ class QrackAceBackend:
                     tot_qubits += 1
                     sim_counts[sim_id] += 1
 
-                    self._lhv_dict[tot_qubits] = LHVQubit(toClone = toClone._lhv_dict[tot_qubits] if toClone else None)
+                    self._lhv_dict[tot_qubits] = LHVQubit(
+                        toClone=toClone._lhv_dict[tot_qubits] if toClone else None
+                    )
                     tot_qubits += 1
 
                     sim_id = (sim_id + 1) % sim_count
@@ -262,7 +276,9 @@ class QrackAceBackend:
             col_len -= 1
         row_len = width // col_len
 
-        self._col_length, self._row_length = (row_len, col_len) if is_transpose else (col_len, row_len)
+        self._col_length, self._row_length = (
+            (row_len, col_len) if is_transpose else (col_len, row_len)
+        )
 
     def _ct_pair_prob(self, q1, q2):
         p1 = self.sim[q1[0]].prob(q1[1]) if isinstance(q1, tuple) else q1.prob()
@@ -356,7 +372,12 @@ class QrackAceBackend:
             b = hq[2]
             self.sim[b[0]].h(b[1])
 
-        p = [self.sim[hq[0][0]].prob(hq[0][1]), hq[1].prob(), self.sim[hq[2][0]].prob(hq[2][1])]
+        # RMS
+        p = [
+            self.sim[hq[0][0]].prob(hq[0][1]),
+            hq[1].prob(),
+            self.sim[hq[2][0]].prob(hq[2][1]),
+        ]
         # Balancing suggestion from Elara (the custom OpenAI GPT)
         prms = math.sqrt((p[0] ** 2 + p[1] ** 2 + p[2] ** 2) / 3)
         qrms = math.sqrt(((1 - p[0]) ** 2 + (1 - p[1]) ** 2 + (1 - p[2]) ** 2) / 3)
@@ -565,14 +586,16 @@ class QrackAceBackend:
         hq2 = self._unpack(lq2)
 
         if lq1_lr and lq2_lr:
-            connected = (lq1_col == lq2_col) or ((self.long_range_columns + 1) >= self._row_length)
+            connected = (lq1_col == lq2_col) or (
+                (self.long_range_columns + 1) >= self._row_length
+            )
             c = (lq1_col - 1) % self._row_length
             while not connected and self._is_col_long_range[c]:
-                connected = (lq2_col == c)
+                connected = lq2_col == c
                 c = (c - 1) % self._row_length
             c = (lq1_col + 1) % self._row_length
             while not connected and self._is_col_long_range[c]:
-                connected = (lq2_col == c)
+                connected = lq2_col == c
                 c = (c + 1) % self._row_length
 
             b1 = hq1[0]
@@ -766,15 +789,29 @@ class QrackAceBackend:
             b = hq[0]
             return self.sim[b[0]].m(b[1])
 
-        result = (random.random() < self.prob(lq))
+        p = self.prob(lq)
+        result = ((p + self._epsilon) >= 1) or (random.random() < p)
+
         b = hq[0]
-        self.sim[b[0]].force_m(b[1], result)
+        p = self.sim[b[0]].prob(b[1]) if result else (1 - self.sim[b[0]].prob(b[1]))
+        if p < self._epsilon:
+            if self.sim[b[0]].m(b[1]):
+                self.sim[b[0]].x(b[1])
+        else:
+            self.sim[b[0]].force_m(b[1], result)
+
         b = hq[1]
         b.reset()
         if result:
             b.x()
+
         b = hq[2]
-        self.sim[b[0]].force_m(b[1], result)
+        p = self.sim[b[0]].prob(b[1]) if result else (1 - self.sim[b[0]].prob(b[1]))
+        if p < self._epsilon:
+            if self.sim[b[0]].m(b[1]):
+                self.sim[b[0]].x(b[1])
+        else:
+            self.sim[b[0]].force_m(b[1], result)
 
         return result
 

{pyqrack_cuda-1.55.0 → pyqrack_cuda-1.55.2}/pyqrack/qrack_simulator.py

@@ -181,12 +181,12 @@ class QrackSimulator:
         self._throw_if_error()
 
     def set_concurrency(self, p):
-        """ Set the CPU parallel thread count"""
+        """Set the CPU parallel thread count"""
         Qrack.qrack_lib.set_concurrency(self.sid, p)
         self._throw_if_error()
 
     def set_device(self, d):
-        """ Set the GPU device ID"""
+        """Set the GPU device ID"""
         Qrack.qrack_lib.set_device(self.sid, d)
         self._throw_if_error()
 

{pyqrack_cuda-1.55.0 → pyqrack_cuda-1.55.2/pyqrack_cuda.egg-info}/PKG-INFO

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

{pyqrack_cuda-1.55.0 → pyqrack_cuda-1.55.2}/setup.py

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