Abstract
We present a computational modeling framework for data-driven simulations and analysis of infectious disease spread in large populations. For the purpose of efficient simulations, we devise a parallel solution algorithm targeting multi-socket shared-memory architectures. The model integrates infectious dynamics as continuous-time Markov chains and available data such as animal movements or aging are incorporated as externally defined events. To bring out parallelism and accelerate the computations, we decompose the spatial domain and optimize cross-boundary communication using dependency-aware task scheduling. Using registered livestock data at a high spatiotemporal resolution, we demonstrate that our approach not only is resilient to varying model configurations but also scales on all physical cores at realistic workloads. Finally, we show that these very features enable the solution of inverse problems on national scales.
References
|
Andreev, K, Räcke, H (2004) Balanced graph partitioning. In: Proceedings of the 16th ACM SPAA, SPAA ‘04. New York, NY, USA, 27–30 June 2004, pp. 120–124. New York: ACM. Google Scholar | Crossref | |
|
Anonymous (2000) Regulation (EC) No 1760/2000 of the European Parliament and of the Council of 17 July 2000 establishing a system for the identification and registration of bovine animals and regarding the labelling of beef and beef products. Official Journal of the European Union L 204: 1–10. Google Scholar | |
|
Anonymous (2004) Commission Regulation (EC) No 911/2004 of 29 April 2004 implementing Regulation (EC) No 1760/2000 of the European Parliament and of the Council as regards eartags, passports and holding registers. Official Journal of the European Union L 163: 65–70. Google Scholar | |
|
Augonnet, C, Thibault, S, Namyst, R. (2011) StarPU: a unified platform for task scheduling on heterogeneous multicore architectures. Concurrency and Computation Practice and Experience 23(2): 187–198. Google Scholar | Crossref | ISI | |
|
Barrett, C, Bisset, K, Eubank, S. (2008) EpiSimdemics: an efficient algorithm for simulating the spread of infectious disease over large realistic social networks. In: Proceedings of Supercomputing 2008, Austin, TX, 15–21 November 2008, pp. 1–12. Piscatway, NJ: IEEE. Google Scholar | Crossref | |
|
Berry, LA, Elwasif, W, Reynolds-Barredo, JM. (2012) Event-based parareal: a data-flow based implementation of parareal. Journal of Computational Physics 231(17): 5945–5954. Google Scholar | Crossref | ISI | |
|
Brooks-Pollock, E, de Jong, M, Keeling, M. (2015) Eight challenges in modelling infectious livestock diseases. Epidemics 10: 1–5. Google Scholar | Crossref | Medline | ISI | |
|
Büttner, K, Krieter, J, Traulsen, A. (2016) Epidemic spreading in an animal trade network - comparison of distance-based and network-based control measures. Transboundary and Emerging Diseases 63: e122–e134. Google Scholar | Crossref | ISI | |
|
Carothers, CD, Perumalla, KS, Fujimoto, RM (1999) Efficient optimistic parallel simulations using reverse computation. ACM Transactions on Modeling and Computer Simulation 9(3): 224–253. Google Scholar | Crossref | |
|
Cassandras, CG, Lafortune, S (2008) Systems and models. In: Introduction to Discrete Event Systems. Boston, MA: Springer, pp. 1–51. Google Scholar | Crossref | |
|
Cray, WC, Moon, HW (1995) Experimental infection of calves and adult cattle with Escherichia coli O157: H7. Applied and Environmental Microbiology 61(4): 1586–1590. Google Scholar | Medline | ISI | |
|
Davis, MA, Rice, DH, Sheng, H. (2006) Comparison of cultures from rectoanal-junction mucosal swabs and feces for detection of Escherichia coli O157 in dairy heifers. Applied and Environmental Microbiology 72(5): 3766–3770. Google Scholar | Crossref | Medline | ISI | |
|
Drawert, B, Engblom, S, Hellander, A (2012) URDME: a modular framework for stochastic simulation of reaction-transport processes in complex geometries. BMCSystems Biology 6(1): 76. Google Scholar | |
|
Drozdowski, M (2009) Parallel tasks. In: Sammes, AJ (ed.) Scheduling for Parallel Processing (Computer Communications and Networks). London: Springer, pp. 87–208. Google Scholar | Crossref | |
|
Duran, A, Ayguadé, E, Badia, RM. (2011) OmpSs: a proposal for programming heterogeneous multi-core architectures. Parallel Processing Letters 21(02): 173–193. Google Scholar | Crossref | |
|
Engblom, S (2015) Strong convergence for split-step methods in stochastic jump kinetics. SIAM Journal on Numerical Analysis 53(6): 2655–2676. Google Scholar | Crossref | ISI | |
|
Engblom, S, Ferm, L, Hellander, A. (2009) Simulation of stochastic reaction-diffusion processes on unstructured meshes. SIAM Journal on Scientific Computing 31(3): 1774–1797. Google Scholar | Crossref | ISI | |
|
Fujimoto, RM (1990) Parallel discrete event simulation. Communications of the ACM 33(10): 30–53. Google Scholar | |
|
Fujimoto, RM (1999) Parallel and Distribution Simulation Systems. 1st ed.New York: John Wiley & Sons, Inc. Google Scholar | |
|
Gerasoulis, A, Yang, T (1993) On the granularity and clustering of directed acyclic task graphs. IEEE Transactions on Parallel and Distributed Systems 4(6): 686–701. Google Scholar | Crossref | ISI | |
|
Gibson, MA, Bruck, J (2000) Efficient exact stochastic simulation of chemical systems with many species and many channels. Journal of Physical Chemistry A 104(9): 1876–1889. Google Scholar | Crossref | ISI | |
|
Gillespie, DT (1977) Exact stochastic simulation of coupled chemical reactions. Journal of Physical Chemistry 81(25): 2340–2361. Google Scholar | Crossref | ISI | |
|
Greenwood, PE, Gordillo, LF (2009) Stochastic epidemic modeling. In: Chowell, G, Hyman, JM, Bettencourt, LMA. (eds) Mathematical and Statistical Estimation Approaches in Epidemiology. Amsterdam, Netherlands: Springer, pp. 31–52. Google Scholar | Crossref | |
|
Haidar, A, Tomov, S, Dongarra, J. (2014) A novel hybrid CPU-GPU generalized eigensolver for electronic structure calculations based on fine-grained memory aware tasks. International Journal of High Performance Computing Applications 28(2): 196–209. Google Scholar | SAGE Journals | ISI | |
|
Hancock, D, Besser, T, Lejeune, J. (2001) The control of VTEC in the animal reservoir. International Journal of Food Microbiology 66(1-2): 71–78. Google Scholar | Crossref | Medline | ISI | |
|
Harvey, N, Reeves, A, Schoenbaum, MA. (2007) The North American animal disease spread model: a simulation model to assist decision making in evaluating animal disease incursions. Preventive Veterinary Medicine 82(3–4): 176–197. Google Scholar | Crossref | Medline | ISI | |
|
Hasonova, L, Pavlik, I (2006) Economic impact of paratuberculosis in dairy cattle herds: a review. Veterinarni Medicina-Czech 51(5): 193–211. Google Scholar | Crossref | ISI | |
|
Heidelberger, P, Nicol, D (1993) Conservative parallel simulation of continuous time Markov chains using uniformization. IEEETransactions on Parallel and Distributed Systems 4(8): 906–921. Google Scholar | Crossref | ISI | |
|
Hooke, R, Jeeves, TA (1961) Direct search: solution of numerical and statistical problems. Journal of the ACM 8(2): 212–229. Google Scholar | Crossref | |
|
Jefferson, DR (1985) Virtual time. ACM Transactions on Programming Languages and Systems 7(3): 404–425. Google Scholar | Crossref | ISI | |
|
Karmali, MA, Petric, M, Lim, C. (1983a) Escherichia coli cytotoxin, haemolytic-uraemic syndrome, and haemorrhagic colitis. Lancet 2(8362): 1299–1300. Google Scholar | Crossref | Medline | ISI | |
|
Karmali, MA, Steele, BT, Petric, M. (1983b) Sporadic cases of haemolytic-uraemic syndrome associated with faecal cytotoxin and cytotoxin-producing escherichia coli in stools. Lancet 1(8325): 619–620. Google Scholar | Crossref | Medline | ISI | |
|
Keeling, MJ (2005) Models of foot-and-mouth disease. Proceedings of the Royal Society B 272(1569): 1195–1202. Google Scholar | Crossref | Medline | ISI | |
|
Kermack, WO, McKendrick, AG (1927) A contribution to the mathematical theory of epidemics. Proceedings of the Royal Society A 115: 700–721. Google Scholar | Crossref | |
|
Kiefer, J (1953) Sequential minimax search for a maximum. Proceedings of the American Mathematical Society 4(3): 502–506. Google Scholar | Crossref | ISI | |
|
Knight-Jones, TJD, Rushton, J (2013) The economic impacts of foot and mouth disease – what are they, how big are they and where do they occur?Preventive Veterinary Medicine 112(3-4): 161–173. Google Scholar | Crossref | Medline | ISI | |
|
Leijen, D, Schulte, W, Burckhardt, S (2009) The design of a task parallel library. In: Proceedings of the 24th ACM SIGPLAN Conference on object oriented programming systems languages and applications, OOPSLA. New York, NY, USA, pp. 227–242. New York: ACM. Google Scholar | Crossref | |
|
Masuda, N, Holme, P (2013) Predicting and controlling infectious disease epidemics using temporal networks. F1000Prime Reports 5: 6. Google Scholar | Crossref | Medline | |
|
Meng, Q, Berzins, M (2014) Scalable large-scale fluid–structure interaction solvers in the uintah framework via hybrid task-based parallelism algorithms. Concurrency and Computation Practice and Experience 26(7): 1388–1407. Google Scholar | Crossref | ISI | |
|
Nicol, D, Liu, J (2002) Composite synchronization in parallel discrete-event simulation. IEEE Transactions on Parallel and Distributed Systems 13(5): 433–446. Google Scholar | Crossref | ISI | |
|
OpenMP Architecture Review Board (2013) OpenMP 4.0 Application Program Interface, p. 117, Lines 14–15. Google Scholar | |
|
Pellis, L, Ball, F, Bansal, S. (2015) Eight challenges for network epidemic models. Epidemics 10: 58–62. Google Scholar | Crossref | Medline | ISI | |
|
Perez, J, Badia, R, Labarta, J (2008) A dependency-aware task-based programming environment for multi-core architectures. In: 2008 IEEE International Conference on Cluster Computing, Tsukuba, 29 September–1 October 2008, pp. 142–151. IEEE. Google Scholar | Crossref | |
|
Riley, LW, Remis, RS, Helgerson, SD. (1983) Hemorrhagic colitis associated with a rare Escherichia coli serotype. New England Journal of Medicine 308(12): 681–685. Google Scholar | Crossref | Medline | ISI | |
|
Shirley, MD, Rushton, SP (2005) The impacts of network topology on disease spread. Ecological Complexity 2(3): 287–299. Google Scholar | Crossref | ISI | |
|
Stevenson, MA, Sanson, RL, Stern, MW. (2013) InterSpread plus: a spatial and stochastic simulation model of disease in animal populations. Preventive Veterinary Medicine 109(1–2): 10–24. Google Scholar | Crossref | Medline | ISI | |
|
Subhlok, J, Stichnoth, JM, O’Hallaron, DR. (1993) Exploiting task and data parallelism on a multicomputer. In: Proceedings of the 4th PPOPP, New York, NY, USA: ACM, pp. 13–22. Google Scholar | Crossref | |
|
Tillenius, M (2015) SuperGlue: a shared memory framework using data version-ing for dependency-aware task-based parallelization. SIAM Journal on Scientific Computing 37(6): C617–C642. Google Scholar | Crossref | ISI | |
|
Widgren, S, Engblom, S, Bauer, P. (2016) Data-driven network modeling of VTEC O157 transmission in Swedish cattle using complete population movement data. Under revision. Google Scholar | |
|
Willeberg, P, Grubbe, T, Weber, S. (2011) The world organisation for animal health and epidemiological modelling: background and objectives. Revue Scientifique et Technique 30(2): 391–405. Google Scholar | Crossref | |
|
Xiao, Z, Unger, B, Simmonds, R. (1999) Scheduling critical channels in conservative parallel discrete event simulation. In: Proceedings of the 13th Workshop of PADS, Atlanta, GA, 1–4 May, 1999, pp. 20–28. Washington, DC: IEEE. Google Scholar | Crossref | |
|
Yeom, JS, Bhatele, A, Bisset, K. (2014) Overcoming the scalability challenges of epidemic simulations on Blue Waters. In: Proceedings of the 28th IEEE International Parallel and Distributed Processing Symposium, Phoenix, AZ, pp. 755–764. Google Scholar | Crossref | |
|
Zafari, A, Tillenius, M, Larsson, E (2012) Programming models based on data versioning for dependency-aware task-based parallelisation. In: Proceedings of the 15th International Conference on Computational Science and Engineering, Nicosia, Cyprus, 5–7 December 2012, pp. 275–280. Los Alamitos: IEEE Computer Society. Google Scholar | Crossref |
