Abstract
In this paper, we survey selected software packages for the numerical solution of boundary value ODEs (BVODEs), time-dependent PDEs in one spatial dimension (1DPDEs), and initial value ODEs (IVODEs). A unifying theme of this paper is our focus on software packages for these problem classes that compute error-controlled, continuous numerical solutions. A continuous numerical solution can be accessed by the user at any point in the domain. We focus on error-control software; this means that the software adapts the computation until it obtains a continuous approximate solution with a corresponding error estimate that satisfies the user tolerance. The second section of the paper will provide an overview of recent work on the development of COLNEWSC, an updated version of the widely used collocation BVODE solver, COLNEW, that returns an error-controlled continuous approximate solution based on the use of a superconvergent interpolant to the underlying collocation solution. The third section of the paper gives a brief review of recent work on the development of a new 1DPDE solver, BACOLIKR, that provides time- and space-dependent event detection for an error-controlled continuous numerical solution. In the fourth section of the paper, we briefly review the state of the art in IVODE software for the computation of error-controlled continuous numerical solutions.
Similar content being viewed by others
References
Test Set for BVP Solvers. https://archimede.uniba.it/~bvpsolvers/testsetbvpsolvers/
Agowun, H., Muir, P.H.: Performance analysis of ODE solvers on Covid-19 models with discontinuities. Saint Mary’s University, Dept. of Mathematics and Computing Science Technical Report Series, Technical Report 2022_001 (2022). http://cs.smu.ca/tech_reports
Agowun, H., Muir, P.H.: Low cost multistep interpolants for Runge Kutta methods. Work in progress, (2024)
Allen, S.M., Cahn, J.W.: Ground state structures in ordered binary alloys with second neighbor interations. Acta Metallurgica 20, 423–433 (1972)
Arsenault, T., Smith, T., Muir, P.H.: Superconvergent interpolants for efficient spatial error estimation in 1D PDE collocation solvers. Can. Appl. Math. Q. 17, 409–431 (2009)
Arsenault, T., Smith, T., Muir, P.H., Pew, J.: Asymptotically correct interpolation-based spatial error estimation for 1D PDE solvers. Can. Appl. Math. Q. 20, 307–328 (2012)
Ascher, U., Christiansen, J., Russell, R.D.: Algorithm 569: COLSYS: Collocation software for boundary-value ODEs [d2]. ACM Trans. Math. Softw. (TOMS) 7(2), 223–229 (1981)
Ascher, U.M., Mattheij, R.M.M. Russell, R.D.: Numerical Solution of Boundary Value Problems for Ordinary Dfferential Euations, vol. 13 of Classics in Applied Mathematics. Society for Industrial and Applied Mathematics (SIAM), Philadelphia, PA, (1995)
Bader, G., Ascher, U.M.: A new basis implementation for a mixed order boundary value ODE solver. SIAM J. Sci. Stat. Comput. 8(4), 483–500 (1987)
Boehme, T.R., Onder, C.H., Guzzella, L.: Code-generator-based software package for defining and solving one-dimensional, dynamic, catalytic reactor models. Comput. Chem. Eng. 32(10), 2445–2454 (2008)
Boisvert, J.J., Muir, P.H., Spiteri, R.J.: A Runge-Kutta BVODE solver with global error and defect control. ACM Transactions on Mathematical Software 39(2), Art. 11, 22 (2013)
Bonaventura, Z., Trunec, D., Meško, M., Vašina, P., Kudrle, V.: Theoretical study of pulsed microwave discharge in nitrogen. Plasma Sources Sci. Technol. 14(4), 751 (2005)
Brenan, K.E., Campbell, S.L., Petzold, L.R.: Numerical Solution of Initial-Value Problems in Differential-Algebraic Equations, vol. 14 of Classics in Applied Mathematics. Society for Industrial and Applied Mathematics (SIAM), Philadelphia, PA, (1996)
Brown, P.N., Hindmarsh, A.C., Petzold, L.R.: Using Krylov methods in the solution of large-scale differential-algebraic systems. SIAM J. Sci. Comput. 15, 1467–1488 (1994)
Brown, P.N., Hindmarsh, A.C., Petzold, L.R.: Consistent initial condition calculation for differential-algebraic systems. SIAM J. Sci. Comput. 19, 1495–1512 (1998)
Burrage, K., Chipman, F.H., Muir, P.H.: Order results for mono-implicit Runge-Kutta methods. SIAM J. Numer. Anal. 31(3), 876–891 (1994)
Cerutti, J.H., Parter, S.V.: Collocation methods for parabolic partial differential equations in one space dimension. Numerische Mathematik 26(3), 227–254 (1976)
Collodel, L.G., Doneva, D.D.: Solitonic Boson Stars: Numerical solutions beyond the thin-wall approximation. Phys. Rev. D 106(8), 084057 (2022)
de Boor, C.: A Practical Guide to Splines, vol. 27 of Applied Mathematical Sciences. Springer-Verlag, New York, revised edition, (2001)
Enright, W.H., Hayes, W.B.: Robust and reliable defect control for Runge-Kutta methods. ACM Transactions on Mathematical Software 33(1), 1-es (2007)
Enright, W.H., Muir, P.H.: Superconvergent interpolants for the collocation solution of boundary value ordinary differential equations. SIAM J. Sci. Comput. 21(1), 227–254 (1999)
Falini, A., Mazzia, F., Sestini, A.: Hermite–Birkhoff spline Quasi-Interpolation with application as dense output for Gauss–Legendre and Gauss–Lobatto Runge–Kutta schemes. Applied Numerical Mathematics, (2023)
Gear, C.W., Osterby, O.: Solving ordinary differential equations with discontinuities. ACM Trans. Math. Softw. 10(1), 23–44 (1984)
Green, K.R., Spiteri, R.J.: Extended BACOLI: Solving one-dimensional multiscale parabolic PDE systems with error control. ACM Trans. Math. Softw. 45(1), 1–19 (2019)
Mazzia, F., Cash, J.R.: A Fortran test set for boundary value problem solvers. In AIP Conference Proceedings, vol. 1648. AIP Publishing (2015)
Mazzia, F., Cash, J.R., Soetaert, K.: A test set for stiff initial value problem solvers in the open source software r: Package detestset. J. Comput. Appl. Math. 236(16), 4119–4131 (2012)
Mazzia, F., Cash, J.R., Soetaert, K.: Solving boundary value problems in the open source software R: Package bvpSolve. Opusc. Math. 34(2), 387–403 (2014)
Mazzia, F., Sestini, A.: The BS class of hermite spline quasi-interpolants on nonuniform knot distributions. BIT Numer. Math. 49, 611–628 (2009)
Muir, P.H., Owren, B.: Order barriers and characterizations for continuous mono-implicit Runge-Kutta schemes. Math. Comput. 61(204), 675–699 (1993)
Petit, N., Sciarretta, A.: Optimal drive of electric vehicles using an inversion-based trajectory generation approach. IFAC Proceedings Volumes 44(1), 14519–14526 (2011)
Pew, J., Li, Z., Muir, P.H.: Algorithm 962: BACOLI: B-spline adaptive collocation software for PDEs with interpolation-based spatial error control. ACM Transactions on Mathematical Software 42(3), 25:1-25:17 (2016)
Pew, J., Tannahill, C., Muir, P.H.: Error control B-spline Gaussian collocation PDE software with time-space event detection. Saint Mary’s University, Dept. of Mathematics and Computing Science Technical Report Series, Technical Report 2020_001 (2020) http://cs.smu.ca/tech_reports
Shampine, L.F.: Efficient Runge-Kutta pairs with local error control equivalent to residual control. Private communication, (2021)
Shampine, L.F.: Error control and ode23. Private communication, (2021)
Soetaert, K., Petzoldt, T., Setzer, R.W.: Solving differential equations in r: package desolve. J. Stat. Softw. 33, 1–25 (2010)
Tannahill, C., Muir, P.: bacoli_py - python package for the error controlled numerical solution of 1D time-dependent PDEs. In International Conference on Applied Mathematics, Modeling and Computational Science, p. 289–300. Springer, (2021)
Trunec, D., Bonaventura, Z., Nečas, D.: Solution of time-dependent Boltzmann equation for electrons in non-thermal plasma. J. Phys. D: Appl. Phys. 39(12), 2544 (2006)
Tsitouras, Ch.: Runge-Kutta interpolants for high precision computations. Numer. Algorithms 44, 291–307 (2007)
Vallino, J.J., Huber, J.A.: Using maximum entropy production to describe microbial biogeochemistry over time and space in a meromictic pond. Front. Environ. Sci. 6, 100 (2018)
Wang, R., Keast, P., Muir, P.H.: BACOL: B-spline Adaptive COLlocation software for 1D parabolic PDEs. ACM Trans. Math. Softw. 30(4), 454–470 (2004)
Wang, R., Keast, P., Muir, P.H.: A high-order global spatially adaptive collocation method for 1-D parabolic PDEs. Appl. Numer. Math. 50(2), 239–260 (2004)
Zong, L., Emami, M.R.: Concurrent base-arm control of space manipulators with optimal rendezvous trajectory. Aerosp. Sci. Technol. 100, 105822 (2020)
Acknowledgements
The authors would like to thank the referees for many helpful suggestions.
Funding
The second author is funded by the Natural Sciences and Engineering Research Council of Canada and Saint Mary’s University.
Author information
Authors and Affiliations
Contributions
Both authors contributed equally to the manuscript.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Springer Nature or its licensor (e.g. a society or other partner) holds exclusive rights to this article under a publishing agreement with the author(s) or other rightsholder(s); author self-archiving of the accepted manuscript version of this article is solely governed by the terms of such publishing agreement and applicable law.
About this article
Cite this article
Adams, M., Muir, P. Differential equation software for the computation of error-controlled continuous approximate solutions. Numer Algor 96, 1021–1044 (2024). https://doi.org/10.1007/s11075-024-01784-1
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s11075-024-01784-1
Keywords
- Initial value ordinary differential equations
- Boundary value ordinary differential equations
- Partial differential equations
- Continuous numerical solutions
- Collocation
- Adaptive methods
- Error estimation
- Error control
- Efficiency
- Reliability