From 1c9ccdc4c62c8c3037408ad5d4a818c867a8f4e6 Mon Sep 17 00:00:00 2001
From: Jacan Chaplais <jacan.chaplais@gmail.com>
Date: Tue, 7 Nov 2023 15:02:01 +0000
Subject: [PATCH] thrust axis domain |0.5 pi| #165

---
 graphicle/calculate.py | 3 ++-
 1 file changed, 2 insertions(+), 1 deletion(-)

diff --git a/graphicle/calculate.py b/graphicle/calculate.py
index 68e9b54..4ad5433 100644
--- 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(