×

Solving initial value problems for ordinary differential equations by two approaches: BDF and piecewise-linearized methods. (English) Zbl 1198.65142

Summary: Many scientific and engineering problems are described using Ordinary Differential Equations (ODEs), where the analytic solution is unknown. Much research has been done by the scientific community on developing numerical methods which can provide an approximate solution of the original ODE. In this work, two approaches have been considered based on BDF and Piecewise-linearized Methods. The approach based on BDF methods uses a Chord-Shamanskii iteration for computing the nonlinear system which is obtained when the BDF schema is used. Two approaches based on piecewise-linearized methods have also been considered. These approaches are based on a theorem proved in this paper which allows to compute the approximate solution at each time step by means of a block-oriented method based on diagonal Padé approximations. The difference between these implementations is in using or not using the scale and squaring technique.
Five algorithms based on these approaches are developed. MATLAB and Fortran versions of the above algorithms have been developed, comparing both precision and computational costs. BLAS and LAPACK libraries have been used in Fortran implementations. In order to compare in equality of conditions all implementations, algorithms with fixed step have been considered. Four of the five case studies analyzed come from biology and chemical kinetics stiff problems. Experimental results show the advantages of the proposed algorithms, especially when they are integrating stiff problems.

MSC:

65L99 Numerical methods for ordinary differential equations
34-04 Software, source code, etc. for problems pertaining to ordinary differential equations
PDFBibTeX XMLCite
Full Text: DOI Link

References:

[1] Hairer, E.; Wanner, G., Solving ordinary differential equations II. Stiff and differential-algebraic problems, (Springer Series in Computational Mathematics, vol. 14 (1996), Springer-Verlag) · Zbl 0729.65051
[2] C.F. Curtiss, J.O. Hirschfelder, Integration of stiff equations, in: Proc. Nat. Acad. Sci., vol. 38, 1952, pp. 235-243; C.F. Curtiss, J.O. Hirschfelder, Integration of stiff equations, in: Proc. Nat. Acad. Sci., vol. 38, 1952, pp. 235-243 · Zbl 0046.13602
[3] Hindmarsh, A. C., Lsode and lsodi, two new initial value ordinary differential equation solvers, ACM-Signum Newslett., 15, 1011 (1980)
[4] Gear, C. W., Algorithm 407-difsub for solution of ordinary differential equations, Comm. ACM, 14, 185-190 (1971)
[5] Brown, P. N.; Byrne, G. D.; Hindmarsh, A. C., Vode: A variable-coefficient ode solver, SIAM J. Sci. Statist. Comput., 10, 1038-1051 (1989) · Zbl 0677.65075
[6] Ramos, J. I.; García-López, C. M., Piecewise-linearized methods for initial-value problems, Applied Mathematics and Computation, 82, 273-302 (1997) · Zbl 0870.65066
[7] C.M. García-López, Métodos de linealización para la resolución numérica de ecuaciones diferenciales, Ph.D. thesis, Departamento de Lenguajes y Ciencias de la Computación, Universidad de Málaga, 1998; C.M. García-López, Métodos de linealización para la resolución numérica de ecuaciones diferenciales, Ph.D. thesis, Departamento de Lenguajes y Ciencias de la Computación, Universidad de Málaga, 1998
[8] Ascher, U. M.; Petzold, L. R., Computer Methods for Ordinary Differential Equations and Differential-Algebraic Equations (1998), SIAM · Zbl 0908.65055
[9] C.T. Kelley, Iterative Methods for Linear and Nonlinear Equations, Philadelphia, 1995; C.T. Kelley, Iterative Methods for Linear and Nonlinear Equations, Philadelphia, 1995 · Zbl 0832.65046
[10] García-López, C. M., Piecewise-linearized and linearized \(θ\)-methods for ordinary and partial differential equation problems, Computer & Mathematics with Applications, 45, 351-381 (2003) · Zbl 1035.65103
[11] Golub, G. H.; Loan, C. V., Matrix Computations, Johns Hopkins Studies in Mathematical Sciences (1996), The Johns Hopkins University Press
[12] Moler, C. B.; Loan, C. V., Nineteen dubious ways to compute the exponential of a matrix, twenty-five years later*, SIAM Review, 45, 3-49 (2003) · Zbl 1030.65029
[13] N.J. Higham, The scaling and squaring method for the matrix exponential revisited, Tech. Rep. 452, Manchester Centre for Computational Mathematics, 2004; N.J. Higham, The scaling and squaring method for the matrix exponential revisited, Tech. Rep. 452, Manchester Centre for Computational Mathematics, 2004 · Zbl 1178.65040
[14] Shampine, L. S.; Gladwell, I.; Thomson, S., Solving ODEs with MAtlab (2003), Cambridge University Press · Zbl 1079.65144
[15] S. Graphics, SGI Altix Applications Development an Optimization, Release August 1st Edition, 2003; S. Graphics, SGI Altix Applications Development an Optimization, Release August 1st Edition, 2003
[16] Dongarra, J.; Croz, J. D.; Hammarling, S.; Hanson, R. J., An extended set of FORTRAN Basic Linear Algebra Subroutines, ACM Trans. Math. Software, 14, 1-17 (1988) · Zbl 0639.65016
[17] Anderson, E.; Bai, Z.; Bischof, C.; Demmel, J.; Dongarra, J.; Croz, J. D.; Greenbaum, A.; Hammarling, S.; McKenney, A.; Ostrouchov, S.; Sorensen, D., LAPACK Users’ Guide (1992), SIAM · Zbl 0843.65018
[19] W.M. Lioen, J.J.B. de Swart, Test set for initial value problems solvers, release 2.0, December 1998; W.M. Lioen, J.J.B. de Swart, Test set for initial value problems solvers, release 2.0, December 1998 · Zbl 0914.65074
[20] Gropp, W.; Lusk, E.; Skjellum, A., Using MPI: Portable Parallel Programming with the Message-Passing Interface (1994), MIT Press
[21] J.J. Dongarra, R.A.V.D. Geijn, Two dimensional basic linear algebra communications subprograms, Tech. rep., Department of Computer Science, University of Tennessee, 1991; J.J. Dongarra, R.A.V.D. Geijn, Two dimensional basic linear algebra communications subprograms, Tech. rep., Department of Computer Science, University of Tennessee, 1991
[22] J. Choi, J. Dongarra, S. Ostrouchov, A. Petitet, D. Walker, A proposal for a set of parallel basic linear algebra subprogram, Tech. Rep. UT-CS-95-292, Department of Computer Science, University of Tennessee, 1995; J. Choi, J. Dongarra, S. Ostrouchov, A. Petitet, D. Walker, A proposal for a set of parallel basic linear algebra subprogram, Tech. Rep. UT-CS-95-292, Department of Computer Science, University of Tennessee, 1995
[23] Blackford, L. S.; Choi, J.; Cleary, A.; D’Azevedo, E.; Demmel, J.; Dhillon, I., ScaLAPACK Users’ Guide (1997), SIAM · Zbl 0886.65022
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. In some cases that data have been complemented/enhanced by data from zbMATH Open. This attempts to reflect the references listed in the original paper as accurately as possible without claiming completeness or a perfect matching.