added the compute discrete geodesics