A Riemannian approach for computing geodesies in elastic shape analysis

We propose a Riemannian quasi-Newton method to compute a geodesic invariant to scaling, translation, rotation and reparameterization and show that it is more efficient than the current state-of-the-art coordinate-descent/path-straightening approach.