 # 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 statistical computing on manifolds. Numerical methods based on Jacobi fields or geodesics parallelograms are currently used in geometric data processing. In this last class, pole ladder is a simplification of Schild's ladder for the parallel transport along geodesics that was shown to be particularly simple and numerically stable in Lie groups. So far, these methods were shown to be first order approximations of the Riemannian parallel transport, but higher order error terms are difficult to establish. In this paper, we build on a BCH-type formula on affine connection spaces to establish the behavior of one pole ladder step up to order 5. It is remarkable that the scheme is of order three in general affine connection spaces with a symmetric connection, much higher than expected. Moreover, the fourth-order term involves the covariant derivative of the curvature only, which is vanishing in locally symmetric space. We show that pole ladder is actually locally exact in these spaces, and even almost surely globally exact in Riemannian symmetric manifolds. These properties make pole ladder a very attractive alternative to other methods in general affine manifolds with a symmetric connection.

## Authors

##### This week in AI

Get the week's most popular data science and artificial intelligence research sent straight to your inbox every Saturday.

## 1 Pole ladder in an affine connection spaces

This section establishes a Taylor expansion of the approximation provided by one step of pole ladder. We first detail the pole ladder algorithm. We then turn to an equivalent on affine manifolds of the BCH formula for Lie group before establishing the approximation order of pole ladder. Figure 2: Schematic of the pole ladder procedure to parallel transports the vector u0=logp0(x0) along the geodesic by arc curve (p0,p1,…pn). Left: First rug of the ladder using an approximate geodesic parallelogram. Right: The method is iterated with rungs at each point sampled along the curve.

Pole ladder is a modification of Schild’s ladder for the parallel transport along geodesic curves which is based on the observation that this geodesic can be taken as one of the diagonals of the geodesic parallelogram. The general procedure is described in Fig.2: the vector (rescaled if needed to be small enough) is identified to the geodesic segment for . One compute the mid-point between and and then extend twice the geodesic from to to obtain the point . The vector is the pole ladder approximation of the parallel transport of at . An alternative is to first compute and then extend twice the geodesic from to to obtain the point . The vector is another approximation of the parallel transport of at .

### 1.2 Algorithms for one step of the ladder

Pole ladder assumes that the successive points of the geodesic curve are in a sufficiently small convex normal neighborhood so that this mid-point is unique. In order to analyze one single step of the ladder along a geodesic segment, we simplify the notations in the sequel according to figure 3. The original pole ladder procedure is summarized in Algorithm 1.

Instead of encoding a geodesic by its initial point and initial tangent vector , we can also encode a geodesic segment by its end-points . We can also rephrase the doubling of the geodesics with a mid-point formulation. This leads to an alternative version of pole ladder using mid-point geodesic symmetries (Algorithm 2). Although the two versions of pole ladder are theoretically completely equivalent, it was experimentally observed that this second version was numerically much more stable when applied to diffeomorphisms acting on images [JDM18]. This is most probably be due to the more symmetric formulation that minimizes numerical errors.

Notice that the mid-point of a geodesic segment is well defined even in affine connection spaces: this is the point of coordinate in arc-length parametrization. This is in particular an exponential barycenter: . In a Riemannian manifold, it would also be a minimizer of the square distance to the two points, necessarily a Fréchet mean if we assumed that all our points belong to a sufficiently small convex neighborhood.

### 1.3 A BCH-type formula on affine connection spaces

In order to analyze the pole ladder approximation of parallel transport in a local coordinate system, we need to compute the Taylor expansion of the parallel transport along a geodesic segment and of the composition of two Riemannian exponentials. Low order Taylor expansion of the Riemannian metric and of geodesic equations are traditional in Riemannian geometry since Gauss to establish the existence of normal coordinate system, Gauss lemma or the infinitesimal change of the volume of a geodesic ball due to the curvature. This is also what is used for the proof of [KMN00] for Schild’s ladder. However, establishing expansions above the order 3 is much more difficult, and it turns out that we need even more for the pole ladder. One solution was provided by [Bre97, Bre09] with the symbolic computer algebra system Cadabra. In this setting, the starting point is the (pseudo)-Riemannian metric, and it is not so easy to prove that the formulas continue to hold for affine connection spaces. More recently, a simpler algebraic formulation of higher order covariant derivatives in affine connection spaces [Gav06] led Gavrilov to obtain a combined Taylor expansion of the parallel transport and the composition of exponentials [Gav07]. After recalling this result in our notations, we adapt this formula to the analysis of pole ladder below.

We consider a small convex normal neighborhood of a point in an affine connection space . In a Riemannian manifold, we can chose a sufficiently small regular geodesic ball (see [Kar77, Ken90]). This means that any two points of the neighborhood can be can be joined by a unique geodesic segment and we can define the parallel transport of a tangent vector along this geodesic. The midpoint of each geodesic segment is also well defined in this neighborhood. The double exponential is defined as the mapping :

 expx(v,u)=expy(Πyxu)withy=expx(v).

For sufficiently small vectors , one can define the log of this expression (see Figure 4 for notations), which is a formal analog of the Baker-Campbell-Hausdorff (BCH) formula for Lie groups. Figure 4: The log of the composition of two exponentials (BCH-type formula) in a normal coordinate system at x.
###### Theorem 1 (Gavrilov’s BCH-type formula on manifolds [Gav06, Gav07]).

Let be an affine connection space with vanishing torsion and Riemannian curvature tensor . In a sufficiently small neighborhood of , the log of the double exponential:

 hx(v,u)=logx(exp(x)(v,u))=logx(expexpx(v)(Πexpx(v)xu)) (1)

has the following series expansion:

 hx(v,u)=v+u+16R(u,v)v+13R(u,v)u+112∇vR(u,v)v+124∇uR(u,v)v+524∇vR(u,v)u+112∇uR(u,v)u+O(∥u+w∥5). (2)

If the connection is not symmetric, then the torsion appears as a second order term and in all higher order terms.

In the above equations, denotes homogeneous terms of degree 5 and higher in and (the norm of can be taken arbitrarily), and the curvature tensor value is taken at . The full series expansion with torsion up to order 5 is provided in [Gav06].

### 1.4 Taylor expansion of pole ladder

We now have the tools to study the second version of pole ladder (algorithm 2). We consider that all points belong to a sufficiently small convex normal coordinate system centered at with the conventions of Figure 5. The geodesic starting at with tangent vector is a straight line in this chart going from point at time -1 to point at time 1. We denote by the parallel transport of the vector along the geodesic segment . This is equivalent to . We denote by the coordinates of in the normal coordinate system at : with . Using the previous BCH formula, we have:

 w=−hm(−v,u)=v−u−16R(u,v)v+13R(u,v)u+112∇vR(u,v)v−124∇uR(u,v)v−524∇vR(u,v)u+112∇uR(u,v)u+O(∥u+w∥5). (3)

Notice that the geodesic from to deviates from the straight line in the normal coordinate system at chart because of curvature.

We have an equivalent construction for the symmetric part of the structure: the symmetric of with respect to is . Taking the log at , we get the opposite of the pole ladder transport . In order to compare this vector to the exact parallel transport in a symmetric way, we now parallel transport it along the geodesic segment to get , that has to be compared with . Rephrasing this combination of transformations in the reverse way, we get that is the solution of with . Using Gavrilov’s formula, we can write the following polynomial expansion of in the variables and :

 w=hm(v,−u′)=v−u′−16R(u′,v)v+13R(u′,v)u′−112∇vR(u′,v)v+124∇′uR(u′,v)v+524∇vR(u′,v)u′−112∇u′R(u′,v)u′+O(∥u′+w∥5).

However, what we want to obtain is the polynomial expansion of with respect to the variables and . Such a polynomial can be written:

 u′=Av+Aw+Bvv+Bvw+Bww+Cvvv+Cvvw+Cvww+Cwww+Dvvvv+Dvvvw+Dvvww+Dvwww+Dwwww+O(∥u′+w∥5),

where each term is a multi-linear map in the indexed variables. Plugging this expression into Gavrilov’s formula above, we can identify each term iteratively for increasing orders: the first order is trivially and the second order terms have to vanish because there are none in the BCH-type formula. Introducing the third order into Gavrilov’s formula, we get:

 Cvvw+Cvww+Cwww=−16R(v−w,v)v+13R(v−w,v)(v−w)+O(∥v+w∥4).

Simplifying the expression thanks to the skew symmetries, we obtain:

 u′=v−w−16R(w,v)v+13R(w,v)w+O(∥v+w∥4).

The forth order is obtained by:

 Dvvvv+Dvvvw+Dvvww+Dvwww+Dwwww=−112∇vR(v−w,v)v+124∇(v−w)R(v−w,v)v+524∇vR(v−w,v)(v−w)−112∇(v−w)R(v−w,v)(v−w)+O(∥v+w∥5),

which gives after simplification:

 u′=v−w−16R(w,v)v+13R(w,v)w−112∇vR(w,v)v−124∇wR(w,v)v+18∇vR(w,v)w+112∇wR(w,v)w+O(∥v+w∥5).

Finally, we can substituting the value of from Equation (3) into this expression. This leads to:

 u′=u+512∇vR(u,v)u−16∇vR(u,v)v+112∇uR(u,v)v−16∇uR(u,v)u+O(∥v+u∥5). (4)
###### Theorem 2 (Pole ladder is a third order scheme for parallel transport).

In a sufficiently small convex normal neighborhood of a point in an affine connection space with symmetric connection, the error on one step of pole ladder to transport the vector along a geodesic segment of tangent vector (all quantities being parallel translated at the mid-point, see Figure 5) is:

 u′−u=112(∇vR(u,v)(5u−2v)+∇uR(u,v)(v−2u))+O(∥v+u∥5).

Instead of performing first a symmetry at and then at , an alternative version of pole ladder is to first perform a symmetry at and then at . One could think of averaging the two versions of the pole ladder to get a more accurate transport. The error on this alternative version is:

 u′′−u=−112(∇vR(u,v)(5u+v)+∇uR(u,v)(v+2u))+O(∥v+u∥5).

Unfortunately, we see that averaging does not kill all the third order terms.

## 2 Pole ladder in symmetric spaces

Theorem 2 shows that the main source of error in the pole ladder procedure is not the curvature itself but its covariant derivative. Since a vanishing covariant derivative is the characteristic of a locally symmetric space, this suggests that the error could vanish also in higher order terms. Thus, we investigate in this section local and global symmetric spaces with affine or Riemannian structure. The description of these spaces closely follows the definitions of [Pos01][Chapter 4]. The interested read should also refer to [Hel62]. We show that pole ladder is actually locally exact in one step, and even almost surely globally exact in Riemannian symmetric manifolds. The key feature is that the differential of the symmetry is the negative of the parallel transport, so that the parallel transport along a geodesic segment can be realized using the composition of two symmetries, which is called a transvection. Even if this result was not really known in the geometric computing community, the construction was already pictured for instance in [Tro94, Fig.5, p.168]).

### 2.1 Locally symmetric affine connection spaces

An affine connection space is said to be locally symmetric if the connection is torsion free and if the the curvature tensor is covariantly constant (). This second condition is equivalent to the preservation of the curvature tensor by the parallel transport along any path. When the manifold is endowed with a (pseudo-) Riemannian metric, one says that it is a locally symmetric (pseudo-) Riemannian space if the covariant derivative of its curvature tensor with respect to the Levi-Civita connection vanishes identically.

Locally symmetric affine connection spaces can be characterized by the geodesic symmetry which maps any point of a normal neighborhood of to the reverse point on the geodesic . The differential of this mapping is clearly a central symmetry on (a symmetric neighborhood of 0 in) the tangent space at : , and can thus be seen as a midpoint of the geodesic segment using the exponential barycenter definition since .

Using this geodesic symmetry, pole ladder can be rewritten as the composition of a symmetry at followed by a symmetry at (or alternatively as a symmetry at followed by a symmetry at ). Such a composition of two symmetries preserves the orientation and is called a transvection (or simply a translation). Moreover, its differential realizes the parallel transport:

###### Theorem 3 (Pole ladder is locally exact in locally symmetric spaces).

In a locally symmetric affine connection space with geodesic symmetry at point , the pole ladder transport along the geodesic segment of the tangent vector (identified to the geodesic segment ) is the tangent vector (identified to the geodesic segment ). Alternatively, one can start by the symmetry at to compute . For a sufficiently small time , both procedure realize exactly the parallel transport.

###### Proof.

We consider two points and that are sufficiently close in a symmetric normal neighborhood of their midpoint . The geodesic starting at with tangent vector remains in for small enough. The geodesic symmetry is an affine mapping in that maps geodesics to geodesics [Pos01][Chapter 3]. Thus, it maps to a geodesic starting at with tangent vector . Now, the key property is that the differential of the central symmetry at is is a linear maps from to which realizes exactly the (negative of) the parallel transport along the geodesic segment : . This property is described in [Hel62, Eq.(6), p.164] and in the proof of Proposition 4.3 of [Pos01, p.53-54]. The second symmetry is well defined for a sufficiently small time and reverses the sign. ∎

Because we only assumed so far a locally symmetric affine connection, the size of the normal neighborhood containing and and the maximal time cannot be specified without an auxiliary metric. More traditional symmetric (pseudo)-Riemannian spaces provide additional tools to extend this result.

### 2.2 Riemannian symmetric spaces

A connected affine connection space is a (globally) symmetric space if there exists at each point a smooth involution for which is an isolated point (i.e. a symmetry) and which is stable by composition: for any two points , . A symmetric space is geodesically complete and is obviously locally symmetric. The converse holds modulo discrete quotients: in essence, for each connected and geodesically complete locally symmetric space, there exists a symmetric covering space, which can be taken to be simply connected. Given a symmetric space, an even number of composition of symmetries generate a translation (or transvection) on the manifold. The Lie group of translations is acting transitively on , and its quotient space by the isotropy subgroup of one point of the manifold can be identified with the symmetric space itself.

Conversely, Cartan theorem shows that a homogeneous space with an involutive automorphism satisfying specific properties detailed below is a symmetric space. Because is an involution, its differential

has eigenvalues

which decompose the Lie algebra into the direct sum

: the +1-eigenspace

is the Lie algebra of , and the -eigenspace is an -invariant complement to in that can be identified to the tangent space at . These two conditions read in the Lie theoretic characterization and , which means that a symmetric space is a reductive homogeneous space. However, reductive homogeneous spaces are not all symmetric: the key feature is that (the composition of two symmetries on has to be a translation on ).

A natural structure on a symmetric space has to be invariant under its symmetry. This condition uniquely define a symmetric connection on each symmetric space. For instance, a connected Lie group is a symmetric space with the symmetry for two points of the group . The associated canonical connection is the mean (or symmetric) Cartan-Schouten connection which was used for defining bi-invariant means on Lie groups even in the absence of an invariant metric [PA12]. On a more general symmetric space , the canonical connection is induced by the symmetric Cartan-Schouten connection on the Lie group of translations , and geodesics of the symmetric space with this canonical connection are realized by the action of one-parameter subgroups of on a point of the manifold.

When is compact, one can moreover obtain an invariant metric on (simply choose a metric on the Lie algebra and average it over ). We obtain in this case a Riemannian symmetric manifold where the symmetry is now an isometry. The Levi-Civita connection of this invariant Riemannian metric coincides with the canonical connection. The invariant metric need not be uniquely define up to a scalar factor. For instance, that there is a one-parameter family of -invariant metrics on SPD matrices [Pen06].

With an invariant Riemannian metric, we can define the maximal injectivity domain of the exponential map: let be the set of all vectors such that is a minimizing geodesic for for some . This is an open star-shape domain containing limited by the tangential cut locus where geodesics stop to be minimizing. The image of the tangential cut locus by the exponential map is the cut locus , which has null measure with respect to the canonical Riemannian measure. The maximal injectivity domain is mapped diffeomorphically by to the normal neighborhood of in the manifold. The shortest distance from a point to its cut locus is the injection radius . Since a Riemannian symmetric space is homogeneous, this radius does not depend on the point and is equal to its infimum over all the points . These definitions allows us to characterize the exactness of pole ladder transport in symmetric spaces.

###### Theorem 4 (Pole ladder is almost surely exact in Riemannian symmetric manifolds).

In a Riemannian symmetric manifold, pole ladder is exact when and . This happens almost surely since the cut locus has null measure. More restrictive metric conditions of exactness are and .

###### Proof.

The first condition is needed to have a well defined mid-point. The second condition ensures that the geodesic does not meet the cut-locus of for , so that the initial tangent vector of the geodesic from to remains . Thanks to the global symmetry , similar conditions hold for and , so that is the exact parallel transport of along the geodesic segment . ∎

## References

• [A.A12] A.A. Schild’s ladder, August 2012.
• [ACPA06] Vincent Arsigny, Olivier Commowick, Xavier Pennec, and Nicholas Ayache. Statistics on Diffeomorphisms in a Log-Euclidean Framework. In Xavier Pennec and Sarang Joshi, editors, 1st MICCAI Workshop on Mathematical Foundations of Computational Anatomy: Geometrical, Statistical and Registration Methods for Modeling Biological Shape Variability, pages 14–15, Copenhagen, Denmark, October 2006.
• [Bre97] Leo Brewin. Riemann normal coordinates. Department of Mathematics Preprint, Monash University, Clayton, Vic. 3168, Australia, 1997.
• [Bre09] Leo Brewin. Riemann normal coordinate expansions using cadabra. Classical and Quantum Gravity, 26(17):175017, 2009.
• [Gav06] A V Gavrilov. Algebraic Properties of Covariant Derivative and Composition of Exponential Maps. Math. Tr., 9(1):3–20, 2006. English version in Siberian Advances in Mathematics, 2006, 16:3, 54–70.
• [Gav07] A. V. Gavrilov. The double exponential map and covariant derivation. Siberian Mathematical Journal, 48(1):56–61, January 2007.
• [Hel62] Sigurdur Helgason. Differential geometry and symmetric spaces, volume 12. Academic press, 1962.
• [HLP13] Søren Hauberg, François Lauze, and Kim Steenstrup Pedersen. Unscented Kalman Filtering on Riemannian Manifolds. Journal of Mathematical Imaging and Vision, 46(1):103–120, May 2013.
• [JDM18] Shuma Jia, Nicolas Duchateau, Pamela Moceri, Maxime Sermesant, and Xavier Pennec. SEPT: Surface-Enhanced Parallel Transport for the Statistical Analysis of Cardiac Sequences. Submitted, March 2018. Submitted.
• [Kar77] Hermann Karcher. Riemannian center of mass and mollifier smoothing. Communications in Pure and Applied Mathematics, 30:509–541, 1977.
• [Ken90] Wilfrid S. Kendall. Probability, convexity, and harmonic maps with small image I: uniqueness and fine existence. Proc. London Math. Soc., 61(2):371–406, 1990.
• [KMN00] Arkady Kheyfets, Warner A. Miller, and Gregory A. Newton. Schild’s Ladder Parallel Transport Procedure for an Arbitrary Connection. International Journal of Theoretical Physics, 39(12):2891–2898, December 2000.
• [LAP11] Marco Lorenzi, Nicholas Ayache, and Xavier Pennec. Schilds Ladder for the parallel transport of deformations in time series of images. In G. Szekely and H. Hahn, editors, IPMI - 22nd International Conference on Information Processing in Medical Images- 2011, volume 6801, pages 463–474, Kloster Irsee, Germany, July 2011. Springer. Honorable Mention (runner-up) for the Erbsmann Award.
• [LBCD17] Maxime Louis, Alexandre Bône, Benjamin Charlier, and Stanley Durrleman. Parallel Transport in Shape Analysis: A Scalable Numerical Scheme. In Frank Nielsen and Frédéric Barbaresco, editors, Geometric Science of Information, volume 10589, pages 29–37. Springer International Publishing, 2017.
• [LCJ17] Maxime Louis, Benjamin Charlier, Paul Jusselin, Susovan Pal, and Stanley Durrleman. A Fanning Scheme for the Parallel Transport Along Geodesics on Riemannian Manifolds. page 21, 2017. Preprint: https://hal.archives-ouvertes.fr/hal-01560787v2hal-01560787.
• [LP13a] Marco Lorenzi and Xavier Pennec. Efficient Parallel Transport of Deformations in Time Series of Images: from Schild’s to Pole Ladder. Journal of Mathematical Imaging and Vision, 50(1-2):5–17, October 2013.
• [LP13b] Marco Lorenzi and Xavier Pennec. Geodesics, Parallel Transport & One-parameter Subgroups for Diffeomorphic Image Registration.

International Journal of Computer Vision

, 105(2):111–127, November 2013.
• [MTW73] Charles W. Misner, Kip S. Thorne, and JohnArchibald Wheeler. Gravitation. W.H. Freeman and Compagny, 1973.
• [MY01] Michael I. Miller and Laurent Younes. Group Actions, Homeomorphisms, and Matching: A General Framework. International Journal of Computer Vision, 41(1-2):61–84, January 2001.
• [PA12] Xavier Pennec and Vincent Arsigny. Exponential Barycenters of the Canonical Cartan Connection and Invariant Means on Lie Groups. In Frederic Barbaresco, Amit Mishra, and Frank Nielsen, editors, Matrix Information Geometry, pages 123–168. Springer, May 2012.
• [Pen06] Xavier Pennec. Statistical Computing on Manifolds for Computational Anatomy. Habilitation à diriger des recherches, Université Nice Sophia-Antipolis, December 2006.
• [Pos01] M. M. Postnikov. Geometry VI: Riemannian Geometry. Encyclopedia of mathematical science. Springer, 2001.
• [Sch70] Alfred Schild. Tearing geometry to pieces: More on conformal geometry. unpublished lecture at Jan. 19 1970 Princeton Univesity relativity seminar, 1970.
• [Tro94] Valerii Vladimirovich Trofimov. Introduction to Geometry of Manifolds with Symmetry. Springer Netherlands, Dordrecht, 1994.
• [You07] Laurent Younes. Jacobi fileds in groups of diffeomorphisms and applciations. Quarterly of Applied Mathematics, 65(1):113–134, 2007.
• [YQWM08] Laurent Younes, Anqi Qiu, Raimond L. Winslow, and Michael I. Miller. Transport of Relational Structures in Groups of Diffeomorphisms. Journal of Mathematical Imaging and Vision, 32(1):41–56, September 2008.