From 86b403b30fbe40ec4acd46def6c2bbef2d32e10b Mon Sep 17 00:00:00 2001
From: Noa Aarts <noa@voorwaarts.nl>
Date: Tue, 27 Jan 2026 19:15:40 +0100
Subject: [PATCH] add hyperspace sampling

---
 src/ga_qas.py | 48 ++++++++++++++++++++++++++++++++++++++++++++----
 1 file changed, 44 insertions(+), 4 deletions(-)

diff --git a/src/ga_qas.py b/src/ga_qas.py
index f933e5d..92e2a99 100755
--- a/src/ga_qas.py
+++ b/src/ga_qas.py
@@ -28,12 +28,52 @@ gate_set: list[GateType] = [
 ]
 
 
-def sample_hyperspace(*args: tuple[int, int] | tuple[float, float]):
-    minimums: list[float | int] = [arg[0] for arg in args]
-    maximums: list[float | int] = [arg[1] for arg in args]
+def sample_hyperspace(
+    *args: tuple[int, int] | tuple[float, float], seed: int = 2010392991
+):
+    minimums: tuple[float | int, ...] = tuple(arg[0] for arg in args)
+    maximums: tuple[float | int, ...] = tuple(arg[1] for arg in args)
+    diffs: tuple[float | int, ...] = tuple(
+        ma - mi for mi, ma in zip(minimums, maximums)
+    )
+    idiffs: tuple[float, ...] = tuple(1.0 / diff for diff in diffs)
+    rng: random.Random = random.Random(seed)
+
+    previous_points: set[tuple[float | int, ...]] = set()
+
+    def dist(point: tuple[float | int, ...], other: tuple[float | int, ...]) -> float:
+        return sum(((p - o) * idiff) ** 2 for p, o, idiff in zip(point, other, idiffs))
 
     while True:
-        yield ()
+        sampled_points: list[tuple[float | int, ...]] = [
+            tuple(
+                min
+                + (
+                    rng.randint(min, min + diff)
+                    if isinstance(min, int) and isinstance(diff, int)
+                    else rng.uniform(min, min + diff)
+                )
+                for min, diff in zip(minimums, diffs)
+            )
+            for _ in range(10)
+        ]
+
+        if len(previous_points) == 0:
+            previous_points.add(sampled_points[0])
+            yield sampled_points[0]
+
+        min_distances: list[float] = [
+            min((dist(point, other) for other in previous_points))
+            for point in sampled_points
+        ]
+
+        mdist: float = max(min_distances)
+        point: tuple[float | int, ...] = sampled_points[
+            [i for i, j in enumerate(min_distances) if j == mdist][0]
+        ]
+
+        previous_points.add(point)
+        yield point
 
 
 def plot_best_circuits(best_circuits: list[QuantumCircuit]) -> None: