## Hamilton–Pontryagin integrators on Lie groups (2007)

Citations: | 21 - 6 self |

### BibTeX

@TECHREPORT{Bou-rabee07hamilton–pontryaginintegrators,

author = {Nawaf Bou-rabee and Jerrold E. Marsden},

title = {Hamilton–Pontryagin integrators on Lie groups},

institution = {},

year = {2007}

}

### OpenURL

### Abstract

Abstract In this paper, structure-preserving time-integrators for rigid body-type mechanical

### Citations

158 | Discrete mechanics and variational integrators
- Marsden, West
(Show Context)
Citation Context ...rategyFound Comput Math yields a Lagrangian analog of a well-known class of symplectic partitioned Runge– Kutta methods including the Lobatto IIIA-IIIB pair which generalize to higher-order accuracy =-=[8, 23, 32]-=-. In the Lie group context, one can generalize this strategy using either constrained or generalized coordinates. To use constrained coordinates, one treats the system as a holonomically constrained m... |

129 | The Euler-Poincarè equations and Semidirect Products with Applications to Continuum Theories - Holm, Marsden, et al. - 1998 |

115 |
Discrete versions of some classical integrable systems and factorization of matrix polynomials
- Moser
- 1991
(Show Context)
Citation Context ... coadjoint action on SO(3) and a numerical estimate of the gradient of the Hamiltonian. Variational integration techniques have been used to derive structure-preserving integrators on Lie groups; see =-=[2, 3, 21, 25, 33]-=-. Moser and Veselov derived a variational integrator for the free rigid body by embedding SO(3) in the linear space of 3 × 3 matrices, R 9 , and using Lagrange multipliers to constrain the matrices to... |

96 | Liegroup methods - Iserles, Munthe-Kaas, et al. - 1994 |

85 | Simulating Hamiltonian Dynamics - LEIMKUHLER, REICH - 2005 |

72 | Variational integrators and the Newmark algorithm for conservative and dissipative mechanical systems - Kane, Marsden, et al. |

71 | Runge–Kutta methods on Lie groups - Munthe-Kaas - 1998 |

61 | Mechanical integrators derived from a discrete variational principle, Physica D 106
- Wendlandt, Marsden
- 1997
(Show Context)
Citation Context ... coadjoint action on SO(3) and a numerical estimate of the gradient of the Hamiltonian. Variational integration techniques have been used to derive structure-preserving integrators on Lie groups; see =-=[2, 3, 21, 25, 33]-=-. Moser and Veselov derived a variational integrator for the free rigid body by embedding SO(3) in the linear space of 3 × 3 matrices, R 9 , and using Lagrange multipliers to constrain the matrices to... |

51 | Computations in a free Lie algebra - Munthe-Kaas, Owren - 1999 |

45 | Discrete Euler-Poincare and Lie-Poisson equations
- Marsden, Pekarsky, et al.
- 1999
(Show Context)
Citation Context ... coadjoint action on SO(3) and a numerical estimate of the gradient of the Hamiltonian. Variational integration techniques have been used to derive structure-preserving integrators on Lie groups; see =-=[2, 3, 21, 25, 33]-=-. Moser and Veselov derived a variational integrator for the free rigid body by embedding SO(3) in the linear space of 3 × 3 matrices, R 9 , and using Lagrange multipliers to constrain the matrices to... |

44 | E and Ratiu T S 1999 Introduction to Mechanics and Symmetry - Marsden |

43 |
Conserving algorithms for the dynamics of Hamiltonian systems on Lie groups
- Lewis, Simo
- 1994
(Show Context)
Citation Context ...ar space using holonomic constraints. Coadjoint-orbit and energy preserving methods of the Simo and Wong type that further preserve the symplectic structure were developed for SO(3) by Lewis and Simo =-=[16, 17]-=-. This was done by defining a one-parameter family of coadjoint and energy-preserving algorithms of the Simo and Wong type in which the free parameter is a functional. The function was specified so th... |

43 | Lie–Butcher theory for Runge–Kutta methods - Munthe-Kaas - 1995 |

43 | Integrable discrete-time systems and difference operators - Veselov - 1988 |

41 | Asynchronous variational integrators - Lew, Marsden, et al. |

33 | The Reduced Euler-Lagrange Equations
- Marsden, Scheurle
- 1993
(Show Context)
Citation Context ... an equivalence is established between critical points of s and S. IftheLagrangian is left-invariant, it is shown that this principle unifies the system’s Lie– Poisson and Euler–Poincaré descriptions =-=[6, 22]-=-. Since the reconstruction equation is a differential equation on a Lie group, one cannot directly discretize it by an RK method. However, one can discretize it using an s-stage Runge–Kutta–Munthe-Kaa... |

31 | Difference schemes for Hamiltonian formalism and symplectic geometry - Feng - 1986 |

30 | Variational time integrators
- LEW, MARSDEN, et al.
- 2003
(Show Context)
Citation Context ...nting that variational integrators correctly computeFound Comput Math other statistical quantities in long-time simulations. For example, in a simulation of a coupled spring-mass lattice, Lew et al. =-=[15]-=- found that variational integrators accurately compute the time-averaged instantaneous temperature (mean kinetic energy over all particles) over long-time intervals, whereas RK4 exhibits a artificial ... |

29 | Geometric numerical integration, volume 31 - Hairer, Lubich, et al. - 2006 |

28 | A Lie group variational integrator for the attitude dynamics of a rigid body with applications to the 3D pendulum
- Lee, Leok, et al.
(Show Context)
Citation Context ...p from the Lie algebra to the Lie group. These techniques were applied to bodies with attitude-dependent potentials, discrete optimal control of rigid bodies, and extended to higher-order accuracy in =-=[13, 14]-=-. Bobenko and Suris [2] considered a more general case where the symmetry group is a subgroup of the Lie group G in the context of semidirect Euler–Poincaré theory (see [9]). They did this by writing ... |

28 |
Unconditionally stable algorithms for rigid body dynamics that exactly preserve nergy and angular
- Simo, Wong
- 1991
(Show Context)
Citation Context ...tegies available to derive structure-preserving Lie group integrators; some of these are discussed here. One strategy involves generalizations of the classical Newmark algorithms to Lie groups due to =-=[30, 31]-=-. It was not apparent that the proposed Lie–Newmark methods had the necessary structure-preserving properties. In fact, Simo and Wong proposed another set of algorithms which preserve momentum by usin... |

27 | Dirac structures in Lagrangian mechanics part II: Variational structures
- Yoshimura, Marsden
- 2006
(Show Context)
Citation Context ...s paper is concerned with efficient, structure-preserving time integrators for mechanical systems whose configuration space is a Lie group, based on the Hamilton–Pontryagin (HP) variational principle =-=[11, 12, 20, 34, 35]-=-. This HP principle has many attractive theoretical properties; for instance, how it handles degenerate Lagrangian systems. The present paper shows that the HP viewpoint also provides a practical way ... |

23 | Lie group variational integrators for the full body problem
- Lee, Leok, et al.
- 2007
(Show Context)
Citation Context ...p from the Lie algebra to the Lie group. These techniques were applied to bodies with attitude-dependent potentials, discrete optimal control of rigid bodies, and extended to higher-order accuracy in =-=[13, 14]-=-. Bobenko and Suris [2] considered a more general case where the symmetry group is a subgroup of the Lie group G in the context of semidirect Euler–Poincaré theory (see [9]). They did this by writing ... |

23 | Equivariant constrained symplectic integration
- McLachlan, Scovel
- 1995
(Show Context)
Citation Context ...They also numerically demonstrated the method’s good performance crediting it to third-order accuracy in the discrete approximation to the Lie–Poisson structure. In related work, Mclachlan and Scovel =-=[24]-=- construct reduced, coadjoint-orbit preserving integrators by reducing G-equivariant integrators on T ∗ G obtained by embedding G in a linear space using holonomic constraints. Coadjoint-orbit and ene... |

22 | Geometric, variational integrators for computer animation
- KHAREVYCH, WEIWEI, et al.
- 2006
(Show Context)
Citation Context ...s paper is concerned with efficient, structure-preserving time integrators for mechanical systems whose configuration space is a Lie group, based on the Hamilton–Pontryagin (HP) variational principle =-=[11, 12, 20, 34, 35]-=-. This HP principle has many attractive theoretical properties; for instance, how it handles degenerate Lagrangian systems. The present paper shows that the HP viewpoint also provides a practical way ... |

22 |
On the dynamics in space of rods undergoing large motions - A geometrically exact approach
- Simo, Vu-Quot
- 1988
(Show Context)
Citation Context ...tegies available to derive structure-preserving Lie group integrators; some of these are discussed here. One strategy involves generalizations of the classical Newmark algorithms to Lie groups due to =-=[30, 31]-=-. It was not apparent that the proposed Lie–Newmark methods had the necessary structure-preserving properties. In fact, Simo and Wong proposed another set of algorithms which preserve momentum by usin... |

22 | Veselov [1991], Discrete versions of some classical integrable systems and factorization of matrix polynomials - Moser, P |

21 | 1998] The Euler–Poincaré equations and semidirect products with applications to continuum theories - Holm, Marsden, et al. |

21 | Symplectic-energy-momentum preserving variational integrators - Kane, Marsden, et al. - 1999 |

20 | Discrete time Lagrangian mechanics on Lie groups, with an application to the Lagrange top
- Bobenko, Suris
- 1999
(Show Context)
Citation Context |

19 | Discrete variational Hamiltonian mechanics
- Lall, West
(Show Context)
Citation Context ...s paper is concerned with efficient, structure-preserving time integrators for mechanical systems whose configuration space is a Lie group, based on the Hamilton–Pontryagin (HP) variational principle =-=[11, 12, 20, 34, 35]-=-. This HP principle has many attractive theoretical properties; for instance, how it handles degenerate Lagrangian systems. The present paper shows that the HP viewpoint also provides a practical way ... |

17 |
Hamiltonian methods of Runge-Kutta type and their variational interpretation
- Suris
- 1990
(Show Context)
Citation Context ...rategyFound Comput Math yields a Lagrangian analog of a well-known class of symplectic partitioned Runge– Kutta methods including the Lobatto IIIA-IIIB pair which generalize to higher-order accuracy =-=[8, 23, 32]-=-. In the Lie group context, one can generalize this strategy using either constrained or generalized coordinates. To use constrained coordinates, one treats the system as a holonomically constrained m... |

16 |
Almost Poisson integration of rigid body systems
- Austin, Krishnaprasad, et al.
- 1993
(Show Context)
Citation Context ... set of algorithms which preserve momentum by using the coadjoint action on SO(3) to advance the flow. Such integrators will be referred to as coadjoint-orbit preserving methods. Later, Austin et al. =-=[1]-=- showed that the midpoint rule member of the Lie–Newmark family with a Cayley reconstruction procedure was, in fact, a coadjoint-orbit preserving method for SO(3). They also numerically demonstrated t... |

16 | Error analysis for a class of methods for stiff non-linear initial value problems - Dahlquist - 1976 |

16 | Symplectic splitting methods for rigid body molecular dynamics - Dullweber, Leimkuhler, et al. - 1997 |

16 | On Cayley-transform methods for the discretization of Liegroup equations - Iserles - 2001 |

16 | Some aspects of the dynamics of Hamiltonian systems - MacKay - 1992 |

14 | Marsden [1997], Mechanical integrators derived from a discrete variational principle - Wendlandt, E |

12 | Discrete Lagrangian reduction, discrete Euler-Poincaré equations, and semidirect products
- Bobenko, Suris
- 1999
(Show Context)
Citation Context |

12 | Variational principles for Lie-Poisson and HamiltonPoincaré equations, Mosc
- Cendra, Marsden, et al.
(Show Context)
Citation Context ... an equivalence is established between critical points of s and S. IftheLagrangian is left-invariant, it is shown that this principle unifies the system’s Lie– Poisson and Euler–Poincaré descriptions =-=[6, 22]-=-. Since the reconstruction equation is a differential equation on a Lie group, one cannot directly discretize it by an RK method. However, one can discretize it using an s-stage Runge–Kutta–Munthe-Kaa... |

11 |
The Moser–Veselov equation
- Cardoso, Leite
- 2001
(Show Context)
Citation Context ...derivative: δ − δ∗ Dskew(g) · δ = 2 = (δg−1g) − (δg−1g) ∗ . 2 By definition of the right-trivialized tangent of τ −1 , it then follows that dskew(x)(y) = yτ(x) − (yτ(x))∗ . (4.27) 2 Cardoso and Leite =-=[5]-=- obtained the following proposition that explicitly determines τ(ξ). Moreover, they give necessary and sufficient conditions for its existence.Found Comput Math Proposition 4.10 Given ξ ∈ so(n), a sp... |

11 | Methods of Integration Which Preserve the Contact Transformation Property of the Hamilton Equations, N7-ONR-43906 - DeVogelaere |

9 |
Dirac structures and Lagrangian mechanics
- YOSHIMURA, MARSDEN
- 2006
(Show Context)
Citation Context |

9 | An overview of variational integrators - Lew, Marsden, et al. - 2004 |

7 |
Discrete mechanics and optimal control for constrained multibody dynamics
- Leyendecker, Ober-Blöbaum, et al.
- 2007
(Show Context)
Citation Context ...size h = 0.05, RK4 corrupts chaotic invariant sets while the lower-order accurate VE method preserves the structure of the benchmark their virtual work [23]), as well as discrete optimal control (see =-=[19]-=- and references therein). Altogether, this description of mechanics stands as a self-contained theory of mechanics akin to Hamiltonian, Lagrangian, or Newtonian mechanics. One of the distinguishing fe... |

7 | 2003], Variational principles for Lie-Poisson and Hamilton-Poincaré equations - Cendra, Marsden, et al. |

5 |
Conserving algorithms for the N-dimensional rigid body
- Lewis, Simo
- 1996
(Show Context)
Citation Context ...ar space using holonomic constraints. Coadjoint-orbit and energy preserving methods of the Simo and Wong type that further preserve the symplectic structure were developed for SO(3) by Lewis and Simo =-=[16, 17]-=-. This was done by defining a one-parameter family of coadjoint and energy-preserving algorithms of the Simo and Wong type in which the free parameter is a functional. The function was specified so th... |

4 | Numerical integration of Lie-Poisson systems while preserving coadjoint orbits and energy
- Engø, Faltinsen
- 2001
(Show Context)
Citation Context ...ng map defined a transformation which preserves the canonical symplectic form. Endowing coadjoint-orbit preserving methods with energy-preserving properties was also the subject of Engø and Faltinsen =-=[7]-=-. Specifically, they introduced integrators of the Runge–Kutta–Munthe-Kaas type that preserved coadjoint orbits and energy using the coadjoint action on SO(3) and a numerical estimate of the gradient ... |

4 |
Geometric Numerical Integration, 2nd edn
- HAIRER, LUBICH, et al.
- 2006
(Show Context)
Citation Context ...rategyFound Comput Math yields a Lagrangian analog of a well-known class of symplectic partitioned Runge– Kutta methods including the Lobatto IIIA-IIIB pair which generalize to higher-order accuracy =-=[8, 23, 32]-=-. In the Lie group context, one can generalize this strategy using either constrained or generalized coordinates. To use constrained coordinates, one treats the system as a holonomically constrained m... |

4 | Stochastic variational partitioned RungeKutta integrators for constraint systems. submitted - Bou-Rabee, Owahdi - 2007 |