[论文解读] Implementation of PDE models of cardiac dynamics on GPUs using OpenCL

Christopher D. Marcotte, Roman O. Grigoriev|arXiv (Cornell University)|Sep 6, 2013

Cardiovascular Function and Risk Factors参考文献 7被引用 6

一句话总结

本文提出了一种基于OpenCL的高性能双精度GPU实现，用于心脏动力学的偏微分方程（PDE）模型，实现了2D心肌组织的精确高效模拟。结果表明，四阶龙格-库塔方法在速度和精度上均优于低阶格式，使得在GPU上实现复杂心脏动力学的实时模拟成为可能，且精度损失极小。

ABSTRACT

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.

研究动机与目标

• 开发一种基于GPU的双精度数值积分器，用于心肌PDE模型，以克服现有单精度GPU代码的局限性。
• 评估并比较高阶显式有限差分方法（如四阶龙格-库塔法）在模拟2D心肌组织动力学时的精度与性能。
• 使用OpenCL实现并优化该方法，以实现跨异构GPU架构的可移植性，避免使用CUDA等厂商专用工具链。
• 解决GPU上双精度计算的数值挑战，包括超越函数的性能损耗以及内存访问模式问题。
• 实现对混沌动力学和无矩阵牛顿-Krylov求解器的精确模拟，这些方法对稳定性与收敛性要求双精度计算。

提出的方法

• 研究采用主机-内核混合架构，由Matlab主机负责控制，OpenCL内核在GPU上执行时间积分循环。
• 对描述心肌组织的反应-扩散PDE模型采用有限差分格式离散化，所有模型均采用诺伊曼（无通量）边界条件。
• 实现了三种心脏模型：Fitzhugh-Nagumo、Karma和Bueno-Orovio-Cherry-Fenton，每种模型具有不同的非线性反应项和扩散结构。
• 通过显式龙格-库塔格式实现高阶时间积分，特别是四阶方法，以在精度上优于低阶方法。
• 性能优化包括通过多项式逼近或替代公式，减少对慢速特殊函数（如tanh、pow）的使用。
• 通过将网格点映射到独立线程，充分利用GPU并行性，结合高效的内存合并与共享内存使用，实现紧致算子的计算。

实验结果

研究问题

• RQ1 使用OpenCL的双精度GPU内核能否实现足够性能，以支持心肌PDE模型的高保真模拟？
• RQ2 在GPU加速环境下，高阶时间积分格式（如四阶龙格-库塔法）与低阶方法相比，在性能与精度方面表现如何？
• RQ3 双精度GPU计算在心肌动力学模拟中的性能瓶颈是什么，特别是超越函数与内存访问方面？
• RQ4 有限差分紧致算子在GPU架构上的映射效率如何，适用于2D心肌组织模拟？
• RQ5 在OpenCL中能否有效实现基于谱变换（如傅里叶或切比雪夫）的算子分裂方法，以提升稳定性和效率？

主要发现

• 尽管四阶龙格-库塔法计算成本更高，但由于允许更大的时间步长，其整体运行速度显著快于低阶方法。
• 使用OpenCL的双精度计算在GPU上是可行的，但特殊函数单元（SFU）性能较慢，尤其影响非整数幂函数和双曲正切函数的计算。
• 用多项式逼近（如2Hk(x)−1）替代tanh等超越函数，可在不损失精度的前提下显著提升性能。
• 拉普拉斯算子的有限差分紧致算子在GPU并行计算中表现良好，支持2D心肌动力学的实时模拟。
• 采用谱方法（如傅里叶或切比雪夫）的半隐式算子分裂方法具有显著的未来性能提升潜力，但当前OpenCL支持限制了其可移植性。
• 三维模拟预计能从GPU加速中获益更大，使此前难以求解的问题变得计算可行。