Geodesics measure the shortest distance (either locally or globally) between two points on a curved surface and serve as a fundamental tool in digital geometry processing. Suppose that we have a parameterized path γ(t)=x(u(t),v(t)) on a surface x=x(u,v) with γ(0)=p and γ(1)=q. We formulate the two-point geodesic problem into a minimization problem [Formula: see text], where H(s) satisfies and H''(s) ≥ 0 for . In our implementation, we choose H(s)=es2-1 and show that it has several unique advantages over other choices such as H(s)=s2 and H(s)=s. It is also a minimizer of the traditional geodesic length variational and able to guarantee the uniqueness and regularity in terms of curve parameterization. In the discrete setting, we construct the initial path by a sequence of moveable points {xi}i=1n and minimize ∑i=1n H(||xi - xi+1||). The resulting points are evenly spaced along the path. It's obvious that our algorithm can deal with parametric surfaces. Considering that meshes, point clouds and implicit surfaces can be transformed into a signed distance function (SDF), we also discuss its implementation on a general SDF. Finally, we show that our method can be extended to solve a general least-cost path problem. We validate the proposed algorithm in terms of accuracy, performance and scalability, and demonstrate the advantages by extensive comparisons.