From: Ian Jackson Date: Sun, 17 Dec 2017 11:34:37 +0000 (+0000) Subject: helixish: debugging output X-Git-Url: https://www.chiark.greenend.org.uk/ucgi/~ian/git?a=commitdiff_plain;h=995951970a5656d0f53efb7a3508fe22320e7862;p=moebius3.git helixish: debugging output Signed-off-by: Ian Jackson --- diff --git a/helixish.py b/helixish.py index d43c622..39d525a 100644 --- 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)),