The efficient solution of linear constant-coefficient systems of differential equations

Abstract

Most numerical techniques that have been proposed for the solution of constant-coefficient linear prob- Zems involve -either directly or indirectly-approxi mating a matrix exponential. For large systems these techniques can be very expensive to implement on digital computers. In this paper we describe how standard stiff ordinary differential equation (ODE) methods can be modified to take advantage of linear ity and thereby efficiently solve large linear prob lems. Methods based on backward differentiation formulas or second derivative formulas are particular ly suitable for this purpose. The implementation of these methods is discussed and numerical results are presented.

Keywords

GEAR linear differential equations,LINGEAR matrix exponential ordinary differential equation software SDLIN,stiff equations

Get full access to this article

View all access options for this article.

References

Addison, C. Implementation Considerations for a Method Based on Second Derivative Multistep Formulas MSc thesis Department of Computer Science University of

Toronto

1978

Brown, R.L. Numerical Integration of Linearized Stiff Ordinary Differential Equations In L. Lapidus and W. Schiesser , editors, Methods for Differential Systems (Academic Press, New York, 1976), pp. 39-45

Davison, E.J.

An Algorithm for the Computer Simulation of Very Large

Dynamic Systems Automatica vol. 9 1973 pp. 665-675

Enright, W.H.

Second Derivative Multistep Methods for Stiff Ordinary Differential Equations SIAM

Journal of Numerical Analysis vol. 11 1974 pp. 321-332

Enright, W.H. Optimal Second Derivative Methods for Stiff Equations In R. Willoughby , editor, Stiff Differential Systems (Plenum Press, New York, 1975), pp. 95-109

Enright, W.H. Hull, T.E. Lindberg, B.

Comparing Numerical Methods for Stiff Systems of Ordinary

Differential Equations BIT vol. 15 1975 pp. 10-48

Hindmarsh, A.C. GEAR: Ordinary Differential Equation System Solver Lawrence Livermore Laboratory Report UCID-30001 revision 3 Livermore, California 1974

Papathomas, T. Wing, O.

Sparse Hessenberg Reduction and the Eigenvalue Problem for Large Sparse Matrices

IEEE Transactions on Circuits and Systems vol. CAS-23 1976 pp. 739-744

10.

Tendler, J.M. Bickart, T.A. Picel, Z. A Stiffly Stable Integration Process Using Cyclic Composite Methods Submitted to Transactions on Mathematical Software 1977

11.

Tewarson, R.P.

on the Reduction of a Sparse Matrix to Hessenberg Form

International Journal of Computer Mathematics vol. 2 1970 pp. 283-295