ef_dirn = unit_v(cp[cp1i] - cp[cp0i])
ef_len = np.linalg.norm(np.array(fc_input[fi]) - np.array(fc_input[ei]))
f = e + ef_dirn * ef_len
+ dbg(repr((end, e,f, ef_dirn, ef_len)))
fc_input[ei] = e
fc_input[fi] = f
+ dbg(repr(fc_input))
+
findcurve_epsilon = 0.01
try: