zbMATH — the first resource for mathematics

Examples
 Geometry Search for the term Geometry in any field. Queries are case-independent. Funct* Wildcard queries are specified by * (e.g. functions, functorial, etc.). Otherwise the search is exact. "Topological group" Phrases (multi-words) should be set in "straight quotation marks". au: Bourbaki & ti: Algebra Search for author and title. The and-operator & is default and can be omitted. Chebyshev | Tschebyscheff The or-operator | allows to search for Chebyshev or Tschebyscheff. "Quasi* map*" py: 1989 The resulting documents have publication year 1989. so: Eur* J* Mat* Soc* cc: 14 Search for publications in a particular source with a Mathematics Subject Classification code (cc) in 14. "Partial diff* eq*" ! elliptic The not-operator ! eliminates all results containing the word elliptic. dt: b & au: Hilbert The document type is set to books; alternatively: j for journal articles, a for book articles. py: 2000-2015 cc: (94A | 11T) Number ranges are accepted. Terms can be grouped within (parentheses). la: chinese Find documents in a given language. ISO 639-1 language codes can also be used.

Operators
 a & b logic and a | b logic or !ab logic not abc* right wildcard "ab c" phrase (ab c) parentheses
Fields
 any anywhere an internal document identifier au author, editor ai internal author identifier ti title la language so source ab review, abstract py publication year rv reviewer cc MSC code ut uncontrolled term dt document type (j: journal article; b: book; a: book article)
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 ODE 34-04 Machine computation, programs (ordinary differential equations)
Software:
Matlab; BLAS; DIFSUB; VODE; ScaLAPACK; KELLEY
References:
 [1] Hairer, E.; Wanner, G.: Solving ordinary differential equations II. Stiff and differential-algebraic problems, Springer series in computational mathematics 14 (1996) [2] C.F. Curtiss, J.O. Hirschfelder, Integration of stiff equations, in: Proc. Nat. Acad. Sci., vol. 38, 1952, pp. 235 – 243 · Zbl 0046.13602 · doi:doi:10.1073/pnas.38.3.235 [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 · doi:doi:10.1137/0910062 [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 · doi:doi:10.1016/S0096-3003(96)00040-9 [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 [8] Ascher, U. M.; Petzold, L. R.: Computer methods for ordinary differential equations and differential-algebraic equations, (1998) · Zbl 0908.65055 [9] C.T. Kelley, Iterative Methods for Linear and Nonlinear Equations, Philadelphia, 1995 [10] García-López, C. M.: Piecewise-linearized and linearized $\theta$-methods for ordinary and partial differential equation problems, Computer & mathematics with applications 45, 351-381 (2003) · Zbl 1035.65103 · doi:doi:10.1016/S0898-1221(03)80023-8 [11] Golub, G. H.; Loan, C. V.: Matrix computations, Johns hopkins studies in mathematical sciences (1996) [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 · doi:doi:10.1137/S00361445024180 [13] N.J. Higham, The scaling and squaring method for the matrix exponential revisited, Tech. Rep. 452, Manchester Centre for Computational Mathematics, 2004 [14] Shampine, L. S.; Gladwell, I.; Thomson, S.: Solving odes with Matlab, (2003) [15] 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 · doi:doi:10.1145/42288.42291 · doi:http://www.acm.org/pubs/contents/journals/toms/1988-14/ [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) [18] [19] W.M. Lioen, J.J.B. de Swart, Test set for initial value problems solvers, release 2.0, December 1998 [20] Gropp, W.; Lusk, E.; Skjellum, A.: Using MPI: portable parallel programming with the message-passing interface, (1994) [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 [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 [23] Blackford, L. S.; Choi, J.; Cleary, A.; D’azevedo, E.; Demmel, J.; Dhillon, I.: Scalapack users’ guide, (1997)