An MEBDF package for the numerical solution of large sparse systems of stiff initial value problems. (English) Zbl 0983.65081

Summary: An efficient algorithm for the numerical integration of large sparse systems of stiff initial value ordinary differential equations and differential-algebraic equations is described. The algorithm is constructed by embedding a standard sparse linear algebraic equation solver into a suitably modified MEBDF code. An important practical application of this algorithm is in the numerical solution of time dependent partial differential equations, particularly in two or more space dimensions, using the method of lines. A code based on this algorithm is illustrated by application to several problems of practical interest and its performance is compared to that of the standard code LSODES.


65L05 Numerical methods for initial value problems involving ordinary differential equations
34A34 Nonlinear ordinary differential equations and systems
34A09 Implicit ordinary differential equations, differential-algebraic equations
65L80 Numerical methods for differential-algebraic equations
65M20 Method of lines for initial value and initial-boundary value problems involving PDEs
35Q53 KdV equations (Korteweg-de Vries equations)
65Y15 Packaged methods for numerical algorithms
Full Text: DOI


[1] Hairer, E.; Wanner, G., ()
[2] Brown, P.N.; Hindmarsh, A.C., Matrix-free methods for stiff systems of odes, SIAM J. numer. anal., 23, 3, 610-638, (1986) · Zbl 0615.65078
[3] Eisenstat, S.C.; Gursky, M.C.; Schultz, M.H.; Sherman, A.H., Yale sparse matrix package II: the nonsymmetric codes, Technical report 114, department of computer science, (1977), Yale University New Haven, CT · Zbl 0492.65012
[4] Duff, I.S., MA28—A set of Fortran subroutines for sparse unsymmetric linear equations, Technical report AERE-R8730, (1977), Harwell
[5] De Swart, J.J.B.; Blom, J.G., Experiences with sparse matrix solvers in parallel ODE software, Computers math. applic., 31, 9, 43-55, (1996) · Zbl 0859.65072
[6] Abdulla, T.J., Modified extended backward differentiation formulae for differential-algebraic equations with applications to time dependent pdes, Ph.D. thesis, (2000), Imperial College, University of London
[7] Lioen, W.M.; de Swart, J.J.B., Test set for IVP solvers, Technical report, (1996), CWI, Department of Mathematics Amsterdam, http://www.cwi.nl/cwi/projects/IVPtestset/
[8] Hindmarsh, A.C., ODEPACK:A systematized collection of ODE solvers, IMACS trans. on scientific computing, 1, 55-64, (1983)
[9] Cash, J.R., Modified extended backward differentiation formula for the numerical solution of stiff initial value problems in ODEs and daes, J. comp. app. math., 125, 117-130, (2000) · Zbl 0971.65063
[10] Thompson, S., The effect of partial pivoting in sparse ordinary differential equation solvers, Computers math. applic., 12A, 12, 1183-1191, (1986) · Zbl 0636.65074
[11] Madsen, N.K.; Sincovec, R.F., Pdecol, ACM toms, 5, 326-351, (1979) · Zbl 0426.35005
[12] Melgaard, D.K.; Sincovec, R.F., Pdetwo, ACM toms, 7, 126-135, (1981)
[13] Schiesser, W.E., The numerical method of lines. integration of partial differential equations, (1991), Academic Press San Diego, CA · Zbl 0763.65076
[14] Schiesser, W.E., Computational mathematics in engineering and applied science: ODEs, DAEs, and pdes, (1993), CRC Press Boca Raton, FL
[15] Finlayson, B.A., Solution of stiff equations resulting from partial differential equations, (), 124-139
[16] Hindmarsh, A.C., Avoiding BDF stability barriers in the MOL solution of advection-dominated problems, Applied numerical mathematics, 17, 311-318, (1995) · Zbl 0832.65102
[17] Whitham, G.B., Linear and nonlinear waves, (1974), Wiley New York · Zbl 0373.76001
[18] Schiesser, W.E., Method of lines solution of the Korteweg-de Vries equation, Computer math. applic., 28, 10-12, 147-154, (1994) · Zbl 0810.65125
[19] Saucez, P.; Vande Wouwer, A.; Schiesser, W.E., An adaptive method of lines solution of the kortewegde Vries equation, Computers math. applic., 35, 12, 13-25, (1998) · Zbl 0999.65098
[20] Fornberg, B., Fast generation of weights in finite difference formulas, (), 114
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.