pyqrack-cuda-complex128 2.2.0__tar.gz → 2.2.2__tar.gz

{pyqrack_cuda_complex128-2.2.0 → pyqrack_cuda_complex128-2.2.2}/Makefile

@@ -30,7 +30,7 @@ build-deps:
 	rm -rf pyqrack/qrack_system/qrack_cl_precompile
 ifneq ($(OS),Windows_NT)
 ifeq ($(QRACK_PRESENT),)
-	git clone https://github.com/unitaryfund/qrack.git; cd qrack; git checkout 4d603e03c8f8dc7ee60e0438042b53a1343443be; cd ..
+	git clone https://github.com/unitaryfund/qrack.git; cd qrack; git checkout ccb9cca0726e48c598696b51a66fb984ff8ab0f1; cd ..
 endif
 	mkdir -p qrack/build
 ifeq ($(UNAME_S),Linux)

{pyqrack_cuda_complex128-2.2.0/pyqrack_cuda_complex128.egg-info → pyqrack_cuda_complex128-2.2.2}/PKG-INFO

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

{pyqrack_cuda_complex128-2.2.0 → pyqrack_cuda_complex128-2.2.2}/pyqrack/qrack_ace_backend.py

@@ -326,14 +326,18 @@ class QrackAceBackend:
                 )
             )
 
-            # You can still "monkey-patch" this, after the constructor.
-            if "QRACK_QUNIT_SEPARABILITY_THRESHOLD" not in os.environ:
-                # (1 - 1 / sqrt(2)) / 4 (but empirically tuned)
-                self.sim[i].set_sdrp(0.073223304703363119)
+            # You can still call "set_sdrp," after the constructor.
+            # if "QRACK_QUNIT_SEPARABILITY_THRESHOLD" not in os.environ:
+            #     # (1 - 1 / sqrt(2)) / 4 (but empirically tuned)
+            #     self.sim[i].set_sdrp(0.073223304703363119)
 
     def clone(self):
         return QrackAceBackend(to_clone=self)
 
+    def set_sdrp(self, sdrp):
+        for sim in self.sim:
+            sim.set_sdrp(sdrp)
+
     def measure_shots_consensus(self, q, s, n_instances=3, threshold=0.1):
         # Consensus measurement across n_instances independent clones.
         #
@@ -590,7 +594,12 @@ class QrackAceBackend:
                 )
                 / 7
             )
-            result = ((prms + (1 - qrms)) / 2) >= 0.5
+            eff_prob = (prms + (1 - qrms)) / 2
+            result = (
+                (random.random() < 0.5)
+                if (eff_prob - 0.5) == 0.0
+                else (eff_prob > 0.5)
+            )
             syndrome = (
                 [1 - p[0], 1 - p[1], 1 - p[2], 1 - p[3], 1 - p[4]]
                 if result
@@ -639,7 +648,12 @@ class QrackAceBackend:
             # 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)
-            result = ((prms + (1 - qrms)) / 2) >= 0.5
+            eff_prob = (prms + (1 - qrms)) / 2
+            result = (
+                (random.random() < 0.5)
+                if (eff_prob - 0.5) == 0.0
+                else (eff_prob > 0.5)
+            )
             syndrome = [1 - p[0], 1 - p[1], 1 - p[2]] if result else [p[0], p[1], p[2]]
             for q in range(3):
                 if syndrome[q] > (0.5 + self._epsilon):
@@ -1125,12 +1139,12 @@ class QrackAceBackend:
             else:
                 self.sim[b[0]].force_m(b[1], result)
 
-        b = hq[1]
+        b = hq[lhv]
         b.reset()
         if result:
             b.x()
 
-        return c
+        return result
 
     def m_all(self):
         # Randomize the order of measurement to amortize error.
@@ -1495,22 +1509,25 @@ class QrackAceBackend:
             return self._coupling_map
 
         coupling_map = set()
-        rows, cols = self._row_length, self._col_length
+        # self._row_length is the number of qubits per row (the row's span,
+        # i.e. the column count); self._col_length is the number of rows.
+        num_rows, num_cols = self._col_length, self._row_length
 
-        # Map each column index to its full list of logical qubit indices
+        # Row-major logical index, strided by the actual row span.
         def logical_index(row, col):
-            return row * cols + col
+            return row * num_cols + col
 
-        for col in range(cols):
-            connected_cols, _ = self._get_connected(col, False)
-            for row in range(rows):
-                connected_rows, _ = self._get_connected(row, False)
+        for row in range(num_rows):
+            connected_rows, _ = self._get_connected(row, True)
+            for col in range(num_cols):
+                connected_cols, _ = self._get_connected(col, False)
                 a = logical_index(row, col)
-                for c in connected_cols:
-                    for r in connected_rows:
+                for r in connected_rows:
+                    for c in connected_cols:
                         b = logical_index(r, c)
                         if a != b:
                             coupling_map.add((a, b))
+                            coupling_map.add((b, a))
 
         self._coupling_map = sorted(coupling_map)
 
@@ -1524,9 +1541,11 @@ class QrackAceBackend:
             )
         noise_model = NoiseModel()
 
-        # Single-qubit depolarizing only on boundary qubits
         boundary_qubits = set()
         for a, b in self.get_logical_coupling_map():
+            if a > b:
+                #Skip duplicates
+                continue
             col_a = a % self._row_length
             col_b = b % self._row_length
             is_long_a = self._is_col_long_range[col_a]
@@ -1537,12 +1556,6 @@ class QrackAceBackend:
                 if not is_long_b:
                     boundary_qubits.add(b)
 
-        for q in boundary_qubits:
-            for gate in ["u", "u1", "u2", "u3", "h", "x", "y", "z",
-                         "s", "sdg", "t", "tdg", "rx", "ry", "rz"]:
-                noise_model.add_quantum_error(
-                    depolarizing_error(x, 1), gate, [q])
-
         # Two-qubit depolarizing on boundary-crossing and boundary-adjacent gates
         for a, b in self.get_logical_coupling_map():
             col_a = a % self._row_length

{pyqrack_cuda_complex128-2.2.0 → pyqrack_cuda_complex128-2.2.2/pyqrack_cuda_complex128.egg-info}/PKG-INFO

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

{pyqrack_cuda_complex128-2.2.0 → pyqrack_cuda_complex128-2.2.2}/setup.py

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