Distributed computing for membrane-based modeling of action potential propagation

Action potential propagation simulations with physiologic membrane currents and macroscopic tissue dimensions are computationally expensive. The authors, therefore, analyzed distributed computing schemes to reduce execution time in workstation clusters by parallelizing solutions with message passing. Four schemes were considered in two-dimensional monodomain simulations with the Peeler-Reuter membrane equations. Parallel speedups measured with each scheme were compared to theoretical speedups, recognizing the relationship between speedup and code portions that executed serially. A data decomposition scheme based on total ionic current provided the best performance. Analysis of communication latencies in that scheme led to a load-balancing algorithm in which measured speedups at 89/spl plusmn/2% and 75/spl plusmn/8% of theoretical speedups were achieved in homogeneous and heterogeneous clusters of workstations. Speedups in this scheme with the Luo-Rudy dynamic membrane equations exceeded 3.0 with 8 distributed workstations. Cluster speedups were comparable to those measured during parallel execution on a shared memory machine.

[1]  C. Henriquez,et al.  Cardiac propagation simulation. , 1992, Critical reviews in biomedical engineering.

[2]  Jack Dongarra,et al.  Message-Passing Performance of Various Computers , 1995 .

[3]  A Garfinkel,et al.  Wave propagation in cardiac tissue and effects of intracellular calcium dynamics (computer simulation study). , 1998, Progress in biophysics and molecular biology.

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

[5]  Volker Strumpen,et al.  Efficient Parallel Computing in Distributed Workstation Environments , 1993, Parallel Comput..

[6]  M. Spach,et al.  Relating the Sodium Current and Conductance to the Shape of Transmembrane and Extracellular Potentials by Simulation: Effects of Propagation Boundaries , 1985, IEEE Transactions on Biomedical Engineering.

[7]  William H. Press,et al.  Numerical recipes in C. The art of scientific computing , 1987 .

[8]  C. Henriquez.  Simulating the electrical behavior of cardiac tissue using the bidomain model. , 1993, Critical reviews in biomedical engineering.

[9]  S Gao,et al.  Numerical analysis of electrical defibrillation. The parallel approach. , 1995, Journal of electrocardiology.

[10]  A. Y. Zomaya.  Parallel and Distributed Computing: The Scene, the Props, the Players , 1996 .

[11]  Denis Noble,et al.  Simulating cardiac sinus and atrial network dynamics on the Connection Machine , 1993 .

[12]  W. Press,et al.  Numerical Recipes: The Art of Scientific Computing , 1987 .

[13]  J W Moore,et al.  On numerical integration of the Hodgkin and Huxley equations for a membrane action potential. , 1974, Journal of theoretical biology.

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

[15]  N V Thakor,et al.  A Massively Parallel Computer Model of Propagation Through a Two‐Dimensional Cardiac Syncytium , 1991, Pacing and clinical electrophysiology : PACE.

[16]  C. Henriquez,et al.  Bipolar stimulation of a three-dimensional bidomain incorporating rotational anisotropy , 1998, IEEE Transactions on Biomedical Engineering.

[17]  R. Barr,et al.  Computer simulations of activation in an anatomically based model of the human ventricular conduction system , 1991, IEEE Transactions on Biomedical Engineering.

[18]  Y Rudy,et al.  Two components of the delayed rectifier K+ current in ventricular myocytes of the guinea pig type. Theoretical formulation and their role in repolarization. , 1995, Circulation research.

[19]  Michael D. Lesh,et al.  Cellular Uncoupling Can Unmask Dispersion of Action Potential Duration in Ventricular Myocardium A Computer Modeling Study , 1989, Circulation research.

[20]  K.T. Ng,et al.  A new three-dimensional finite-difference bidomain formulation for inhomogeneous anisotropic cardiac tissues , 1998, IEEE Transactions on Biomedical Engineering.

[21]  David J. Jackson,et al.  A Simple Yet Effective Load Balancing Extension to the PVM Software System , 1997, Parallel Comput..

[22]  C. Luo,et al.  A dynamic model of the cardiac ventricular action potential. I. Simulations of ionic currents and concentration changes. , 1994, Circulation research.

[23]  R. Barr,et al.  Propagation of excitation in idealized anisotropic two-dimensional tissue. , 1984, Biophysical journal.

[24]  Weilun Quan,et al.  Efficient integration of a realistic two-dimensional cardiac tissue model by domain decomposition , 1998, IEEE Transactions on Biomedical Engineering.

[25]  R. Winslow,et al.  Cardiac Ca2+ dynamics: the roles of ryanodine receptor adaptation and sarcoplasmic reticulum load. , 1998, Biophysical journal.

[26]  Jean-Pierre Drouhard,et al.  A Simulation Study of the Ventricular Myocardial Action Potential , 1982, IEEE Transactions on Biomedical Engineering.

[27]  D. Rose,et al.  Linear algebraic transformations of the bidomain equations: implications for numerical methods. , 1994, Mathematical biosciences.

[28]  H. I. Saleheen,et al.  Three-dimensional finite-difference bidomain modeling of homogeneous cardiac tissue on a data-parallel computer , 1997, IEEE Transactions on Biomedical Engineering.

[29]  Y. Rudy,et al.  Unidirectional block and reentry of cardiac excitation: a model study. , 1990, Circulation research.

[30]  B. Victorri,et al.  Numerical integration in the reconstruction of cardiac action potentials using Hodgkin-Huxley-type models. , 1985, Computers and biomedical research, an international journal.

[31]  D DiFrancesco,et al.  A model of cardiac electrical activity incorporating ionic pumps and concentration changes. , 1985, Philosophical transactions of the Royal Society of London. Series B, Biological sciences.

[32]  J E Saffitz,et al.  The Molecular Basis of Anisotropy: Role of Gap Junctions , 1995, Journal of cardiovascular electrophysiology.

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

[34]  R. W. Joyner,et al.  Effects of the Discrete Pattern of Electrical Coupling on Propagation through an Electrical Syncytium , 1982, Circulation research.