pyqrack-cuda 1.49.5__tar.gz → 1.49.8__tar.gz

{pyqrack_cuda-1.49.5/pyqrack_cuda.egg-info → pyqrack_cuda-1.49.8}/PKG-INFO

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

{pyqrack_cuda-1.49.5 → pyqrack_cuda-1.49.8}/pyqrack/qrack_ace_backend.py

@@ -19,6 +19,12 @@ try:
 except ImportError:
     _IS_QISKIT_AVAILABLE = False
 
+_IS_QISKIT_AER_AVAILABLE = True
+try:
+    from qiskit_aer.noise import NoiseModel, depolarizing_error
+except ImportError:
+    _IS_QISKIT_AER_AVAILABLE = False
+
 
 class QrackAceBackend:
     """A back end for elided quantum error correction
@@ -62,11 +68,12 @@ class QrackAceBackend:
         self.alternating_codes = alternating_codes
         self.long_range_columns = long_range_columns
         self._is_init = [False] * qubit_count
+        self._coupling_map = None
 
         # If there's only one or zero "False" columns,
         # the entire simulator is connected, anyway.
         if (long_range_columns + 1) >= self.row_length:
-            self._is_col_long_range = [True] * long_range_columns
+            self._is_col_long_range = [True] * self.row_length
         else:
             col_seq = [True] * long_range_columns + [False]
             len_col_seq = len(col_seq)
@@ -333,12 +340,15 @@ class QrackAceBackend:
 
         if not math.isclose(ph, -lm) and not math.isclose(abs(ph), math.pi / 2):
             # Produces/destroys superposition
-            self._correct_if_like_h(th, lq)
-            self._decode(lq, hq)
+            if self._is_init[lq]:
+                self._correct_if_like_h(th, lq)
+                self.sim.mcx([hq[0]], hq[1])
             self.sim.u(hq[0], th, ph, lm)
-            if not self._is_init[lq]:
-                self.sim.u(hq[2], th, ph, lm)
-            self._encode(lq, hq)
+            self.sim.u(hq[2], th, ph, lm)
+            if self._is_init[lq]:
+                self.sim.mcx([hq[0]], hq[1])
+            else:
+                self._encode(lq, hq)
         else:
             # Shouldn't produce/destroy superposition
             for b in hq:
@@ -354,7 +364,7 @@ class QrackAceBackend:
             th -= 2 * math.pi
         while th <= -math.pi:
             th += 2 * math.pi
-        if p == Pauli.PauliY:
+        if self._is_init[lq] and (p == Pauli.PauliY):
             self._correct_if_like_h(th, lq)
 
         if (p == Pauli.PauliZ) or math.isclose(abs(th), math.pi):
@@ -363,11 +373,14 @@ class QrackAceBackend:
                 self.sim.r(p, th, b)
         else:
             # Produces/destroys superposition
-            self._decode(lq, hq)
+            if self._is_init[lq]:
+                self.sim.mcx([hq[0]], hq[1])
             self.sim.r(p, th, hq[0])
-            if not self._is_init[lq]:
-                self.sim.r(p, th, hq[2])
-            self._encode(lq, hq)
+            self.sim.r(p, th, hq[2])
+            if self._is_init[lq]:
+                self.sim.mcx([hq[0]], hq[1])
+            else:
+                self._encode(lq, hq)
 
     def h(self, lq):
         hq = self._unpack(lq)
@@ -375,11 +388,15 @@ class QrackAceBackend:
             self.sim.h(hq[0])
             return
 
-        self._decode(lq, hq)
+        self._correct(lq)
+        if self._is_init[lq]:
+            self.sim.mcx([hq[0]], hq[1])
         self.sim.h(hq[0])
-        if not self._is_init[lq]:
-            self.sim.h(hq[2])
-        self._encode(lq, hq)
+        self.sim.h(hq[2])
+        if self._is_init[lq]:
+            self.sim.mcx([hq[0]], hq[1])
+        else:
+            self._encode(lq, hq)
 
     def s(self, lq):
         hq = self._unpack(lq)
@@ -1093,6 +1110,9 @@ class QrackAceBackend:
 
     # Mostly written by Dan, but with a little help from Elara (custom OpenAI GPT)
     def get_logical_coupling_map(self):
+        if self._coupling_map:
+            return self._coupling_map
+
         coupling_map = set()
         rows, cols = self.row_length, self.col_length
 
@@ -1123,4 +1143,53 @@ class QrackAceBackend:
                         if a != b:
                             coupling_map.add((a, b))
 
-        return sorted(coupling_map)
+        self._coupling_map = sorted(coupling_map)
+
+        return self._coupling_map
+
+    # Designed by Dan, and implemented by Elara:
+    def create_noise_model(self, x=0.25, y=0.25):
+        if not _IS_QISKIT_AER_AVAILABLE:
+            raise RuntimeError(
+                "Before trying to run_qiskit_circuit() with QrackAceBackend, you must install Qiskit Aer!"
+            )
+        noise_model = NoiseModel()
+
+        for a, b in self.get_logical_coupling_map():
+            col_a, col_b = a % self.row_length, b % self.row_length
+            row_a, row_b = a // self.row_length, b // self.row_length
+            is_long_a = self._is_col_long_range[col_a]
+            is_long_b = self._is_col_long_range[col_b]
+
+            if is_long_a and is_long_b:
+                continue  # No noise on long-to-long
+
+            same_col = col_a == col_b
+            even_odd = (row_a % 2) != (row_b % 2)
+
+            if same_col and not even_odd:
+                continue  # No noise for even-even or odd-odd within a boundary column
+
+            if same_col:
+                x_cy = 1 - (1 - x)**2
+                x_swap = 1 - (1 - x)**3
+                noise_model.add_quantum_error(depolarizing_error(x, 2), 'cx', [a, b])
+                noise_model.add_quantum_error(depolarizing_error(x_cy, 2), 'cy', [a, b])
+                noise_model.add_quantum_error(depolarizing_error(x_cy, 2), 'cz', [a, b])
+                noise_model.add_quantum_error(depolarizing_error(x_swap, 2), 'swap', [a, b])
+            elif is_long_a or is_long_b:
+                y_cy = 1 - (1 - y)**2
+                y_swap = 1 - (1 - y)**3
+                noise_model.add_quantum_error(depolarizing_error(y, 2), 'cx', [a, b])
+                noise_model.add_quantum_error(depolarizing_error(y_cy, 2), 'cy', [a, b])
+                noise_model.add_quantum_error(depolarizing_error(y_cy, 2), 'cz', [a, b])
+                noise_model.add_quantum_error(depolarizing_error(y_swap, 2), 'swap', [a, b])
+            else:
+                y_cy = 1 - (1 - y)**2
+                y_swap = 1 - (1 - y)**3
+                noise_model.add_quantum_error(depolarizing_error(y_cy, 2), 'cx', [a, b])
+                noise_model.add_quantum_error(depolarizing_error(y_cy, 2), 'cy', [a, b])
+                noise_model.add_quantum_error(depolarizing_error(y_cy, 2), 'cz', [a, b])
+                noise_model.add_quantum_error(depolarizing_error(y_swap, 2), 'swap', [a, b])
+
+        return noise_model

{pyqrack_cuda-1.49.5 → pyqrack_cuda-1.49.8/pyqrack_cuda.egg-info}/PKG-INFO

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

{pyqrack_cuda-1.49.5 → pyqrack_cuda-1.49.8}/setup.py

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