Computational modelling of cardiac electrophysiology: explanation of the variability of results from different numerical solvers

A recent verification study compared 11 large-scale cardiac electrophysiology solvers on an unambiguously defined common problem. An unexpected amount of variation was observed between the codes, including significant error in conduction velocity in the majority of the codes at certain spatial resolutions. In particular, the results of the six finite element codes varied considerably despite each using the same order of interpolation. In this present study, we compare various algorithms for cardiac electrophysiological simulation, which allows us to fully explain the differences between the solvers. We identify the use of mass lumping as the fundamental cause of the largest variations, specifically the combination of the commonly used techniques of mass lumping and operator splitting, which results in a slightly different form of mass lumping to that supported by theory and leads to increased numerical error. Other variations are explained through the manner in which the ionic current is interpolated. We also investigate the effect of different forms of mass lumping in various types of simulation.