Efficient Numerical Schemes for Computing Cardiac Electrical Activation over Realistic Purkinje Networks: Method and Verification

We present a numerical solver for the fast conduction system in the heart using both a CPU and a hybrid CPU/GPU implementation. To verify both implementations, we construct analytical solutions and show that the L-2-error is similar in both implementations and decreases linearly with the spatial step size. Finally, we test the performance of the implementations with networks of varying complexity, where the hybrid implementation is, on average, 5.8 times faster.