Applying stiff integrators for ODEs and DDEs to problems with distributed delays
Abstract: There exist excellent codes for an efficient numerical treatment of stiff and differential-algebraic problems. Let us mention {\sc Radau5} which is based on the $3$-stage Radau IIA collocation method, and its extension to problems with discrete delays {\sc Radar5}. The aim of the present work is to present a technique that permits a direct application of these codes to problems having a right-hand side with an additional distributed delay term (which is a special case of an integro-differential equation). Models with distributed delays are of increasing importance in pharmacodynamics and pharmacokinetics for the study of the interaction between drugs and the body. The main idea is to approximate the distribution kernel of the integral term by a sum of exponential functions or by a quasi-polynomial expansion, and then to transform the distributed (integral) delay term into a set of ordinary differential equations. This set is typically stiff and, for some distribution kernels (e.g., Pareto distribution), it contains discrete delay terms with constant delay. The original equations augmented by this set of ordinary differential equations can have a very large dimension, and a careful treatment of the solution of the arising linear systems is necessary. The use of the codes {\sc Radau5} and {\sc Radar5} is illustrated at three examples (two test equations and one problem taken from pharmacodynamics). The driver programs for these examples are publicly available from the homepages of the authors.
- G. Beylkin and L. Monzón. Approximation by exponential sums revisited. Appl. Comput. Harmon. Anal., 28:131–149, 2010.
- D. Braess and W. Hackbusch. On the efficient computation of high-dimensional integrals and the approximation by exponential sums. In A. Kunoth R. DeVore, editor, Multiscale, Nonlinear and Adaptive Approximation, volume 56, pages 39 – 74. Springer Verlag, 2009.
- H. Brunner. Collocation methods for Volterra integral and related functional differential equations, volume 15 of Cambridge Monographs on Applied and Computational Mathematics. Cambridge University Press, Cambridge, 2004.
- The numerical solution of Volterra equations, volume 3 of CWI Monographs. North-Holland Publishing Co., Amsterdam, 1986.
- Transit and lifespan in neutrophil production: implications for drug intervention. Journal of Pharmacokinetics and Pharmacodynamics, 45(1):59 – 77, 2018.
- N. Guglielmi and E. Hairer. Implementing Radau IIA methods for stiff delay differential equations. Computing, 67(1):1–12, 2001.
- N. Guglielmi and E. Hairer. Computing breaking points in implicit delay differential equations. Adv. Comput. Math., 29(3):229–247, 2008.
- E. Hairer and P. Maass. Numerical methods for singular nonlinear integro-differential equations. Appl. Numer. Math., 3(3):243–256, 1987.
- E. Hairer and G. Wanner. Solving Ordinary Differential Equations II. Stiff and Differential-Algebraic Problems. Springer Series in Computational Mathematics 14. Springer-Verlag, Berlin, 2nd edition, 1996.
- W. Krzyzanski. Ordinary differential equation approximation of gamma distributed delay model. Journal of Pharmacokinetics and Pharmacodynamics, 46(1):53 – 63, 2019.
- Evaluation of performance of distributed delay model for chemotherapy-induced myelosuppression. Journal of Pharmacokinetics and Pharmacodynamics, 45(2):329 – 337, 2018.
- Adaptive, fast, and oblivious convolution in evolution equations with memory. SIAM J. Sci. Comput., 30(2):1015–1037, 2008.
- C. Lubich. Discretized fractional calculus. SIAM J. Math. Anal., 17(3):704–719, 1986.
- Fast and oblivious convolution quadrature. SIAM J. Sci. Comput., 28(2):421–438, 2006.
- The Exponentially Convergent Trapezoidal Rule. SIAM Review, 56(3):385–458, 2014.
- Comparison of the gamma-pareto convolution with conventional methods of characterising metformin pharmacokinetics in dogs. J. Pharmacokinet. Pharmacodyn., 47:19–45, 2020.
Paper Prompts
Sign up for free to create and run prompts on this paper using GPT-5.
Top Community Prompts
Collections
Sign up for free to add this paper to one or more collections.