Geometric integration algorithms on homogeneous manifolds

Abstract

Given an ordinary differential equation on a homogeneous manifold, one can construct a "geometric integrator'' by determining a compatible ordinary differential equation on the associated Lie group, using a Lie group integration scheme to construct a discrete time approximation of the solution curves in the group, and then mapping the discrete trajectories onto the homogeneous manifold using the group action. If the points of the manifold have continuous isotropy, a vector field on the manifold determines a continuous family of vector fields on the group, typically with distinct discretizations. If sufficient isotropy is present, an appropriate choice of vector field can yield improved capture of key features of the original system. In particular, if the algebra of the group is "full,'' then the order of accuracy of orbit capture (i.e. approximation of trajectories modulo time reparametrization) within a specified family of integration schemes can be increased by an appropriate choice of isotropy element. We illustrate the approach developed here with comparisons of several integration schemes for the reduced rigid body equations on the sphere.