helixish: debugging output

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

diff --git a/helixish.py b/helixish.py
index d43c622d9015d7b96b9c6082317d2a28071ea525..39d525a4e9e19cf8652b49af0f66aff6fe625ee8 100644 (file)
--- a/helixish.py
+++ b/helixish.py
@@ -61,8 +61,9 @@ class HelixishCurve():
     dPQplane_basis = np.vstack((dPQplane_basis, [0,0,0,1]))
     dbg(dPQplane_basis)
 
-    #for ax in range(0,3):
-    #  vdbg().arrow(dPQplane_basis[0:3,3], dPQplane_basis[0:3,ax])
+    for ax in range(0,3):
+      vdbg().arrow(dPQplane_basis[0:3,3], dPQplane_basis[0:3,ax],
+                   colour= [ c * (1.0 - ax*0.2) for c in (1,1,1) ])
 
     dPQplane_into = np.linalg.inv(dPQplane_basis)
     dbg(dPQplane_into)
@@ -87,6 +88,7 @@ class HelixishCurve():
     railway_phi   = atan2(*dq_plane[0:2]) - railway_angleoffset
     railway_cos_theta = cos(railway_theta)
     railway_cos_phi   = cos(railway_phi)
+    dbg('railway:', railway_theta, railway_phi, railway_angleoffset)
     if railway_cos_theta**2 + railway_cos_phi**2 > 1E-6:
       railway_roots = np.roots([
         2 * (1 + cos(railway_theta - railway_phi)),