This paper is concerned with adaptive stiff solvers at low accuracy and complexity for systems of ordinary differential equations. The considered stiff solvers are: two second order Rosenbrock methods with low complexity, and the BDF method of the same order. For the adaptive algorithm we propose to use a monitor function defined by comparing a measure of the local variability of the solution times the used step size and the order of magnitude of the solution instead of the classical approach based on some local error estimation. This simple step-size selection procedure is implemented in order to control the behavior of the numerical solution. It is easily used to automatically adjust the step size, as the calculation progresses, until user-specified tolerance bounds for the introduced monitor function are fulfilled. This leads to important advantages in accuracy, efficiency and general ease-of-use. At the end of the paper we present two numerical tests which show the performance of the implementation of the stiff solvers, with the proposed adaptive procedure.

Adaptive stiff solvers at low accuracy and complexity

JANNELLI, Alessandra;FAZIO, Riccardo
2006-01-01

Abstract

This paper is concerned with adaptive stiff solvers at low accuracy and complexity for systems of ordinary differential equations. The considered stiff solvers are: two second order Rosenbrock methods with low complexity, and the BDF method of the same order. For the adaptive algorithm we propose to use a monitor function defined by comparing a measure of the local variability of the solution times the used step size and the order of magnitude of the solution instead of the classical approach based on some local error estimation. This simple step-size selection procedure is implemented in order to control the behavior of the numerical solution. It is easily used to automatically adjust the step size, as the calculation progresses, until user-specified tolerance bounds for the introduced monitor function are fulfilled. This leads to important advantages in accuracy, efficiency and general ease-of-use. At the end of the paper we present two numerical tests which show the performance of the implementation of the stiff solvers, with the proposed adaptive procedure.
File in questo prodotto:
Non ci sono file associati a questo prodotto.
Pubblicazioni consigliate

I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.

Utilizza questo identificativo per citare o creare un link a questo documento: https://hdl.handle.net/11570/1889763
 Attenzione

Attenzione! I dati visualizzati non sono stati sottoposti a validazione da parte dell'ateneo

Citazioni
  • ???jsp.display-item.citation.pmc??? ND
  • Scopus 20
  • ???jsp.display-item.citation.isi??? 25
social impact