Implementation of PDE models of cardiac dynamics on GPUs using OpenCL

Graphical processing units (GPUs) promise to revolutionize scientific computing in the near future. Already, they allow almost real-time integration of simplified numerical models of cardiac tissue dynamics. However, the integration methods that have been developed so far are typically of low order and use single precision arithmetics. In this work, we describe numerical implementation of double precision integrators required by, e.g., matrix-free Newton-Krylov solvers and compare several higher order, fully explicit numerical methods using finite-difference discretization of a range of models of two-dimensional cardiac tissue.