+ -1,
+ ]
+ railway_roots = np.roots(railway_polynomial)
+ dbg('railway poly, roots:', railway_polynomial, railway_roots)
+
+ #vdbg_railway.circle([0,0,0], [0,0, dist_pq_plane], color=(.5,0,0))
+ #vdbg_railway.circle([1,0,0], [0,0, 0.05], color=(.5,0,0))
+ #vdbg().circle(p, dPQplane_normal * dist_pq_plane, color=(.5,.5,0))
+
+ for railway_r_pq1 in railway_roots:
+ # roots for r are calculated based on coordinates where
+ # Q is at (1,0) but our PQ distance is different
+ railway_r = railway_r_pq1 * dist_pq_plane
+ dbg(' twoarcs root r_pq1=', railway_r_pq1, 'r=',railway_r,
+ railway_polynomial[0] * railway_r_pq1 * railway_r_pq1 +
+ railway_polynomial[1] * railway_r_pq1 +
+ railway_polynomial[2]
+ )
+
+ vdbg_railway.circle([0,0,0], [0,0, railway_r], color=(1,0,0))
+ #vdbg().circle(p, dPQplane_normal * railway_r, color=(1,1,0))