A Time Adaptive Scheme Based On The Rush-Larsen Method For Solving Cardiac Cell Models

This paper presents a comparison between numerical schemes for solving electrophysiology models of cardiac myocytes and a new time adaptive scheme based on the Rush-Larsen method. The most recent models for the electrophysiology of cardiac cells are based on highly stiff nonlinear systems of ordinary differential equations. Several explicit numerical schemes for cardiac membrane models have been developed for increasing the computational performance in the simulation of the electrical activity of the heart. In this work we propose a time adaptive method based on the widely used Rush-Larsen method, which is known to have high numerical stability for solving models based on the Hodgkin-Huxley formulation for ionic channels. The new method could improve the computational performance up to six times compared to the performance of Rush-Larsen method.