ncp3 = np.linalg.norm(cp[3])
cpt = ncp3 * ncp3 / 4
cp2scale = m._cp2b.point_at_t(ncp3/2)
- cp1scale = 0.75 * cpt + 0.33 * (1-cpt)
+ cp1scale = 1.0 * cpt + 0.33 * (1-cpt)
#print('u=%d ncp3=%f cp2scale=%f' % (u, ncp3, cp2scale),
# file=sys.stderr)
cp[1] = cp[0] + cp1scale * m.edge .dirn (theta)