pyqrack-cuda 1.45.1__tar.gz → 1.46.0__tar.gz

{pyqrack_cuda-1.45.1/pyqrack_cuda.egg-info → pyqrack_cuda-1.46.0}/PKG-INFO

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

{pyqrack_cuda-1.45.1 → pyqrack_cuda-1.46.0}/pyqrack/qrack_ace_backend.py

@@ -41,20 +41,41 @@ class QrackAceBackend:
     def __init__(
         self,
         qubit_count=1,
+        recursive_stack_depth=1,
         alternating_codes=True,
         isTensorNetwork=False,
+        isStabilizerHybrid=False,
+        isBinaryDecisionTree=False,
         toClone=None,
     ):
-        self.sim = (
-            toClone.sim.clone()
-            if toClone
-            else QrackSimulator(3 * qubit_count + 1, isTensorNetwork=isTensorNetwork)
-        )
+        if recursive_stack_depth < 1:
+            recursive_stack_depth = 1
+        if toClone:
+            qubit_count = toClone.num_qubits()
+        if recursive_stack_depth > 1:
+            recursive_stack_depth -= 1
+            self.sim = (
+                toClone.sim
+                if toClone
+                else QrackAceBackend(3 * qubit_count + 1, recursive_stack_depth=recursive_stack_depth, alternating_codes=alternating_codes, isTensorNetwork=isTensorNetwork, isStabilizerHybrid=isStabilizerHybrid, isBinaryDecisionTree=isBinaryDecisionTree)
+            )
+        else:
+            self.sim = (
+                toClone.sim.clone()
+                if toClone
+                else QrackSimulator(3 * qubit_count + 1, isTensorNetwork=isTensorNetwork, isStabilizerHybrid=isStabilizerHybrid, isBinaryDecisionTree=isBinaryDecisionTree)
+            )
         self._ancilla = 3 * qubit_count
         self._factor_width(qubit_count)
         self.alternating_codes = alternating_codes
         self._is_init = [False] * qubit_count
 
+    def clone(self):
+        return QrackAceBackend(toClone=self)
+
+    def num_qubits(self):
+        return self.sim.num_qubits() // 3
+
     def _factor_width(self, width):
         col_len = math.floor(math.sqrt(width))
         while ((width // col_len) * col_len) != width:
@@ -244,7 +265,7 @@ class QrackAceBackend:
         if not math.isclose(th, 0):
             self._correct(lq)
 
-    def u(self, th, ph, lm, lq):
+    def u(self, lq, th, ph, lm):
         while ph > math.pi:
             ph -= 2 * math.pi
         while ph <= -math.pi:
@@ -404,6 +425,36 @@ class QrackAceBackend:
     def acz(self, lq1, lq2):
         self._cpauli(lq1, lq2, True, Pauli.PauliZ)
 
+    def mcx(self, lq1, lq2):
+        if len(lq1) > 1:
+            raise RuntimeError("QrackAceBackend.mcx() is provided for syntax convenience and only supports 1 control qubit!")
+        self._cpauli(lq1[0], lq2, False, Pauli.PauliX)
+
+    def mcy(self, lq1, lq2):
+        if len(lq1) > 1:
+            raise RuntimeError("QrackAceBackend.mcy() is provided for syntax convenience and only supports 1 control qubit!")
+        self._cpauli(lq1[0], lq2, False, Pauli.PauliY)
+
+    def mcz(self, lq1, lq2):
+        if len(lq1) > 1:
+            raise RuntimeError("QrackAceBackend.mcz() is provided for syntax convenience and only supports 1 control qubit!")
+        self._cpauli(lq1[0], lq2, False, Pauli.PauliZ)
+
+    def macx(self, lq1, lq2):
+        if len(lq1) > 1:
+            raise RuntimeError("QrackAceBackend.macx() is provided for syntax convenience and only supports 1 control qubit!")
+        self._cpauli(lq1[0], lq2, True, Pauli.PauliX)
+
+    def macy(self, lq1, lq2):
+        if len(lq1) > 1:
+            raise RuntimeError("QrackAceBackend.macy() is provided for syntax convenience and only supports 1 control qubit!")
+        self._cpauli(lq1[0], lq2, True, Pauli.PauliY)
+
+    def macz(self, lq1, lq2):
+        if len(lq1) > 1:
+            raise RuntimeError("QrackAceBackend.macz() is provided for syntax convenience and only supports 1 control qubit!")
+        self._cpauli(lq1[0], lq2, True, Pauli.PauliZ)
+
     def swap(self, lq1, lq2):
         self.cx(lq1, lq2)
         self.cx(lq2, lq1)
@@ -437,7 +488,7 @@ class QrackAceBackend:
             syndrome += self.sim.m(hq[q])
         # The two separable parts of the code are correlated,
         # but not non-locally, via entanglement.
-        # Prefer to collapse the analytical part toward agreement.
+        # Collapse the other separable part toward agreement.
         if syndrome == 0:
             self.sim.force_m(hq[single_bit], False)
         elif syndrome == 2:
@@ -448,12 +499,22 @@ class QrackAceBackend:
 
         return True if (syndrome > 1) else False
 
+    def force_m(self, lq, c):
+        hq = self._unpack(lq)
+        self._correct(lq)
+        for q in hq:
+            self.sim.force_m(q, c)
+        self._is_init[lq] = False
+
+        return c
+
     def m_all(self):
         result = 0
         # Whenever a nonzero syndrome occurs (so the code has an error),
         # we insert the more-probable results and collapse towards it.
         # Randomize the order of post-selection to amortize error.
-        lqubits = random.shuffle(list(range(self.sim.num_qubits() // 3)))
+        lqubits = list(range(self.sim.num_qubits() // 3))
+        random.shuffle(lqubits)
         for lq in lqubits:
             if self.m(lq):
                 result |= 1 << lq
@@ -463,12 +524,14 @@ class QrackAceBackend:
     def measure_shots(self, q, s, high_accuracy=True):
         if high_accuracy:
             samples = []
+            _q = q.copy()
             for _ in range(s):
                 clone = self.sim.clone()
                 sample = 0
-                for i in range(len(q)):
-                    if clone.m(q[i]):
-                        sample |= 1 << i
+                random.shuffle(_q)
+                for b in _q:
+                    if clone.m(b):
+                        sample |= 1 << q.index(b)
                 samples.append(sample)
 
             return samples
@@ -497,6 +560,20 @@ class QrackAceBackend:
 
         return results
 
+    def prob(self, lq):
+        self._correct(lq)
+        hq = self._unpack(lq)
+        even_row = not ((lq // self.row_length) & 1)
+        if not self.alternating_codes or even_row:
+            other_bits = [0, 1]
+        else:
+            other_bits = [1, 2]
+        self.sim.mcx([other_bits[0]], other_bits[1])
+        result = self.sim.prob(hq[1])
+        self.sim.mcx([other_bits[0]], other_bits[1])
+
+        return result
+
     def _apply_op(self, operation):
         name = operation.name
 
@@ -520,27 +597,27 @@ class QrackAceBackend:
                     return
 
         if (name == "u1") or (name == "p"):
-            self._sim.u(0, 0, float(operation.params[0]), operation.qubits[0]._index)
+            self._sim.u(operation.qubits[0]._index, 0, 0, float(operation.params[0]))
         elif name == "u2":
             self._sim.u(
+                operation.qubits[0]._index,
                 math.pi / 2,
                 float(operation.params[0]),
                 float(operation.params[1]),
-                operation.qubits[0]._index,
             )
         elif (name == "u3") or (name == "u"):
             self._sim.u(
+                operation.qubits[0]._index,
                 float(operation.params[0]),
                 float(operation.params[1]),
                 float(operation.params[2]),
-                operation.qubits[0]._index,
             )
         elif name == "r":
             self._sim.u(
+                operation.qubits[0]._index,
                 float(operation.params[0]),
                 float(operation.params[1]) - math.pi / 2,
                 (-1 * float(operation.params[1])) + math.pi / 2,
-                operation.qubits[0]._index,
             )
         elif name == "rx":
             self._sim.r(

{pyqrack_cuda-1.45.1 → pyqrack_cuda-1.46.0/pyqrack_cuda.egg-info}/PKG-INFO

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

{pyqrack_cuda-1.45.1 → pyqrack_cuda-1.46.0}/setup.py

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