This reverts commit
5aba1dbd37afe0078ba09dbf6cf142b559d82d2b.
OK, so we'll make the railway basis have a scale in the (x,y) plane
so that Q is at [1,0,0]. This is used for debugging output, only.
Signed-off-by: Ian Jackson <ijackson@chiark.greenend.org.uk>
railway_cos_theta = cos(railway_theta)
railway_cos_phi = cos(railway_phi)
- railway_inplane_basis_x = unit_v(np.hstack((q_plane[0:2], [0])))
+ railway_inplane_basis_x = np.hstack((q_plane[0:2], [0]))
railway_inplane_basis = np.column_stack((
railway_inplane_basis_x,
np.cross([0,0,1], railway_inplane_basis_x),