Implementation of PDE models of cardiac dynamics on GPUs using OpenCL
School of Physics, Georgia Institute of Technology, Atlanta, GA 30332-0430, USA
arXiv:1309.1720 [nlin.CD], (6 Sep 2013)
@ARTICLE{2013arXiv1309.1720M,
author={Marcotte}, C.~D. and {Grigoriev}, R.~O.},
title={"{Implementation of PDE models of cardiac dynamics on GPUs using OpenCL}"},
journal={ArXiv e-prints},
archivePrefix={"arXiv"},
eprint={1309.1720},
primaryClass={"nlin.CD"},
keywords={Nonlinear Sciences – Chaotic Dynamics, Physics – Computational Physics},
year={2013},
month={sep},
adsurl={http://adsabs.harvard.edu/abs/2013arXiv1309.1720M},
adsnote={Provided by the SAO/NASA Astrophysics Data System}
}
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.
September 9, 2013 by hgpu