On the Linear Stability of Splitting Methods

  • Authors:

  • Sergio Blanes;Fernando Casas;Ander Murua

  • Affiliations:

  • Universitat Politécnica de Valencia, Instituto de Matemática Multidisciplinar, 46022, Valencia, Spain;Universitat Jaume I, Departament de Matemàtiques, 12071, Castellón, Spain;EHU/UPV, Konputazio Zientziak eta A.A. saila, Informatika Fakultatea, 12071, Donostia/San Sebastián, Spain

  • Venue:
  • Foundations of Computational Mathematics
  • Year:
  • 2008

Quantified Score

Hi-index 0.00

Visualization

Abstract

A comprehensive linear stability analysis of splitting methods is carried out by means of a 2×2 matrix K(x) with polynomial entries (the stability matrix) and the stability polynomial p(x) (the trace of K(x) divided by two). An algorithm is provided for determining the coefficients of all possible time-reversible splitting schemes for a prescribed stability polynomial. It is shown that p(x) carries essentially all the information needed to construct processed splitting methods for numerically approximating the evolution of linear systems. By conveniently selecting the stability polynomial, new integrators with processing for linear equations are built which are orders of magnitude more efficient than other algorithms previously available.