Improving the Efficiency of Matrix Operations in the Numerical Solution of Stiff Ordinary Differential Equations

In the numerical solution of large stiff systems of ordinary differential equations, matrix operations associated with the solution of linear equations often dominate the solution time. A matrix factorization is suggested that will allow efficient updating after a change in stepsize or order. This updating technique is shown to be applicable to a wide variety of methods for stiff systems including multistep methods, Runge-Kutta methods, and methods using a rational function of a matrix The technique is particularly useful if the system is large and the Jacobian is dense Numerical results are included to illustrate the use of the technique.

[1]  W. H. Enright,et al.  Second Derivative Multistep Methods for Stiff Ordinary Differential Equations , 1974 .

[2]  J. H. Wilkinson,et al.  Similarity reduction of a general matrix to Hessenberg form , 1968 .

[3]  B. S. Garbow,et al.  Matrix Eigensystem Routines — EISPACK Guide , 1974, Lecture Notes in Computer Science.

[4]  J. Butcher.  Implicit Runge-Kutta processes , 1964 .

[5]  A. Hindmarsh,et al.  GEAR: ORDINARY DIFFERENTIAL EQUATION SYSTEM SOLVER. , 1971 .

[6]  W. H. Enright.  Optimal Second Derivative Methods for Stiff Systems , 1974 .

[7]  Edward J. Davison.  An algorithm for the computer simulation of very large dynamic systems , 1973 .

[8]  T. E. Hull,et al.  Comparing numerical methods for stiff systems of O.D.E:s , 1975 .

[9]  Alan C. Hindmarsh,et al.  A Polyalgorithm for the Numerical Solution of Ordinary Differential Equations , 1975, TOMS.

[10]  Jr. Roy Leonard Brown,et al.  Multi-Derivative Numerical Methods for the Solution of Stiff Ordinary Differential Equations , 1975 .

[11]  C. W. Gear,et al.  Algorithm 407: DIFSUB for solution of ordinary differential equations [D2] , 1971, Commun. ACM.

[12]  J. H. Wilkinson.  The algebraic eigenvalue problem , 1966 .

[13]  J. L. Blue,et al.  Rational approximations to matrix exponential for systems of stiff differential equations , 1970 .

[14]  A. C. Hindmarsh,et al.  Applications of EPISODE: an experimental package for the integration of systems of ordinary differential equations. [In FORTRAN for CDC 7600] , 1975 .