pyqrack 1.60.0__py3-none-macosx_15_0_arm64.whl → 1.61.1__py3-none-macosx_15_0_arm64.whl

pyqrack/qrack_ace_backend.py

@@ -332,7 +332,8 @@ class QrackAceBackend:
 
             # You can still "monkey-patch" this, after the constructor.
             if "QRACK_QUNIT_SEPARABILITY_THRESHOLD" not in os.environ:
-                self.sim[i].set_sdrp(0.02375)
+                # (1 - 1 / sqrt(2)) / 4 (but empirically tuned)
+                self.sim[i].set_sdrp(0.073223304703363119)
 
     def clone(self):
         return QrackAceBackend(toClone=self)
@@ -442,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
@@ -575,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:
@@ -621,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 /= 5
+                i_target /= 5
+                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-1.60.0.dist-info → pyqrack-1.61.1.dist-info}/METADATA

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

{pyqrack-1.60.0.dist-info → pyqrack-1.61.1.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=_M03YHvczhw65gtg0ly8qlen7aOPDLC1K5gtM-cTgbI,48741
+pyqrack/qrack_ace_backend.py,sha256=2O5JxL76w-QzrF1-D7yY1sa1MXjw6VKvATR9s3DZoj4,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=V2_marSne0NGSonHzK5
 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-1.60.0.dist-info/LICENSE,sha256=HxB-7SaWTuewAk1nz-3_3FUD6QhgX73kNT_taKVUTq8,1069
-pyqrack-1.60.0.dist-info/METADATA,sha256=PVCNBbl44Lun0SlVUj1KJkFQRzDEkcVFD-rtR6ye4yI,5892
-pyqrack-1.60.0.dist-info/WHEEL,sha256=uKwoxROtjcnx4Gbh9Fbsan8Fm8wY5xdQpS0Bc4agBqM,106
-pyqrack-1.60.0.dist-info/top_level.txt,sha256=YE_3q9JTGRLMilNg2tGP1y7uU-Dx8PDao2OhwoIbv8E,8
-pyqrack-1.60.0.dist-info/RECORD,,
+pyqrack-1.61.1.dist-info/LICENSE,sha256=HxB-7SaWTuewAk1nz-3_3FUD6QhgX73kNT_taKVUTq8,1069
+pyqrack-1.61.1.dist-info/METADATA,sha256=oOcJTndhIrtyj6SIXXSgyO6jU_ETUcyk-Emv0h6rtSY,5892
+pyqrack-1.61.1.dist-info/WHEEL,sha256=uKwoxROtjcnx4Gbh9Fbsan8Fm8wY5xdQpS0Bc4agBqM,106
+pyqrack-1.61.1.dist-info/top_level.txt,sha256=YE_3q9JTGRLMilNg2tGP1y7uU-Dx8PDao2OhwoIbv8E,8
+pyqrack-1.61.1.dist-info/RECORD,,