Constant coefficient linear multistep methods with step density control

Research output: Contribution to journalArticle


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. (c) 2006 Elsevier B.V. All rights reserved.


Research areas and keywords

Subject classification (UKÄ) – MANDATORY

  • Mathematics


  • density control, non-uniform grid, step, step size selection, multistep method, variable step size method
Original languageEnglish
Pages (from-to)891-900
JournalJournal of Computational and Applied Mathematics
Issue number2
Publication statusPublished - 2007
Publication categoryResearch

Bibliographic note

The information about affiliations in this record was updated in December 2015. The record was previously connected to the following departments: Numerical Analysis (011015004)