helixish: debugging output

Signed-off-by: Ian Jackson <ijackson@chiark.greenend.org.uk>

diff --git a/helixish.py b/helixish.py
index 1f48fa891034cf0b8f8da57400e05e997c1ccc66..17922d6e3affe25955f0760830867b97f8cd37f9 100644 (file)
--- a/helixish.py
+++ b/helixish.py
@@ -114,11 +114,13 @@ class HelixishCurve():
     vdbg_railway_plane_angle([1, 0, 0.1], railway_phi,   color=(1, 0.5, 0))
 
     if railway_cos_theta**2 + railway_cos_phi**2 > 1E-6:
-      railway_roots = np.roots([
+      railway_polynomial = [
         2 * (1 + cos(railway_theta - railway_phi)),
         2 * (railway_cos_theta - railway_cos_phi),
         -1
-        ])
+        ]
+      railway_roots = np.roots(railway_polynomial)
+      dbg('railway poly, roots:', railway_polynomial, railway_roots)
       for railway_r in railway_roots:
         dbg(' twoarcs root r=',railway_r)