Collocation and Relaxed Collocation for the Fer and the Magnus Expansions

We consider the Fer and the Magnus expansions for the numerical solution of the nonlinear matrix Lie-group ODE $y' = \gamma(t,y) y$, $y(0)= y_0,$ where y evolves in a matrix Lie group G and $\gamma(t,y)$ is in the Lie algebra %$\hbox{\Fr{g}}$. $\ghe$. Departing from a geometrical approach that distinguishes between those operations performed in the group and those performed in the tangent space, we construct Lie-group invariant methods based on collocation. We prove that, as long as the two expansions are correctly truncated, the collocation nodes $c_1, c_2, \ldots, c_{\nu}$ yield numerical methods whose order is the same as in the classical setting. We also relax the collocation conditions, thereby devising explicit methods of order 3. To conclude, we discuss the proposed methods in a numerical experiment that arises in Hamiltonian mechanics, comparing the results with projection techniques.

[1]  Arieh Iserles.  Multistep methods on manifolds , 2001 .

[2]  Simulation of Ordinary Di erentialEquations on Manifolds | Some Numerical Experiments and Veri cations , 1996 .

[3]  Arne Marthinsen,et al.  Simulation of ordinary differential equations on manifolds : some numerical experiments and verifications , 1997 .

[4]  D. Higham.  Time-stepping and preserving orthonormality , 1997 .

[5]  Antonella Zanna.  The method of iterated commutators for ordinary differential equations on Lie groups , 1996 .

[6]  V. Varadarajan.  Lie groups, Lie algebras, and their representations , 1974 .

[7]  C. T. Mullis,et al.  Inverse eigenvalue problem for real symmetric Toeplitz matrices , 1988, ICASSP-88., International Conference on Acoustics, Speech, and Signal Processing.

[8]  M. Chu.  A List of Matrix Flows with Applications , 1994 .

[9]  Antonella Zanna,et al.  Iterated commutators, Lie's reduction method and ordinary differential equations on matrix Lie groups , 1997 .

[10]  P. Olver.  Applications of Lie Groups to Differential Equations , 1986 .

[11]  Antonella Zanna,et al.  Numerical integration of differential equations on homogeneous manifolds , 1997 .

[12]  W. Magnus.  On the exponential solution of differential equations for a linear operator , 1954 .

[13]  P. Lax.  INTEGRALS OF NONLINEAR EQUATIONS OF EVOLUTION AND SOLITARY WAVES. , 1968 .

[14]  G. Wanner,et al.  Perturbed collocation and Runge-Kutta methods , 1981 .

[15]  Integration Methods Based on Rigid Frames , 1999 .

[16]  A. Perelomov.  Integrable systems of classical mechanics and Lie algebras , 1989 .

[17]  P. Crouch,et al.  Numerical integration of ordinary differential equations on manifolds , 1993 .

[18]  R. Russell,et al.  Unitary integrators and applications to continuous orthonormalization techniques , 1994 .

[19]  J. C. Simo,et al.  Unconditionally stable algorithms for rigid body dynamics that exactly preserve energy and momentum , 1991 .

[20]  J. M. Sanz-Serna,et al.  Numerical Hamiltonian Problems , 1994 .

[21]  F. Casas.  Solution of linear partial differential equations by Lie algebraic methods , 1996 .

[22]  B. Owren,et al.  The Newton Iteration on Lie Groups , 2000 .

[23]  戸田 盛和.  Theory of nonlinear lattices , 1981 .

[24]  A. Iserles,et al.  University of Cambridge Qualitative Numerical Analysis of Ordinary Diierential Equations , 1995 .

[25]  Arieh Iserles,et al.  Solving linear ordinary differential equations by exponentials of iterated commutators , 1984 .

[26]  Antonella Zanna,et al.  Numerical Integration of Di erential Equations on Homoge-neous Manifolds , 1996 .

[27]  H. Munthe-Kaas.  High order Runge-Kutta methods on manifolds , 1999 .

[28]  Antonella Zanna,et al.  Numerical solution of isospectral flows , 1997, Math. Comput..