摘要：Cyclic calcium release and uptake in heart cells has an important role in heart rhythm and contraction, and it is known that the malregulation of calcium release is a predictor for cardiac arrhythmia. A model of this calcium release process was proposed, which consists of a large number of discrete calcium release sites, each involving stiff, stochastic, and nonlinear systems. In this paper we have developed a simulation of this calcium release model, which is parallel across the problem. The simulation is developed in a CUDA framework to be solved using GPUs. Computational efficiency is enhanced by using a DIIRK solver and taking advantage of the sparsity of the Jacobian. The output is shown to display behaviour similar to empirical observations, in particular displaying behaviour known as calcium waves.
