pyqrack-cuda 1.58.7__tar.gz → 1.58.9__tar.gz

{pyqrack_cuda-1.58.7/pyqrack_cuda.egg-info → pyqrack_cuda-1.58.9}/PKG-INFO

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

{pyqrack_cuda-1.58.7 → pyqrack_cuda-1.58.9}/pyqrack/qrack_ace_backend.py

@@ -443,6 +443,20 @@ 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]
@@ -471,6 +485,24 @@ class QrackAceBackend:
 
         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
     ):
@@ -493,17 +525,6 @@ class QrackAceBackend:
 
         sim.mtrx([m00, m01, m10, m11], q)
 
-    def logit(self, x):
-        # Theoretically, these limit points are "infinite,"
-        # but precision caps out between 36 and 37:
-        if 5e-17 > (1 - x):
-            return 37
-        # For the negative limit, the precision caps out
-        # between -37 and -38
-        elif x < 1e-17:
-            return -38
-        return max(-38, min(37, math.log(x / (1 - x))))
-
 
     def _correct(self, lq, phase=False):
         hq = self._unpack(lq)
@@ -559,6 +580,7 @@ class QrackAceBackend:
             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[2], a[2], i[2] = self.get_lhv_bloch_angles(hq[2])
             p[3], a[3], i[3] = self.get_bloch_angles(hq[3])
             p[4], a[4], i[4] = self.get_bloch_angles(hq[4])
 
@@ -570,7 +592,7 @@ class QrackAceBackend:
                 if p[x] < 0.5:
                     continue
                 indices.append(x)
-                w = self.logit(1 - p[x])
+                w = (1.5 - p[x])
                 a_target += w * a[x]
                 i_target += w * i[x]
                 weight += w
@@ -579,7 +601,10 @@ class QrackAceBackend:
                 a_target /= weight
                 i_target /= weight
                 for x in indices:
-                    self.rotate_to_bloch(hq[x], a[x], i[x], a_target, i_target)
+                    if x == 2:
+                        self.rotate_lhv_to_bloch(hq[x], a[x], i[x], a_target, i_target)
+                    else:
+                        self.rotate_to_bloch(hq[x], a[x], i[x], a_target, i_target)
         else:
             # RMS
             p = [
@@ -599,9 +624,10 @@ class QrackAceBackend:
                     else:
                         self.sim[hq[q][0]].x(hq[q][1])
 
-            p, a, i = [0, 0], [0, 0], [0, 0]
+            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])
+            p[2], a[2], i[2] = self.get_lhv_bloch_angles(hq[2])
 
             indices = []
             a_target = 0
@@ -611,8 +637,8 @@ class QrackAceBackend:
                 if p[x] < 0.5:
                     continue
                 indices.append(x)
-                w = self.logit(1.0 - p[x])
-                a_target += x * a[x]
+                w = (1.5 - p[x])
+                a_target += w * a[x]
                 i_target += w * i[x]
                 weight += w
 
@@ -620,7 +646,10 @@ class QrackAceBackend:
                 a_target /= weight
                 i_target /= weight
                 for x in indices:
-                    self.rotate_to_bloch(hq[x], a[x], i[x], a_target, i_target)
+                    if x == 2:
+                        self.rotate_lhv_to_bloch(hq[x], a[x], i[x], a_target, i_target)
+                    else:
+                        self.rotate_to_bloch(hq[x], a[x], i[x], a_target, i_target)
 
         if phase:
             for q in qb:

{pyqrack_cuda-1.58.7 → pyqrack_cuda-1.58.9/pyqrack_cuda.egg-info}/PKG-INFO

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

{pyqrack_cuda-1.58.7 → pyqrack_cuda-1.58.9}/setup.py

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