A nonlinear optimization approach to the construction of general linear methods of high order

We describe the construction of diagonally implicit multistage integration methods of order and stage order p = q = 7 and p = q = 8 for ordinary differential equations. These methods were obtained using state-of-the-art optimization methods, particularly variable-model trust-region least-squares algorithms.

[1]  John C. Butcher,et al.  Diagonally-implicit multi-stage integration methods , 1993 .

[2]  J. Butcher,et al.  Generalized Padé approximations to the exponential function , 1992 .

[3]  John E. Dennis,et al.  An Adaptive Nonlinear Least-Squares Algorithm , 1977, TOMS.

[4]  Roy E. Welsch,et al.  Algorithm 717: Subroutines for maximum likelihood and quasi-likelihood estimation of parameters in nonlinear regression models , 1993, TOMS.

[5]  Ernst Hairer,et al.  Solving Ordinary Differential Equations I: Nonstiff Problems , 2009 .

[6]  Zdzislaw Jackiewicz,et al.  Implementation of Diagonally Implicit Multistage Integration Methods for Ordinary Differential Equations , 1997 .

[7]  John E. Dennis,et al.  Numerical methods for unconstrained optimization and nonlinear equations , 1983, Prentice Hall series in computational mathematics.

[8]  S. Oren Self-scaling variable metric algorithms without line search for unconstrained minimization , 1973 .

[9]  Zdzislaw Jackiewicz,et al.  Construction of diagonally implicit general linear methods of type 1 and 2 for ordinary differential equations , 1996 .

[10]  J. Lambert Computational Methods in Ordinary Differential Equations , 1973 .

[11]  Zdzislaw Jackiewicz,et al.  Diagonally implicit general linear methods for ordinary differential equations , 1993 .

[12]  Jorge J. Moré,et al.  The Levenberg-Marquardt algo-rithm: Implementation and theory , 1977 .

[13]  A. D. Hall,et al.  The PORT Mathematical Subroutine Library , 1978, TOMS.

[14]  A. Iserles On the Generalized Pad Approximations to the Exponential Function , 1979 .

[15]  John E. Dennis,et al.  Algorithm 573: NL2SOL—An Adaptive Nonlinear Least-Squares Algorithm [E4] , 1981, TOMS.