High Order Optimized Geometric Integrators for Linear Differential Equations

In this paper new integration algorithms based on the Magnus expansion for linear differential equations up to eighth order are obtained. These methods are optimal with respect to the number of commutators required. Starting from Magnus series, integration schemes based on the Cayley transform an the Fer factorization are also built in terms of univariate integrals. The structure of the exact solution is retained while the computational cost is reduced compared to similar methods. Their relative performance is tested on some illustrative examples.

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

[2]  Luciano Lopez,et al.  The Cayley transform in the numerical solution of unitary differential systems , 1998, Adv. Comput. Math..

[3]  Elena Celledoni,et al.  Complexity Theory for Lie-Group Solvers , 2002, J. Complex..

[4]  Arieh Iserles,et al.  Geometric integration: numerical solution of differential equations on manifolds , 1999, Philosophical Transactions of the Royal Society of London. Series A: Mathematical, Physical and Engineering Sciences.

[5]  Arne Marthinsen,et al.  Quadrature methods based on the Caylay transform , 2001 .

[6]  H. Yoshida.  Construction of higher order symplectic integrators , 1990 .

[7]  A. Iserles,et al.  Lie-group methods , 2000, Acta Numerica.

[8]  Arieh Iserles,et al.  On Cayley-Transform Methods for the Discretization of Lie-Group Equations , 2001, Found. Comput. Math..

[9]  G. Bodenhausen,et al.  Principles of nuclear magnetic resonance in one and two dimensions , 1987 .

[10]  Fernando Casas,et al.  Improved High Order Integrators Based on the Magnus Expansion , 2000 .

[11]  Antonella Zanna,et al.  Collocation and Relaxed Collocation for the Fer and the Magnus Expansions , 1999 .

[12]  H. Munthe-Kaas,et al.  Computations in a free Lie algebra , 1999, Philosophical Transactions of the Royal Society of London. Series A: Mathematical, Physical and Engineering Sciences.

[13]  A. Iserles,et al.  On the solution of linear differential equations in Lie groups , 1999, Philosophical Transactions of the Royal Society of London. Series A: Mathematical, Physical and Engineering Sciences.

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

[15]  Antonella Zanna,et al.  On the Dimension of Certain Graded Lie Algebras Arising in Geometric Integration of Differential Equations , 2000 .

[16]  J Rosyx.  Magnus and Fer expansions for matrix differential equations: the convergence problem , 1998 .

[17]  N. Bourbaki.  Lie groups and Lie algebras , 1998 .

[18]  A. Zanna.  The Fer expansion and time-symmetry: a Strang-type approach , 2001 .