Implementation of DIMSIMs for stiff differential systems

Authors:
Z. Jackiewicz
Affiliations:
Department of Mathematics, Arizona State University, Tempe, AZ
Venue:
Applied Numerical Mathematics
Year:
2002

Citing 9
Cited 3

Automatic selection of the initial step size for an ODE solver

Journal of Computational and Applied Mathematics
Solving ordinary differential equations I (2nd revised. ed.): nonstiff problems

Solving ordinary differential equations I (2nd revised. ed.): nonstiff problems
Diagonally-implicit multi-stage integration methods

Applied Numerical Mathematics
Construction of diagonally implicit general linear methods of type 1 and 2 for ordinary differential equations

Applied Numerical Mathematics
The MATLAB ODE Suite

SIAM Journal on Scientific Computing
A nonlinear optimization approach to the construction of general linear methods of high order

Journal of Computational and Applied Mathematics
Implementation of Diagonally Implicit Multistage Integration Methods for Ordinary Differential Equations

SIAM Journal on Numerical Analysis
Construction of high order diagonally implicit multistage integration methods for ordinary differential equations

Applied Numerical Mathematics
Exploiting structure in the construction of DIMSIMs

Journal of Computational and Applied Mathematics

Construction and Implementation of General Linear Methods for Ordinary Differential Equations: A Review

Journal of Scientific Computing
Error propagation of general linear methods for ordinary differential equations

Journal of Complexity
Perturbed MEBDF methods

Computers & Mathematics with Applications

Quantified Score

Hi-index	0.00

Visualization

Abstract

Some issues related to the implementation of diagonally implicit multistage integration methods for stiff differential systems are discussed. They include reliable estimation of the local discretization error, construction of continuous interpolants, solution of nonlinear systems of equations by simplified Newton iterations, choice of initial stepsize and order, and step and order changing strategy. Numerical results are presented which indicate that an experimental Matlab code based on type 2 methods of order one, two and three outperforms ode15s code from Matlab ODE suite on problems whose Jacobian has eigenvalues which are close to the imaginary axis.