Numerical Accuracy of Ladder Schemes for Parallel Transport on Manifolds

by   Nicolas Guigui, et al.

Parallel transport is a fundamental tool to perform statistics on Rie-mannian manifolds. Since closed formulae don't exist in general, practitioners often have to resort to numerical schemes. Ladder methods are a popular class of algorithms that rely on iterative constructions of geodesic parallelograms. And yet, the literature lacks a clear analysis of their convergence performance. In this work, we give Taylor approximations of the elementary constructions of Schild's ladder and the pole ladder with respect to the Riemann curvature of the underlying space. We then prove that these methods can be iterated to converge with quadratic speed, even when geodesics are approximated by numerical schemes. We also contribute a new link between Schild's ladder and the Fanning Scheme which explains why the latter naturally converges only linearly. The extra computational cost of ladder methods is thus easily compensated by a drastic reduction of the number of steps needed to achieve the requested accuracy. Illustrations on the 2-sphere, the space of symmetric positive definite matrices and the special Euclidean group show that the theoretical errors we have established are measured with a high accuracy in practice. The special Euclidean group with an anisotropic left-invariant metric is of particular interest as it is a tractable example of a non-symmetric space in general , which reduces to a Riemannian symmetric space in a particular case. As a secondary contribution, we compute the covariant derivative of the curvature in this space.


page 1

page 2

page 3

page 4


Riemannian Geometry of Symmetric Positive Definite Matrices via Cholesky Decomposition

We present a new Riemannian metric, termed Log-Cholesky metric, on the m...

Parallel Transport with Pole Ladder: a Third Order Scheme in Affine Connection Spaces which is Exact in Affine Symmetric Spaces

Parallel transport is an important step in many discrete algorithms for ...

Additive Models for Symmetric Positive-Definite Matrices, Riemannian Manifolds and Lie groups

In this paper an additive regression model for a symmetric positive-defi...

Numerical approximation of boundary value problems for curvature flow and elastic flow in Riemannian manifolds

We present variational approximations of boundary value problems for cur...

Symmetric Algorithmic Components for Shape Analysis with Diffeomorphisms

In computational anatomy, the statistical analysis of temporal deformati...

Principal Autoparallel Analysis: Data Analysis in Weitzenböck Space

The statistical analysis of data lying on a differentiable, locally Eucl...

Geometric subdivision and multiscale transforms

Any procedure applied to data, and any quantity derived from data, is re...