pyqrack-complex128 2.1.1__tar.gz → 2.2.0__tar.gz

{pyqrack_complex128-2.1.1 → pyqrack_complex128-2.2.0}/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 7df813ab81978d06decac14b1744b19a9b0cadee; cd ..
+	git clone https://github.com/unitaryfund/qrack.git; cd qrack; git checkout 4d603e03c8f8dc7ee60e0438042b53a1343443be; cd ..
 endif
 	mkdir -p qrack/build
 ifeq ($(UNAME_S),Linux)

{pyqrack_complex128-2.1.1/pyqrack_complex128.egg-info → pyqrack_complex128-2.2.0}/PKG-INFO

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

{pyqrack_complex128-2.1.1 → pyqrack_complex128-2.2.0}/pyqrack/qrack_ace_backend.py

@@ -334,6 +334,50 @@ class QrackAceBackend:
     def clone(self):
         return QrackAceBackend(to_clone=self)
 
+    def measure_shots_consensus(self, q, s, n_instances=3, threshold=0.1):
+        # Consensus measurement across n_instances independent clones.
+        #
+        # For each shot, run n_instances clones of the current state and
+        # compare their marginal probabilities on every logical qubit via
+        # prob(). For each qubit, take the majority vote over all instances:
+        # if the average marginal >= 0.5, force_m to |1>; else to |0>.
+        #
+        # This resolves the parity ambiguity that causes XEB flip-flopping:
+        # a bit-flip on a boundary qubit appears as disagreement between
+        # instances (one sees p~0.9, another sees p~0.1). The majority vote
+        # across n_instances corrects isolated odd-parity errors without
+        # requiring a second full circuit run — just cheap prob() queries.
+        #
+        # threshold: minimum spread in marginals to trigger consensus
+        # correction (below threshold, instances agree and no correction needed).
+        n_qubits = self.num_qubits()
+        samples = []
+        for _ in range(s):
+            # Run n_instances clones, collect marginals for every qubit
+            clones = [self.clone() for _ in range(n_instances)]
+            marginals = [
+                [c.prob(lq) for lq in range(n_qubits)]
+                for c in clones
+            ]
+            # Majority vote: average marginal across instances per qubit
+            avg_marginals = [
+                sum(marginals[i][lq] for i in range(n_instances)) / n_instances
+                for lq in range(n_qubits)
+            ]
+            # Force all qubits in first clone to majority-vote outcome
+            primary = clones[0]
+            for lq in range(n_qubits):
+                result = avg_marginals[lq] >= 0.5
+                primary.force_m(lq, result)
+            # Read out requested qubits
+            _sample = primary.m_all()
+            sample = 0
+            for i in range(len(q)):
+                if (_sample >> q[i]) & 1:
+                    sample |= 1 << i
+            samples.append(sample)
+        return samples
+
     def num_qubits(self):
         return self._row_length * self._col_length
 
@@ -1473,38 +1517,48 @@ class QrackAceBackend:
         return self._coupling_map
 
     # Designed by Dan, and implemented by Elara:
-    def create_noise_model(self, x=0.25, y=0.25):
+    def create_noise_model(self, x=0.125, 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()
 
+        # Single-qubit depolarizing only on boundary qubits
+        boundary_qubits = set()
+        for a, b in self.get_logical_coupling_map():
+            col_a = a % self._row_length
+            col_b = 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 not (is_long_a and is_long_b):
+                if not is_long_a:
+                    boundary_qubits.add(a)
+                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, col_b = a % self._row_length, b % self._row_length
-            row_a, row_b = a // self._row_length, b // self._row_length
+            col_a = a % self._row_length
+            col_b = 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
-
-            if (col_a == col_b) or (row_a == row_b):
-                continue  # No noise for same column
-
-            if 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])
+                continue
+
+            p2 = 1 - (1 - y) ** 2 if (is_long_a or is_long_b) else y
+            p3 = 1 - (1 - y) ** 3 if (is_long_a or is_long_b) else 1 - (1 - y) ** 2
+
+            for gate in ["cx", "cy", "cz"]:
+                noise_model.add_quantum_error(depolarizing_error(p2, 2), gate, [a, b])
+            for gate in ["swap", "iswap"]:
+                noise_model.add_quantum_error(depolarizing_error(p3, 2), gate, [a, b])
 
         return noise_model

{pyqrack_complex128-2.1.1 → pyqrack_complex128-2.2.0}/pyqrack/qrack_simulator.py

@@ -554,8 +554,8 @@ class QrackSimulator:
         Qrack.qrack_lib.MCAdjT(self.sid, len(c), QrackSimulator._ulonglong_byref(c), q)
         self._throw_if_error()
 
-    def mcu(self, c, q, th, ph, la):
-        """Multi-controlled arbitraty unitary
+    def mcu(self, c, q, th, ph, la, gm=0.0):
+        """Multi-controlled arbitrary unitary
 
         If all controlled qubits are `|1>` then the unitary gate described by
         parameters is applied to the target qubit.
@@ -566,6 +566,7 @@ class QrackSimulator:
             th: theta
             ph: phi
             la: lambda
+            gm: gamma
 
         Raises:
             RuntimeError: QrackSimulator raised an exception.
@@ -578,6 +579,7 @@ class QrackSimulator:
             ctypes.c_double(th),
             ctypes.c_double(ph),
             ctypes.c_double(la),
+            ctypes.c_double(gm),
         )
         self._throw_if_error()
 
@@ -736,8 +738,8 @@ class QrackSimulator:
         Qrack.qrack_lib.MACAdjT(self.sid, len(c), QrackSimulator._ulonglong_byref(c), q)
         self._throw_if_error()
 
-    def macu(self, c, q, th, ph, la):
-        """Anti multi-controlled arbitraty unitary
+    def macu(self, c, q, th, ph, la, gm=0.0):
+        """Anti multi-controlled arbitrary unitary
 
         If all controlled qubits are `|0>` then the unitary gate described by
         parameters is applied to the target qubit.
@@ -748,6 +750,7 @@ class QrackSimulator:
             th: theta
             ph: phi
             la: lambda
+            gm: gamma
 
         Raises:
             RuntimeError: QrackSimulator raised an exception.
@@ -760,11 +763,12 @@ class QrackSimulator:
             ctypes.c_double(th),
             ctypes.c_double(ph),
             ctypes.c_double(la),
+            ctypes.c_double(gm),
         )
         self._throw_if_error()
 
     def macmtrx(self, c, m, q):
-        """Anti multi-controlled arbitraty operator
+        """Anti multi-controlled arbitrary operator
 
         If all controlled qubits are `|0>` then the arbitrary operation by
         parameters is applied to the target qubit.
@@ -4194,7 +4198,7 @@ class QrackSimulator:
                 float(operation.params[0]),
                 float(operation.params[1]),
             )
-        elif (name == "cu3") or (name == "cu"):
+        elif name == "cu3":
             self._sim.mcu(
                 [q._index for q in operation.qubits[0:1]],
                 operation.qubits[1]._index,
@@ -4202,6 +4206,15 @@ class QrackSimulator:
                 float(operation.params[1]),
                 float(operation.params[2]),
             )
+        elif name == "cu":
+            self._sim.mcu(
+                [q._index for q in operation.qubits[0:1]],
+                operation.qubits[1]._index,
+                float(operation.params[0]),
+                float(operation.params[1]),
+                float(operation.params[2]),
+                float(operation.params[3]),
+            )
         elif name == "cx":
             self._sim.mcx([q._index for q in operation.qubits[0:1]], operation.qubits[1]._index)
         elif name == "cy":

{pyqrack_complex128-2.1.1 → pyqrack_complex128-2.2.0}/pyqrack/qrack_system/qrack_system.py

@@ -629,6 +629,7 @@ class QrackSystem:
             c_double,
             c_double,
             c_double,
+            c_double,
         ]
 
         self.qrack_lib.MCMtrx.restype = None
@@ -715,6 +716,7 @@ class QrackSystem:
             c_double,
             c_double,
             c_double,
+            c_double,
         ]
 
         self.qrack_lib.MACMtrx.restype = None

{pyqrack_complex128-2.1.1 → pyqrack_complex128-2.2.0/pyqrack_complex128.egg-info}/PKG-INFO

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

{pyqrack_complex128-2.1.1 → pyqrack_complex128-2.2.0}/setup.py

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