pyqrack-cpu 1.61.0__py3-none-macosx_13_0_x86_64.whl → 1.61.2__py3-none-macosx_13_0_x86_64.whl

pyqrack/qrack_ace_backend.py

@@ -443,49 +443,42 @@ class QrackAceBackend:
     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
+        return azimuth, inclination
 
     def _get_bloch_angles(self, hq):
-        sim = self.sim[hq[0]]
+        sim = self.sim[hq[0]].clone()
         q = hq[1]
+        sim.separate([q])
 
         # 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
+        return azimuth, inclination
 
     def _rotate_to_bloch(
         self, hq, delta_azimuth, delta_inclination
@@ -576,30 +569,25 @@ class QrackAceBackend:
                         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])
+                a, i = [0, 0, 0, 0, 0], [0, 0, 0, 0, 0]
+                a[0], i[0] = self._get_bloch_angles(hq[0])
+                a[1], i[1] = self._get_bloch_angles(hq[1])
+                a[2], i[2] = self._get_lhv_bloch_angles(hq[2])
+                a[3], i[3] = self._get_bloch_angles(hq[3])
+                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:
+                    a_target += a[x]
+                    i_target += i[x]
+
+                a_target /= 5
+                i_target /= 5
+                for x in range(5):
+                    if x == 2:
+                        self._rotate_lhv_to_bloch(hq[x], a_target - a[x], i_target - i[x])
+                    else:
                         self._rotate_to_bloch(hq[x], a_target - a[x], i_target - i[x])
 
         else:
@@ -622,28 +610,23 @@ class QrackAceBackend:
                         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])
+                a, i = [0, 0, 0], [0, 0, 0]
+                a[0], i[0] = self._get_bloch_angles(hq[0])
+                a[1], i[1] = self._get_bloch_angles(hq[1])
+                a[2], i[2] = self._get_lhv_bloch_angles(hq[2])
 
-                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:
+                    a_target += a[x]
+                    i_target += i[x]
+
+                a_target /= 3
+                i_target /= 3
+                for x in range(3):
+                    if x == 2:
+                        self._rotate_lhv_to_bloch(hq[x], a_target - a[x], i_target - i[x])
+                    else:
                         self._rotate_to_bloch(hq[x], a_target - a[x], i_target - i[x])
 
         if phase:

{pyqrack_cpu-1.61.0.dist-info → pyqrack_cpu-1.61.2.dist-info}/METADATA

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

{pyqrack_cpu-1.61.0.dist-info → pyqrack_cpu-1.61.2.dist-info}/RECORD

@@ -1,7 +1,7 @@
 pyqrack/__init__.py,sha256=3tBwfCCD-zQjQ2g1EUZdggKdn-3b2uSFTbT7LS0YLaU,785
 pyqrack/neuron_activation_fn.py,sha256=fQTTFfsvwcot_43Vopacot47IV2Rxk8pelUyuzwpXPs,593
 pyqrack/pauli.py,sha256=wg500wDOwdIU4lEVJoMmjtbAdmtakZYzLPjdzC2rwUQ,654
-pyqrack/qrack_ace_backend.py,sha256=YNrG-Yte5pq7L81G0GzLyL1uEhuOMNkXq8C8eLxo_FY,48818
+pyqrack/qrack_ace_backend.py,sha256=x1s_I_YijLsKkhkQPWZUUQNHiJom_uTjGFoqz1cgDQQ,48358
 pyqrack/qrack_circuit.py,sha256=vDCKGbcEHJDFUKprjCpWgit8lXFnMrPimKHURD2_Hj4,19538
 pyqrack/qrack_neuron.py,sha256=UiJdjAGB6usjAGHWSosSFCUUeIkhh3MtZbsaxfsIsNw,9043
 pyqrack/qrack_neuron_torch_layer.py,sha256=Bs5BLC2GFevfSpo_jSJ2AZl-hfDRJmzlGN9pFw1CtoQ,6160
@@ -16,8 +16,8 @@ pyqrack/qrack_system/qrack_lib/libqrack_pinvoke.dylib,sha256=Cdb_wa2cG_V5lzswBH4
 pyqrack/stats/__init__.py,sha256=Hla85my2fY_roR9lIjGBVpEG7ySOTMwjWa8D6-kgCnY,276
 pyqrack/stats/load_quantized_data.py,sha256=z12u9F7Nt3P-i44nY1xxvso_klS6WIHS3iqq7R2_lqE,1184
 pyqrack/stats/quantize_by_range.py,sha256=UM0_7jJDdQ7g30cR3UQAxkbzkqrmsy1oUfqg0h11FUY,2270
-pyqrack_cpu-1.61.0.dist-info/LICENSE,sha256=HxB-7SaWTuewAk1nz-3_3FUD6QhgX73kNT_taKVUTq8,1069
-pyqrack_cpu-1.61.0.dist-info/METADATA,sha256=cvRYRJwOGuX9K_HsD55ALKOKt-CCmx6eFmdljhJC4AE,5951
-pyqrack_cpu-1.61.0.dist-info/WHEEL,sha256=nZx8s83OrgdDmpcWX-8FgI0sEAjdQimt4SdYsdcCaC8,107
-pyqrack_cpu-1.61.0.dist-info/top_level.txt,sha256=YE_3q9JTGRLMilNg2tGP1y7uU-Dx8PDao2OhwoIbv8E,8
-pyqrack_cpu-1.61.0.dist-info/RECORD,,
+pyqrack_cpu-1.61.2.dist-info/LICENSE,sha256=HxB-7SaWTuewAk1nz-3_3FUD6QhgX73kNT_taKVUTq8,1069
+pyqrack_cpu-1.61.2.dist-info/METADATA,sha256=BzwhQQ7tAt2oQ4IqtDrLcmbP00EybPVvg690dg7mhIg,5951
+pyqrack_cpu-1.61.2.dist-info/WHEEL,sha256=nZx8s83OrgdDmpcWX-8FgI0sEAjdQimt4SdYsdcCaC8,107
+pyqrack_cpu-1.61.2.dist-info/top_level.txt,sha256=YE_3q9JTGRLMilNg2tGP1y7uU-Dx8PDao2OhwoIbv8E,8
+pyqrack_cpu-1.61.2.dist-info/RECORD,,