Cardiac simulation on multi-GPU platform

The cardiac bidomain model is a popular approach to study electrical behavior of tissues and simulate interactions between the cells by solving partial differential equations. The iterative and data parallel model is an ideal match for the parallel architecture of Graphic Processing Units (GPUs). In this study, we evaluate the effectiveness of architecture-specific optimizations and fine grained parallelization strategies, completely port the model to GPU, and evaluate the performance of single-GPU and multi-GPU implementations. Simulating one action potential duration (350 msec real time) for a 256×256×256 tissue takes 453 hours on a high-end general purpose processor, while it takes 664 seconds on a four-GPU based system including the communication and data transfer overhead. This drastic improvement (a factor of 2460×) will allow clinicians to extend the time-scale of simulations from milliseconds to seconds and minutes; and evaluate hypotheses in a shorter amount of time that was not feasible previously.

[1]  B. Roth,et al.  Action potential propagation in a thick strand of cardiac muscle. , 1991, Circulation research.

[2]  D. Noble,et al.  Reconstruction of the electrical activity of cardiac Purkinje fibres. , 1975, The Journal of physiology.

[3]  Daisuke Sato,et al.  Acceleration of cardiac tissue simulation with graphic processing units , 2009, Medical & Biological Engineering & Computing.

[4]  A. Garfinkel,et al.  An advanced algorithm for solving partial differential equation in cardiac conduction , 1999, IEEE Transactions on Biomedical Engineering.

[5]  C. Nicholson.  Electric current flow in excitable cells J. J. B. Jack, D. Noble &R. W. Tsien Clarendon Press, Oxford (1975). 502 pp., £18.00 , 1976, Neuroscience.

[6]  K. T. ten Tusscher,et al.  Alternans and spiral breakup in a human ventricular tissue model. , 2006, American journal of physiology. Heart and circulatory physiology.

[7]  David Kaeli,et al.  Proceedings of 2nd Workshop on General Purpose Processing on Graphics Processing Units , 2009 .

[8]  G. W. Beeler,et al.  Reconstruction of the action potential of ventricular myocardial fibres , 1977, The Journal of physiology.

[9]  G. Smith,et al.  Numerical Solution of Partial Differential Equations: Finite Difference Methods , 1978 .

[10]  A. Noma,et al.  Reconstruction of sino-atrial node pacemaker potential based on the voltage clamp experiments. , 1980, The Japanese journal of physiology.

[11]  Gernot Plank,et al.  Near-real-time simulations of biolelectric activity in small mammalian hearts using graphical processing units , 2009, 2009 Annual International Conference of the IEEE Engineering in Medicine and Biology Society.

[12]  D. Noble.  A modification of the Hodgkin—Huxley equations applicable to Purkinje fibre action and pacemaker potentials , 1962, The Journal of physiology.

[13]  K.H.W.J. ten Tusscher,et al.  Comments on 'A model for human ventricular tissue' : reply , 2005 .

[14]  Paulius Micikevicius,et al.  3D finite difference computation on GPUs using CUDA , 2009, GPGPU-2.

[15]  C. Luo,et al.  A model of the ventricular cardiac action potential. Depolarization, repolarization, and their interaction. , 1991, Circulation research.