This paper presents a single shooting method with an approximate Fréchet derivative for efficiently computing geodesics on the Stiefel manifold.
