Constant coefficient linear multistep methods with step density control

  • Authors:
  • Carmen Arévalo;Gustaf Söderlind;José Diaz López

  • Affiliations:
  • Centre for Mathematical Sciences, Lund University, Box 118, SE-221 00 Lund, Sweden;Centre for Mathematical Sciences, Lund University, Box 118, SE-221 00 Lund, Sweden;Dynasim AB, Research Park Ideon, SE-223 70 Lund, Sweden

  • Venue:
  • Journal of Computational and Applied Mathematics
  • Year:
  • 2007

Quantified Score

Hi-index 7.29

Visualization

Abstract

In linear multistep methods with variable step size, the method's coefficients are functions of the step size ratios. The coefficients therefore need to be recomputed on every step to retain the method's proper order of convergence. An alternative approach is to use step density control to make the method adaptive. If the step size sequence is smooth, the method can use constant coefficients without losing its order of convergence. The paper introduces this new adaptive technique and demonstrates its feasibility with a few test problems. The technique works in perfect agreement with theory for a given step density function. For practical use, however, the density must be generated with data computed from the numerical solution. We introduce a local error tracking controller, which automatically adapts the density to computed data, and demonstrate in computational experiments that the technique works well at least up to fourth-order methods.