pyqrack-cuda 1.57.1__tar.gz → 1.58.5__tar.gz

{pyqrack_cuda-1.57.1 → pyqrack_cuda-1.58.5}/Makefile

@@ -18,7 +18,7 @@ help:
 build-deps:
 ifneq ($(OS),Windows_NT)
 ifeq ($(QRACK_PRESENT),)
-	git clone https://github.com/unitaryfund/qrack.git; cd qrack; git checkout f0738286f49fb0ffaeb46d6e5d2294dbfc7f6138; cd ..
+	git clone https://github.com/unitaryfund/qrack.git; cd qrack; git checkout 8c3f33832d334196836b10afe6765bb44326802d; cd ..
 endif
 	mkdir -p qrack/build
 ifeq ($(UNAME_S),Linux)

{pyqrack_cuda-1.57.1/pyqrack_cuda.egg-info → pyqrack_cuda-1.58.5}/PKG-INFO

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

{pyqrack_cuda-1.57.1 → pyqrack_cuda-1.58.5}/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,6 +443,88 @@ class QrackAceBackend:
 
         return qb, lhv
 
+    def get_lhv_bloch_angles(self, q):
+        z = (1 - q.bloch[2]) / 2
+        prob = z**2
+        x = (1 - q.bloch[0]) / 2
+        prob += x**2
+        y = (1 - q.bloch[1]) / 2
+        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_lhv_to_bloch(
+        self, q, azimuth_curr, inclination_curr, azimuth_target, inclination_target
+    ):
+        delta_azimuth = azimuth_target - azimuth_curr
+        delta_inclination = inclination_target - inclination_curr
+
+        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)
+
+        q.mtrx([m00, m01, m10, m11])
+
+    def rotate_to_bloch(
+        self, hq, azimuth_curr, inclination_curr, azimuth_target, inclination_target
+    ):
+        sim = self.sim[hq[0]]
+        q = hq[1]
+
+        delta_azimuth = azimuth_target - azimuth_curr
+        delta_inclination = inclination_target - inclination_curr
+
+        # 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 _correct(self, lq, phase=False):
         hq = self._unpack(lq)
 
@@ -454,6 +575,28 @@ class QrackAceBackend:
                         hq[q].x()
                     else:
                         self.sim[hq[q][0]].x(hq[q][1])
+
+            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[2], i[2] = self.get_bloch_angles(hq[3])
+            p[4], a[3], i[3] = self.get_bloch_angles(hq[4])
+
+            indices = []
+            a_target = 0
+            i_target = 0
+            for x in range(4):
+                if p[x] < 0.5:
+                    continue
+                indices.append(x)
+                a_target += a[x]
+                i_target += i[x]
+
+            if len(indices) > 1:
+                a_target /= len(indices)
+                i_target /= len(indices)
+                for x in indices:
+                    self.rotate_to_bloch(hq[x], a[x], i[x], a_target, i_target)
         else:
             # RMS
             p = [
@@ -473,6 +616,26 @@ class QrackAceBackend:
                     else:
                         self.sim[hq[q][0]].x(hq[q][1])
 
+            p, a, i = [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
+            for x in range(2):
+                if p[x] < 0.5:
+                    continue
+                indices.append(x)
+                a_target += a[x]
+                i_target += i[x]
+
+            if len(indices) > 1:
+                a_target /= len(indices)
+                i_target /= len(indices)
+                for x in indices:
+                    self.rotate_to_bloch(hq[x], a[x], i[x], a_target, i_target)
+
         if phase:
             for q in qb:
                 b = hq[q]

{pyqrack_cuda-1.57.1 → pyqrack_cuda-1.58.5}/pyqrack/qrack_simulator.py

@@ -122,6 +122,9 @@ class QrackSimulator:
     def _ulonglong_byref(self, a):
         return (ctypes.c_ulonglong * len(a))(*a)
 
+    def _longlong_byref(self, a):
+        return (ctypes.c_longlong * len(a))(*a)
+
     def _double_byref(self, a):
         return (ctypes.c_double * len(a))(*a)
 
@@ -190,6 +193,11 @@ class QrackSimulator:
         Qrack.qrack_lib.set_device(self.sid, d)
         self._throw_if_error()
 
+    def set_device_list(self, d):
+        """Set the GPU device ID"""
+        Qrack.qrack_lib.set_device_list(self.sid, len(d), self._longlong_byref(d))
+        self._throw_if_error()
+
     def clone(self):
         return QrackSimulator(cloneSid=self.sid)
 

{pyqrack_cuda-1.57.1 → pyqrack_cuda-1.58.5}/pyqrack/qrack_system/qrack_system.py

@@ -144,7 +144,14 @@ class QrackSystem:
         self.qrack_lib.set_concurrency.argtypes = [c_ulonglong, c_ulonglong]
 
         self.qrack_lib.set_device.restype = None
-        self.qrack_lib.set_device.argtypes = [c_ulonglong, c_ulonglong]
+        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),
+        ]
 
         # pseudo-quantum
 

{pyqrack_cuda-1.57.1 → pyqrack_cuda-1.58.5/pyqrack_cuda.egg-info}/PKG-INFO

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

{pyqrack_cuda-1.57.1 → pyqrack_cuda-1.58.5}/setup.py

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