pyqrack-cpu-complex128 1.58.0__py3-none-win_amd64.whl → 1.58.10__py3-none-win_amd64.whl

pyqrack/qrack_ace_backend.py

@@ -102,6 +102,45 @@ class LHVQubit:
         self.ry(theta)
         self.rz(phi)
 
+    # Provided verbatim by Elara (the custom OpenAI GPT):
+    def mtrx(self, matrix):
+        """
+        Apply a 2x2 unitary matrix to the LHV Bloch vector using only standard math/cmath.
+        Matrix format: [a, b, c, d] for [[a, b], [c, d]]
+        """
+        a, b, c, d = matrix
+
+        # Current Bloch vector
+        x, y, z = self.bloch
+
+        # Convert to density matrix ρ = ½ (I + xσx + yσy + zσz)
+        rho = [[(1 + z) / 2, (x - 1j * y) / 2], [(x + 1j * y) / 2, (1 - z) / 2]]
+
+        # Compute U * ρ
+        u_rho = [
+            [a * rho[0][0] + b * rho[1][0], a * rho[0][1] + b * rho[1][1]],
+            [c * rho[0][0] + d * rho[1][0], c * rho[0][1] + d * rho[1][1]],
+        ]
+
+        # Compute (U * ρ) * U†
+        rho_prime = [
+            [
+                u_rho[0][0] * a.conjugate() + u_rho[0][1] * b.conjugate(),
+                u_rho[0][0] * c.conjugate() + u_rho[0][1] * d.conjugate(),
+            ],
+            [
+                u_rho[1][0] * a.conjugate() + u_rho[1][1] * b.conjugate(),
+                u_rho[1][0] * c.conjugate() + u_rho[1][1] * d.conjugate(),
+            ],
+        ]
+
+        # Extract Bloch components: Tr(ρ'σi) = 2 * Re[...]
+        new_x = 2 * rho_prime[0][1].real + 2 * rho_prime[1][0].real
+        new_y = 2 * (rho_prime[0][1].imag - rho_prime[1][0].imag)
+        new_z = 2 * rho_prime[0][0].real - 1  # since Tr(ρ') = 1
+
+        self.bloch = [new_x, new_y, new_z]
+
     def prob(self, basis=Pauli.PauliZ):
         """Sample a classical outcome from the current 'quantum' state"""
         if basis == Pauli.PauliZ:
@@ -404,7 +443,91 @@ class QrackAceBackend:
 
         return qb, lhv
 
-    def _correct(self, lq, phase=False):
+    def _get_lhv_bloch_angles(self, sim):
+        # Z axis
+        z = 1 - 2 * sim.prob(Pauli.PauliZ)
+        prob = z**2
+
+        # X axis
+        x = 1 - 2 * sim.prob(Pauli.PauliX)
+        prob += x**2
+
+        # Y axis
+        y = 1 - 2 * sim.prob(Pauli.PauliY)
+        prob += y**2
+
+        prob = math.sqrt(prob)
+        inclination = math.atan2(math.sqrt(x**2 + y**2), z)
+        azimuth = math.atan2(y, x)
+
+        return prob, azimuth, inclination
+
+    def _get_bloch_angles(self, hq):
+        sim = self.sim[hq[0]]
+        q = hq[1]
+
+        # Z axis
+        z = 1 - 2 * sim.prob(q)
+        prob = z**2
+
+        # X axis
+        sim.h(q)
+        x = 1 - 2 * sim.prob(q)
+        prob += x**2
+        sim.h(q)
+
+        # Y axis
+        sim.adjs(q)
+        sim.h(q)
+        y = 1 - 2 * sim.prob(q)
+        prob += y**2
+        sim.h(q)
+        sim.s(q)
+
+        prob = math.sqrt(prob)
+        inclination = math.atan2(math.sqrt(x**2 + y**2), z)
+        azimuth = math.atan2(y, x)
+
+        return prob, azimuth, inclination
+
+    def _rotate_to_bloch(
+        self, hq, delta_azimuth, delta_inclination
+    ):
+        sim = self.sim[hq[0]]
+        q = hq[1]
+
+        # Apply rotation as "Azimuth, Inclination" (AI)
+        cosA = math.cos(delta_azimuth)
+        sinA = math.sin(delta_azimuth)
+        cosI = math.cos(delta_inclination / 2)
+        sinI = math.sin(delta_inclination / 2)
+
+        m00 = complex(cosI, 0)
+        m01 = complex(-cosA, sinA) * sinI
+        m10 = complex(cosA, sinA) * sinI
+        m11 = complex(cosI, 0)
+
+        sim.mtrx([m00, m01, m10, m11], q)
+
+
+    def _rotate_lhv_to_bloch(
+        self, sim, delta_azimuth, delta_inclination
+    ):
+        # Apply rotation as "Azimuth, Inclination" (AI)
+        cosA = math.cos(delta_azimuth)
+        sinA = math.sin(delta_azimuth)
+        cosI = math.cos(delta_inclination / 2)
+        sinI = math.sin(delta_inclination / 2)
+
+        m00 = complex(cosI, 0)
+        m01 = complex(-cosA, sinA) * sinI
+        m10 = complex(cosA, sinA) * sinI
+        m11 = complex(cosI, 0)
+
+        sim.mtrx([m00, m01, m10, m11])
+
+
+    def _correct(self, lq, phase=False, skip_rotation=False):
         hq = self._unpack(lq)
 
         if len(hq) == 1:
@@ -454,6 +577,34 @@ class QrackAceBackend:
                         hq[q].x()
                     else:
                         self.sim[hq[q][0]].x(hq[q][1])
+
+            if not skip_rotation:
+                p, a, i = [0, 0, 0, 0, 0], [0, 0, 0, 0, 0], [0, 0, 0, 0, 0]
+                p[0], a[0], i[0] = self._get_bloch_angles(hq[0])
+                p[1], a[1], i[1] = self._get_bloch_angles(hq[1])
+                p[3], a[3], i[3] = self._get_bloch_angles(hq[3])
+                p[4], a[4], i[4] = self._get_bloch_angles(hq[4])
+
+                indices = []
+                a_target = 0
+                i_target = 0
+                weight = 0
+                for x in range(5):
+                    if p[x] < 0.5:
+                        continue
+                    indices.append(x)
+                    w = (1.5 - p[x])
+                    w *= w
+                    a_target += w * a[x]
+                    i_target += w * i[x]
+                    weight += w
+
+                if len(indices) > 1:
+                    a_target /= weight
+                    i_target /= weight
+                    for x in indices:
+                        self._rotate_to_bloch(hq[x], a_target - a[x], i_target - i[x])
+
         else:
             # RMS
             p = [
@@ -473,6 +624,31 @@ class QrackAceBackend:
                     else:
                         self.sim[hq[q][0]].x(hq[q][1])
 
+            if not skip_rotation:
+                p, a, i = [0, 0, 0], [0, 0, 0], [0, 0, 0]
+                p[0], a[0], i[0] = self._get_bloch_angles(hq[0])
+                p[1], a[1], i[1] = self._get_bloch_angles(hq[1])
+
+                indices = []
+                a_target = 0
+                i_target = 0
+                weight = 0
+                for x in range(3):
+                    if p[x] < 0.5:
+                        continue
+                    indices.append(x)
+                    w = (1.5 - p[x])
+                    w *= w
+                    a_target += w * a[x]
+                    i_target += w * i[x]
+                    weight += w
+
+                if len(indices) > 1:
+                    a_target /= weight
+                    i_target /= weight
+                    for x in indices:
+                        self._rotate_to_bloch(hq[x], a_target - a[x], i_target - i[x])
+
         if phase:
             for q in qb:
                 b = hq[q]
@@ -496,8 +672,8 @@ class QrackAceBackend:
         b = hq[lhv]
         b.u(th, ph, lm)
 
-        self._correct(lq, False)
-        self._correct(lq, True)
+        self._correct(lq, False, True)
+        self._correct(lq, True, False)
 
     def r(self, p, th, lq):
         hq = self._unpack(lq)
@@ -521,7 +697,7 @@ class QrackAceBackend:
             b.rz(th)
 
         if p != Pauli.PauliZ:
-            self._correct(lq, False)
+            self._correct(lq, False, p != Pauli.PauliX)
         if p != Pauli.PauliX:
             self._correct(lq, True)
 
@@ -745,7 +921,7 @@ class QrackAceBackend:
 
         self._correct(lq1, True)
         if pauli != Pauli.PauliZ:
-            self._correct(lq2, False)
+            self._correct(lq2, False, pauli != Pauli.PauliX)
         if pauli != Pauli.PauliX:
             self._correct(lq2, True)
 

pyqrack/qrack_system/qrack_system.py

@@ -147,7 +147,11 @@ class QrackSystem:
         self.qrack_lib.set_device.argtypes = [c_ulonglong, c_longlong]
 
         self.qrack_lib.set_device_list.restype = None
-        self.qrack_lib.set_device_list.argtypes = [c_ulonglong, c_ulonglong, POINTER(c_longlong)]
+        self.qrack_lib.set_device_list.argtypes = [
+            c_ulonglong,
+            c_ulonglong,
+            POINTER(c_longlong),
+        ]
 
         # pseudo-quantum
 

{pyqrack_cpu_complex128-1.58.0.dist-info → pyqrack_cpu_complex128-1.58.10.dist-info}/METADATA

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

{pyqrack_cpu_complex128-1.58.0.dist-info → pyqrack_cpu_complex128-1.58.10.dist-info}/RECORD

@@ -1,7 +1,7 @@
 pyqrack/__init__.py,sha256=9fGCYdEUg_AsENCDRylbrZvUMVbUajRTn3CyuFGuFyM,805
 pyqrack/neuron_activation_fn.py,sha256=GOqcCkiEB60jCojTrcuHuZMDP5aTOy0adtR8SY_2EZc,614
 pyqrack/pauli.py,sha256=TUm1SN_HLz3eMW9gD_eg-IYXcMCyr36mYSytq_ExZIU,673
-pyqrack/qrack_ace_backend.py,sha256=FN4Dqs0yMLy60nLpEhBh0L2LcWjtDRRGhQaVzthU5F0,44721
+pyqrack/qrack_ace_backend.py,sha256=NH8O8i-3_TI73-d2WAz6KzWN1GFlu9iZcCmXtpfJjcQ,50441
 pyqrack/qrack_circuit.py,sha256=QK2dtPYurdXuw-efVq7lYz40_480NxCejXhOVwBGkXs,20122
 pyqrack/qrack_neuron.py,sha256=a3F0hduVvXZi9aXsY0du5hBFXpMq_R5UHHDOv2-YtFM,9305
 pyqrack/qrack_neuron_torch_layer.py,sha256=OhNSldzaqLaMoNBkin68j8QWOOiuZCQJDZPgSDRI2Fk,6330
@@ -9,13 +9,13 @@ pyqrack/qrack_simulator.py,sha256=uETnc4eNOg7w0f9aO_BbRMnQHr-XrA2Doy1kkD7cW9Q,14
 pyqrack/qrack_stabilizer.py,sha256=AJe7dfFcxFKyig3tjWXw0UKhXer5Wl9QNvjNNqlOL5M,2182
 pyqrack/quimb_circuit_type.py,sha256=iC0CCpZBGhziFC8-uBCH43Mi29uvVUrtBG6W9YBlyps,638
 pyqrack/qrack_system/__init__.py,sha256=PUterej-xpA4BqFmiBrQCMeTQlsRf-K8Dxnwp-iVvUQ,343
-pyqrack/qrack_system/qrack_system.py,sha256=57tk3CKihY0x_eYnwcGthcDLIA02YrJU2H3dk7Yw8zw,44230
-pyqrack/qrack_system/qrack_lib/qrack_pinvoke.dll,sha256=n0aHTA9AeVD9wDogofAd3U2qQEyrXzNHmvNLGDX1Qj0,1746432
+pyqrack/qrack_system/qrack_system.py,sha256=WY4gLHfkHrNKoYOkHMqBE4wUxwGgCmAk_GgujGRZms0,44281
+pyqrack/qrack_system/qrack_lib/qrack_pinvoke.dll,sha256=9llAt7eJzXW4F3qsBufp-ZQVl0MbagCKCVP9lG-120c,1746432
 pyqrack/stats/__init__.py,sha256=hI715MGW7D4mDYhUFpRI4ZLsynYDO4tN-rjsuuYbG6Q,282
 pyqrack/stats/load_quantized_data.py,sha256=_1w9BPrZNreP0wOAyaAZHdEGKoGiI7tMeFD9P3eyJC0,1219
 pyqrack/stats/quantize_by_range.py,sha256=0eBIqByIxa4vfm4fQGZLAMGR9y8raxde3e5rvUWJ_dQ,2326
-pyqrack_cpu_complex128-1.58.0.dist-info/LICENSE,sha256=IdAVedmFOPQtHi_XeEI9OhJwUuwlT6tCJwrT55zAn3w,1090
-pyqrack_cpu_complex128-1.58.0.dist-info/METADATA,sha256=pzUwt7WHjHCNTgRFGx1GtZ45x0hnH9VklDo-ePQV3Zc,6005
-pyqrack_cpu_complex128-1.58.0.dist-info/WHEEL,sha256=JMWfR_Dj7ISokcwe0cBhCfK6JKnIi-ZX11L6d_ntt6o,98
-pyqrack_cpu_complex128-1.58.0.dist-info/top_level.txt,sha256=YE_3q9JTGRLMilNg2tGP1y7uU-Dx8PDao2OhwoIbv8E,8
-pyqrack_cpu_complex128-1.58.0.dist-info/RECORD,,
+pyqrack_cpu_complex128-1.58.10.dist-info/LICENSE,sha256=IdAVedmFOPQtHi_XeEI9OhJwUuwlT6tCJwrT55zAn3w,1090
+pyqrack_cpu_complex128-1.58.10.dist-info/METADATA,sha256=psjp_ya6xhbcGMUkX2IdN_PyajF1boAbW268lRlsqvM,6006
+pyqrack_cpu_complex128-1.58.10.dist-info/WHEEL,sha256=JMWfR_Dj7ISokcwe0cBhCfK6JKnIi-ZX11L6d_ntt6o,98
+pyqrack_cpu_complex128-1.58.10.dist-info/top_level.txt,sha256=YE_3q9JTGRLMilNg2tGP1y7uU-Dx8PDao2OhwoIbv8E,8
+pyqrack_cpu_complex128-1.58.10.dist-info/RECORD,,