A: (-8, -13, -9)
B: (12, 7, 1)
u = Vector((2, 2, 1))
g = IntersectPath((x - 12) / 2 + (0 * z) = (y - 7) / 2, (x - 12) / 2 = (z - 1) / 1)