collision returns VectorD with (sin(arcsin(y1y2.x/y1y2.abs), cos(arccos(y1y2.y/y1y2.abs))