From: Ian Jackson <ijackson@chiark.greenend.org.uk>
Date: Wed, 29 Nov 2017 17:51:15 +0000 (+0000)
Subject: helixish: better debug
X-Git-Url: https://www.chiark.greenend.org.uk/ucgi/~ian/git?a=commitdiff_plain;h=3c403f828eff980a2425952fa5addc681ee5b220;p=moebius3.git

helixish: better debug

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

diff --git a/helixish.py b/helixish.py
index 14df461..19c592b 100644
--- a/helixish.py
+++ b/helixish.py
@@ -73,16 +73,20 @@ class HelixishCurve():
                                       dp,
                                       dPQplane_normal,
                                       p));
-    dbg(dPQplane_basis)
-    dPQplane_basis = np.vstack((dPQplane_basis, [0,0,0,1]))
     #dbg(dPQplane_basis)
+    dPQplane_basis = np.vstack((dPQplane_basis, [0,0,0,1]))
+    dbg(dPQplane_basis)
     dPQplane_into = np.linalg.inv(dPQplane_basis)
+    dbg(dPQplane_into)
 
+    p_plane_check = augmatmultiply(dPQplane_into, p)
     dp_plane = augmatmultiply(dPQplane_into, dp, augwith=0)
     dq_plane = augmatmultiply(dPQplane_into, dq, augwith=0)
     q_plane  = augmatmultiply(dPQplane_into, q)
     dist_pq_plane = np.linalg.norm(q_plane)
 
+    dbg('plane:', p_plane_check, dp_plane, dq_plane, q_plane)
+
     # two circular arcs of equal maximum possible radius
     # algorithm courtesy of Simon Tatham (`Railway problem',
     # pers.comm. to ijackson@chiark 23.1.2004)