ABSTRACT
Simulating the human heart is a challenging problem, with simulations being very time consuming---some can take days to compute even on high performance computing resources. There is considerable interest in optimisation techniques, with a view to making whole-heart simulations tractable. Reliability of heart model simulations is also of great concern, particularly considering clinical applications. Simulation software should be easily testable (against empirical data) and maintainable, which is often not the case with extensively hand-optimised software. Automating any optimisations will greatly improve this situation. This paper presents a framework for automatically optimising cardiac ionic cell models. An abstract format for such models, CellML [9], has been developed at Auckland University and is gaining in popularity. We utilise this format and investigate robust transformations of models that lead to reduced simulation times. In particular, we demonstrate that partial evaluation [13] is a promising technique for this purpose, and that it combines well with a lookup table technique, commonly used in cardiac modelling, which we have automated. In our tests, the technique of partial evaluation gives a speedup of 1.2 times. Further, applying such transformations prior to the lookup table optimisation results in additional speedups from the latter technique. When both optimisations are used we obtain nearly a 5-fold speedup, compared with a speedup of 1.7 times when lookup tables alone are used.
- P. H. Andersen. Partial evaluation applied to ray tracing. Technical Report 95/2, DIKU, University of Copenhagen, January 1995.Google Scholar
- A. A. Berlin. Partial evaluation applied to numerical computation. In 1990 ACM Conference on Lisp and Functional Programming, Nice, France, pages 139--150, 1990. Google ScholarDigital Library
- A. A. Berlin and D. Weise. Compiling scientific code using partial evaluation. IEEE Computer, 23(12):25--37, December 1990. Google ScholarDigital Library
- H. Boehm and R. Cartwright. Exact real arithmetic: Formulating real numbers as functions. In D. A. Turner, editor, Research Topics in Functional Programming, chapter 3. Addison-Wesley, 1990. Google ScholarDigital Library
- F. Dexter, G. M. Saidel, M. N. Levy, and Y. Rudy. Mathematical model of dependence of heart rate on tissue concentration of acetylcholine. Am J Physiol Heart Circ Physiol, 256(2):H520--526, 1989.Google ScholarCross Ref
- A. Garny, P. Kohl, P. J. Hunter, M. R. Boyett, and D. Noble. One-Dimensional Rabbit Sinoatrial Node Models: Benefits and Limitations. J Cardiovasc Electrophysiol, 14(s10):S121--S132, 2003.Google Scholar
- A. Garny, P. Kohl, and D. Noble. Cellular Open Resource (COR): a public CellML based environment for modelling biological function. Int J Bif Chaos, 13(12):3579--3590, 2003.Google ScholarCross Ref
- E. Hairer, S. P. Norsett, and G. Wanner. Solving Ordinary Differential Equations I. Nonstiff Problems. Springer-Verlag, 2 edition, 1993. Google ScholarDigital Library
- W. J. Hedley, M. R. Nelson, D. P. Bullivant, and P. F. Nielsen. A short introduction to CellML. Philosophical Transactions: Mathematical, Physical and Engineering Sciences, 359(1783):1073--1089, June 2001.Google ScholarCross Ref
- A. L. Hodgkin and A. F. Huxley. A quantitative description of membrane current and its application to conduction and excitation in the nerve. J Physiol, 117:500--544, 1952.Google ScholarCross Ref
- J. Hughes. Generalising Monads to Arrows. Science of Computer Programming, 37:67--111, May 2000. Google ScholarDigital Library
- P. J. Hunter, A. J. Pullen, and B. H. Smaill. Modelling total heart function. Annu. Rev. Biomed. Eng., 5:147--177, 2003.Google ScholarCross Ref
- N. D. Jones, C. K. Gomard, and P. Sestoft. Partial Evaluation and Automatic Program Generation. Prentice Hall, 1993. Google ScholarDigital Library
- L. J. Leon and F. A. Roberge. A new cable model based on Green's theorem. Ann Biomed Eng, 18(1):1--17, 1990.Google ScholarCross Ref
- L. M. Loew and J. C. Schaff. The Virtual Cell: a software environment for computational cell biology. Trends Biotechnol, 19:401--406, 2001.Google ScholarCross Ref
- C. Luo and Y. Rudy. Dynamic Model of the Cardiac Ventricular Action Potential---Simulations of Ionic Currents and Concentration Changes. Circulation Research, 74:1071--1097, 1994.Google ScholarCross Ref
- D. McNamee, J. Walpole, C. Cowan, C. Pu, C. Krasic, P. Wagle, C. Consel, G. Muller, and R. Marlet. Specialization Tools and Techniques for Systematic Optimization of Systems Software. ACM Transactions on Computer Systems, 19(2):217--251, May 2001. Google ScholarDigital Library
- D. Michie. Memo functions and machine learning. Nature, 218:19--22, April 1968.Google ScholarCross Ref
- S. S. Muchnick. Advanced Compiler Design and Implementation. Morgan Kaufmann, 1997. Google ScholarDigital Library
- D. Noble. Cardiac action and pacemaker potentials based on the hodgkin-huxley equations. Nature, 188:495--497, 1960.Google ScholarCross Ref
- D. Noble. Modelling the heart: insights, failures and progress. BioEssays, 24:1155--1163, 2002.Google ScholarCross Ref
- E. Süli and D. Mayers. An Introduction to Numerical Analysis. Cambridge University Press, 2003.Google ScholarCross Ref
- P. Thiemann. Towards partial evaluation of full scheme. In G. Kiczales, editor, Reflection'96, pages 95--106, April 1996.Google Scholar
- L. Tung. A bidomain model for describing ischemic myocardial D.C. potentials. PhD thesis, Massachusettes Institute of Technology, Cambridge, MA, 1978.Google Scholar
- E. J. Vigmond, M. Hughes, G. Plank, and L. J. Leon. Computational tools for modeling electrical activity in cardiac tissue. Journal of Electrocardiology, 36:69--74, Dec 2003.Google ScholarCross Ref
- D. Weise and S. Seligman. Accelerating Object-Oriented Simulation via Automatic Program Specialization. Technical Report CSL-TR-92-519, Stanford University, April 1992.Google Scholar
Index Terms
- On the application of partial evaluation to the optimisation of cardiac electrophysiological simulations
Recommendations
Automatic Segmentation of Cardiac CTs - Personalized Atrial Models Augmented with Electrophysiological Structures
FIMH 2011: Proceedings of the 6th International Conference on Functional Imaging and Modeling of the Heart - Volume 6666Electrophysiological simulations of the atria could improve diagnosis and treatment of cardiac arrhythmia, like atrial fibrillation or flutter. For this purpose, a precise segmentation of both atria is needed. However, the atrial epicardium and the ...
Automatic segmentation of cardiac CTs - personalized atrial models augmented with electrophysiological structures
FIMH'11: Proceedings of the 6th international conference on Functional imaging and modeling of the heartElectrophysiological simulations of the atria could improve diagnosis and treatment of cardiac arrhythmia, like atrial fibrillation or flutter. For this purpose, a precise segmentation of both atria is needed. However, the atrial epicardium and the ...
Optimization of cardiac resynchronization therapy based on a cardiac electromechanics-perfusion computational model
AbstractCardiac resynchronization therapy (CRT) is an established treatment for left bundle branch block (LBBB) resulting in mechanical dyssynchrony. Approximately 1/3 of patients with CRT, however, are non-responders. To understand factors ...
Highlights- A novel electromechanics-perfusion computational model based on animal-specific geometry and measurements is developed.
Comments