thrust axis domain |0.5 pi| #165

jacanchaplais · Nov 7, 2023 · 1c9ccdc · 1c9ccdc
1 parent 1ad61fe
commit 1c9ccdc
Showing 1 changed file with 2 additions and 1 deletion.
diff --git a/graphicle/calculate.py b/graphicle/calculate.py
@@ -806,7 +806,8 @@ def thrust(
     axis : ndarray[float64], optional
         The x, y, and z components of the thrust axis, respectively.
     """
-    domain = (0.0, 0.5 * math.pi)
+    half_pi = 0.5 * math.pi
+    domain = (-half_pi, half_pi)
     rng = np.random.default_rng(seed=rng_seed)
     guess = rng.uniform(*domain, size=2)
     optim = spo.minimize(