pyqrackising 9.5.1__tar.gz → 9.5.3__tar.gz

{pyqrackising-9.5.1/pyqrackising.egg-info → pyqrackising-9.5.3}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: pyqrackising
-Version: 9.5.1
+Version: 9.5.3
 Summary: Fast MAXCUT, TSP, and sampling heuristics from near-ideal transverse field Ising model (TFIM)
 Home-page: https://github.com/vm6502q/PyQrackIsing
 Author: Dan Strano

{pyqrackising-9.5.1 → pyqrackising-9.5.3}/pyproject.toml

@@ -10,7 +10,7 @@ build-backend = "setuptools.build_meta"
 
 [project]
 name = "pyqrackising"
-version = "9.5.1"
+version = "9.5.3"
 requires-python = ">=3.8"
 description = "Fast MAXCUT, TSP, and sampling heuristics from near-ideal transverse field Ising model (TFIM)"
 readme = {file = "README.txt", content-type = "text/markdown"}

{pyqrackising-9.5.1 → pyqrackising-9.5.3}/pyqrackising/otoc.py

@@ -43,27 +43,26 @@ def get_otoc_hamming_distribution(J=-1.0, h=2.0, z=4, theta=0.0, t=5, n_qubits=6
         rev = probability_by_hamming_weight(h, J, z, phi + np.pi, t, n_qubits + 1)
         diff_phi = rev - fwd
 
-        # Lambda (Y-axis) is at a right angle to both J and h,
-        # so there is no difference in this dimension.
+        diff_lam = diff_theta + diff_phi
+
+        diff_theta[0] += 1.0
+        diff_phi[0] += 1.0
+        diff_lam[0] += 1.0
+
+        diff_theta += max_entropy
+        diff_phi += max_entropy
+        diff_lam += max_entropy
 
-        diff_z[0] += n_qubits
-        entropy_frac = 0
         for b in pauli_string:
             match b:
                 case 'X':
                     diff_z += diff_theta
-                    entropy_frac += 1
                 case 'Z':
                     diff_z += diff_phi
-                    entropy_frac += 1
                 case 'Y':
-                    diff_z += diff_theta + diff_phi
-                    entropy_frac += 1
+                    diff_z += diff_lam
                 case _:
-                    pass
-
-        entropy_frac /= n_qubits
-        diff_z = ((1 - entropy_frac) / n_qubits) * diff_z + entropy_frac * max_entropy
+                    diff_z[0] += 1.0
 
     # Normalize:
     diff_z /= diff_z.sum()
@@ -185,6 +184,12 @@ def get_inv_dist(butterfly_idx_x, butterfly_idx_z, n_qubits, row_len, col_len):
 def generate_otoc_samples(J=-1.0, h=2.0, z=4, theta=0.0, t=5, n_qubits=65, pauli_strings = ['X' + 'I' * 64], shots=100, is_orbifold=True):
     thresholds = fix_cdf(get_otoc_hamming_distribution(J, h, z, theta, t, n_qubits, pauli_strings))
 
+    entropy_frac = 0.0
+    for pauli_string in pauli_strings:
+        pauli_string = list(pauli_string)
+        entropy_frac += pauli_string.count('X') + pauli_string.count('Y') + pauli_string.count('Z')
+    entropy_frac /= n_qubits * len(pauli_strings)
+
     row_len, col_len = factor_width(n_qubits)
     inv_dist = np.zeros(n_qubits, dtype=np.float64)
     for pauli_string in pauli_strings:
@@ -196,6 +201,7 @@ def generate_otoc_samples(J=-1.0, h=2.0, z=4, theta=0.0, t=5, n_qubits=65, pauli
             inv_dist += get_willow_inv_dist(butterfly_idx_x, butterfly_idx_z, n_qubits, row_len, col_len)
         inv_dist /= 2.0
 
+    qubit_pows = [1 << q for q in range(n_qubits)]
     samples = []
     for _ in range(shots):
         # First dimension: Hamming weight
@@ -208,6 +214,13 @@ def generate_otoc_samples(J=-1.0, h=2.0, z=4, theta=0.0, t=5, n_qubits=65, pauli
             continue
 
         # Second dimension: permutation within Hamming weight
-        samples.append(take_sample(n_qubits, 0, m, inv_dist))
+        if np.random.random() < entropy_frac:
+            bit_pows = np.random.choice(qubit_pows, size=m, replace=False)
+            sample = 0
+            for bit_pow in bit_pows:
+                sample |= bit_pow
+            samples.append(sample)
+        else:
+            samples.append(take_sample(n_qubits, 0, m, inv_dist))
 
     return samples

{pyqrackising-9.5.1 → pyqrackising-9.5.3/pyqrackising.egg-info}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: pyqrackising
-Version: 9.5.1
+Version: 9.5.3
 Summary: Fast MAXCUT, TSP, and sampling heuristics from near-ideal transverse field Ising model (TFIM)
 Home-page: https://github.com/vm6502q/PyQrackIsing
 Author: Dan Strano

{pyqrackising-9.5.1 → pyqrackising-9.5.3}/setup.py

@@ -7,7 +7,7 @@ with open(README_PATH) as readme_file:
 
 setup(
     name='pyqrackising',
-    version='9.5.1',
+    version='9.5.3',
     author='Dan Strano',
     author_email='stranoj@gmail.com',
     description='Fast MAXCUT, TSP, and sampling heuristics from near-ideal transverse field Ising model (TFIM)',