From 3c403f828eff980a2425952fa5addc681ee5b220 Mon Sep 17 00:00:00 2001 From: Ian Jackson Date: Wed, 29 Nov 2017 17:51:15 +0000 Subject: [PATCH] helixish: better debug Signed-off-by: Ian Jackson --- helixish.py | 8 ++++++-- 1 file changed, 6 insertions(+), 2 deletions(-) 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) -- 2.30.2