# zbMATH — the first resource for mathematics

Exploiting the separability in the solution of systems of linear ordinary differential equations. (English) Zbl 0683.65052
The first order system: $$y'=A(t)y+b(t),\quad y(a)=y_ 0,$$ is called separable if the number of stiff eigenvalues of A(t) is k with $$k<<s$$. In this paper the method of the first author for stiff systems [Lect. Notes Math. 1066, 30-43 (1984; Zbl 0542.65037) and BIT 23, 329-345 (1983; Zbl 0523.65053)] is implemented supposing k known. The efficiency of the code as function of k/s is systematically analyzed and a family of test examples is presented.
Reviewer: A.de Castro

##### MSC:
 65L05 Numerical methods for initial value problems 34A30 Linear ordinary differential equations and systems, general
LINPACK; ODEPACK
Full Text:
##### References:
 [1] Lambert, J.D., Computational methods in ordinary differential equations, (1973), Wiley London · Zbl 0258.65069 [2] Söderlind, G., On the efficient solution of nonlinear equations in numerical methods for stiff differential systems, () · Zbl 0771.65036 [3] Enright, W.H.; Hull, T.E.; Lindberg, B., Comparing methods for stiff systems of odes, Bit, 15, 10-48, (1975) · Zbl 0301.65040 [4] Alexander, R., Diagonally implicit Runge-Kutta methods for stiff odes, SIAM J. numer. analysis, 14, 1006-1021, (1977) · Zbl 0374.65038 [5] Dahlquist, G.; Söderlind, G., Some problems related to stiff nonlinear differential systems, (), 57-74 [6] Robertson, H.H., Some factors affecting the efficiency of stiff integration routines, (), 279-301 [7] Enright, W.H.; Kamel, M.S., Automatic partitioning of stiff systems and exploiting the resulting structures, ACM trans. math. software, 5, 374-385, (1979) · Zbl 0429.65061 [8] Watkins, D.S.; HansonSmith, R.W., The numerical solution of separably stiff systems by precise partitioning, ACM trans. math. software, 9, 293-301, (1983) · Zbl 0517.65052 [9] Björck, Å., A block QR algorithm for partitioning stiff differential equations, Bit, 23, 329-345, (1983) [10] Björck, Å., Some methods for separating stiff components in initial value problems, (), 30-43 · Zbl 0542.65037 [11] Brown, P.N.; Hindmarsh, A.C., Matrix-free methods for stiff systems of odes, SIAM J. numer. analysis, 23, 610-638, (1986) · Zbl 0615.65078 [12] Chan, T.F.; Jackson, K.R., The use of iterative linear-equation solvers in codes for large systems of stiff IVPs for odes, SIAM J. scient. statist. comput., 7, 378-417, (1986) · Zbl 0611.65050 [13] Gear, C.W.; Saad, Y., Iterative solution of linear equations in ODE codes, SIAM J. scient. statist. comput., 4, 583-601, (1983) · Zbl 0541.65051 [14] Miranker, W.L.; Chern, I-L., Dichotomy and conjugate gradients in the stiff initial value problem, Linear algebra applic., 36, 57-77, (1981) · Zbl 0488.65026 [15] Zlatev, Z.; Wasniewski, J.; Schaumburg, K., Subroutine DENS1 for solving stiff systems of linear ordinary differential equations (basic algorithms, documentation, demonstration programs), () [16] Zlatev, Z., Modified diagonally implicit Runge-Kutta methods, SIAM J. scient. statist. comput., 2, 321-334, (1981) · Zbl 0475.65040 [17] Zlatev, Z.; Wasniewski, J.; Schaumburg, K., Numerical treatment of models arising in nuclear magnetic resonance spectroscopy, Adv. eng software, 8, 223-233, (1986) [18] Stewart, G.W., Introduction to matrix computations, (1973), Academic Press London · Zbl 0302.65021 [19] Wilkinson, J.H., The algebraic eigenvalue problem, (1965), Oxford University Press London · Zbl 0258.65037 [20] Lancaster, P., Theory of matrices, (1969), Academic Press London · Zbl 0186.05301 [21] Dahlquist, G., On transformation of graded matrices, with applications to stiff odes, Numer. math., 47, 363-385, (1985) · Zbl 0599.65020 [22] Zlatev, Z.; Wasniewski, J.; Schaumburg, K., Exploiting the sparsity in the solution of linear ordinary differential equations, Comput. math. applic., 11, 1069-1087, (1985) · Zbl 0576.65062 [23] Dongarra, J.J.; Eisenstat, S.C., Squeezing the most of an algorithm in cray FORTRAN, ACM trans. math. software, 10, 221-230, (1984) [24] Dongarra, J.J.; Bunch, J.R.; Moler, C.B.; Stewart, G.W., LINPACK users’ guide, (1979), SIAM Philadelphia · Zbl 0476.68025 [25] Lawson, C.; Hanson, R.; Kincaid, D.; Krogh, F., Basic linear algebra subprograms for FORTRAN usage, ACM trans. math. software, 5, 308-323, (1979) · Zbl 0412.65022 [26] Hindmarsh, A.C., ODEPACK, a systemized collection of ODE solvers, () [27] Østerby, O.; Zlatev, Z., Direct methods for sparse matrices, (1983), Springer Berlin · Zbl 0516.65011 [28] Curtis, A.R., Jacobian matrix properties and their impact on choice of software for stiff ODE systems, IMA J. numer. analysis, 3, 397-415, (1983) · Zbl 0534.65039 [29] McRae, G.J.; Goodin, W.R.; Seinfeld, J.H., Numerical solution of the atmospheric diffusion equation for chemically reacting flows, J. comput. phys., 45, 1-42, (1982) · Zbl 0502.76098
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. It attempts to reflect the references listed in the original paper as accurately as possible without claiming the completeness or perfect precision of the matching.