Efficient numerical integration of Nth-order non-autonomous linear differential equations

• Philipp Bader ^[1] ; Sergio Blanes ^[2] ; Fernando Casas ^[3] ; Enrique Ponsoda ^[2]
  1. [1] La Trobe University

     La Trobe University

     Australia

     
  2. [2] Universidad Politécnica de Valencia

     Universidad Politécnica de Valencia

     Valencia, España

     
  3. [3] Universitat Jaume I

     Universitat Jaume I

     Castellón, España

     

  Mostrar afiliaciones +
• Localización: Journal of computational and applied mathematics, ISSN 0377-0427, Vol. 291, Nº 1 (1 January 2016), 2016, págs. 380-390
• Idioma: inglés
• DOI: 10.1016/j.cam.2015.02.052
• Enlaces
  • Texto completo
• Resumen

  • We consider the numerical integration of high-order linear non-homogeneous differential equations, written as first order homogeneous linear equations, and using exponential methods. Integrators like Magnus expansions or commutator-free methods belong to the class of exponential methods showing high accuracy on stiff or oscillatory problems, but the computation of the exponentials or their action on vectors can be computationally costly. The first order differential equations to be solved present a special algebraic structure (associated with the companion matrix) which allows to build new methods (hybrid methods between Magnus and commutator-free methods). The new methods are of similar accuracy as standard exponential methods with a reduced complexity. Additional parameters can be included into the scheme for optimization purposes. We illustrate how these methods can be obtained and present several sixth-order methods which are tested in several numerical experiments.