In this paper, a curve evolution approach for the computation of geodesic curves on 3D surfaces is presented. The algorithm is based on deforming, via the curve shortening flow, an arbitrary initial curve ending at two given surface points. The 3D curve shortening flow is first transformed into an equivalent 2D one. This 2D flow is implemented, using an efficient numerical algorithm for curve evolution with fixed end points. © 1995.