labs.py

Author: WrathfulSpatula

scripts/labs.py

@@ -0,0 +1,45 @@
+# Spin Glass Ground State (considered NP-complete)
+# Produced by Dan Strano, Elara (the OpenAI custom GPT)
+
+from pyqrackising import spin_glass_solver
+import networkx as nx
+import numpy as np
+
+import sys
+
+
+# LABS spin glass
+def generate_labs_maxcut(N, seed=None):
+    if seed is not None:
+        np.random.seed(seed)
+
+    # Initialize random ±1 sequence
+    s = np.random.choice([-1, 1], size=N)
+
+    # Initialize weight matrix
+    W = np.zeros((N, N))
+
+    # Build LABS autocorrelation-derived weights
+    for k in range(1, N):
+        for i in range(N - k):
+            j = i + k
+            # Each off-peak correlation contributes to (i,j)
+            W[i, j] += 1
+            W[j, i] += 1  # symmetry
+
+    # Normalize or scale weights (optional)
+    W /= np.max(W)
+
+    return W, s
+
+
+if __name__ == "__main__":
+    n_nodes = int(sys.argv[1]) if len(sys.argv) > 1 else 64
+    quality = int(sys.argv[2]) if len(sys.argv) > 2 else None
+    seed = int(sys.argv[3]) if len(sys.argv) > 3 else None
+    G_m, s = generate_labs_maxcut(n_nodes, seed)
+    best_bitstring, best_cut_value, best_cut, best_energy = spin_glass_solver(G_m, quality=quality)
+
+    print(f"Ground State Energy: {best_energy}")
+    if not (seed is None):
+        print(f"Seed {seed} solution: {best_bitstring}")

scripts/maxcut_approx.py

@@ -77,8 +77,8 @@ def generate_spin_glass_graph(n_nodes=16, degree=3, seed=None):
     # NP-complete spin glass
     # G = generate_spin_glass_graph(seed=42)
 
-    # bitstring, cut_value, cut = maxcut_tfim_sparse(G)
-    bitstring, cut_value, cut, energy = spin_glass_solver_sparse(G, is_spin_glass=False)
+    bitstring, cut_value, cut = maxcut_tfim_sparse(G)
+    # bitstring, cut_value, cut, energy = spin_glass_solver_sparse(G, is_spin_glass=False)
     # bitstring, cut_value, cut, energy = spin_glass_solver_sparse(G, best_guess=maxcut_tfim_sparse(G, quality=8)[0])
     print((bitstring, cut_value, cut))
 

scripts/spin_glass.py

@@ -1,7 +1,7 @@
 # Spin Glass Ground State (considered NP-complete)
 # Produced by Dan Strano, Elara (the OpenAI custom GPT)
 
-from pyqrackising import spin_glass_solver
+from pyqrackising import spin_glass_solver_sparse
 import networkx as nx
 import numpy as np
 
@@ -19,7 +19,7 @@ def generate_spin_glass_graph(n_nodes=64, degree=3, seed=None):
 if __name__ == "__main__":
     # NP-complete spin glass
     G = generate_spin_glass_graph(n_nodes=64, seed=42)
-    best_bitstring, best_cut_value, best_cut, best_energy = spin_glass_solver(G)
+    best_bitstring, best_cut_value, best_cut, best_energy = spin_glass_solver_sparse(G)
 
     print((best_bitstring, best_cut_value, best_cut))
     print(f"Ground State Energy: {best_energy}")