[論文レビュー] 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上で2次元心筋組織の高精度かつ高効率なシミュレーションを可能にする、高精度な二重精度GPU実装を提示している。偏微分方程式（PDE）モデルを用いた心臓動態のシミュレーションにおいて、4次ルンゲ・クッタ法が低次のスキームよりも速度と精度に優れ、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コードの制限を克服し、GPU上で二重精度の数値積分器を開発すること。
• 2次元心筋組織動態のシミュレーションにおける精度とパフォーマンスを評価・比較する、高次明示的有限差分法（例：4次ルンゲ・クッタ）の検討。
• OpenCLを用いた実装と最適化により、ベンダー固有のツールチェイン（例：CUDA）に依存しない、異種GPUアーキテクチャ間での移植性を確保すること。
• トランセンドental関数の性能ペナルティやメモリアクセスパターンの影響を含む、GPU上での二重精度計算における数値的課題の解決。
• 不安定性や収束性を確保するため二重精度を必要とする、混沌とした動的挙動や行列を用いないニュートン・クリロフ法の正確なシミュレーションを可能にすること。

提案手法

• 制御はMatlabホストで行い、時間積分ループのGPU実行はOpenCLカーネルを用いるハイブリッドホスト・カーネルアーキテクチャを採用している。
• 心筋組織をモデル化する反応拡散PDEに対して、有限差分法による離散化を適用し、すべてのモデルでノイマン（不透過）境界条件を設定している。
• 3つの心臓モデルを実装した：フィッツフーリュ・ナグumo、カルマ、ブルーノ・オロビオ・チェリー・フェントン。各モデルは非線形反応項と拡散構造が異なる。
• 時間積分の高次化は、特に4次ルンゲ・クッタ法を用いた明示的ルンゲ・クッタスキームにより実現し、低次の手法に比べて精度を向上させている。
• 性能最適化として、tanh や pow などの遅延する特殊関数の使用を最小限に抑え、多項式近似や代替式に置き換える手法を採用している。
• GPUの並列性を活用し、グリッド点を個々のスレッドにマッピングし、ステンシル計算に適したメモリコalescingと共有メモリの使用を実装している。

実験結果

リサーチクエスチョン

• RQ1 OpenCLを用いた二重精度GPUカーネルは、心臓PDEモデルの高精度なシミュレーションに十分なパフォーマンスを発揮できるか？
• RQ2 GPUアクセcelレーション環境下で、高次時間積分スキーム（例：4次ルンゲ・クッタ）は低次の手法に比べて、性能と精度の点でどのように差がつくか？
• RQ3 特にトランセンドental関数とメモリアクセスの点で、心臓動態のための二重精度GPU計算におけるパフォーマンスボトルネックは何か？
• RQ4 2次元心筋組織シミュレーションにおいて、有限差分ステンシルはGPUアーキテクチャに効率的にマッピング可能か？
• RQ5 スペクトル変換（例：フーリエまたはチェビシェフ）を用いたオペレータ分割法は、OpenCLで効果的に実装可能か？安定性と効率性の向上が見込まれる。

主な発見

• 4次ルンゲ・クッタ法は、計算コストが高めでも、許容可能な時間ステップが大きくとれるため、低次の手法よりも著しく高速である。
• OpenCLを用いたGPU上での二重精度計算は実現可能であるが、特に非整数乗数関数や双曲線正接関数の特殊関数ユニット（SFU）の遅延によりパフォーマンスが低下する。
• tanh などのトランセンドental関数を多項式近似（例：2Hk(x)−1）に置き換えることで、精度を損なわずパフォーマンスが向上する。
• ラプラシアン演算子の有限差分ステンシルはGPUの並列性に適しており、2次元心臓動態のリアルタイムシミュレーションを可能にしている。
• スペクトル法（例：フーリエまたはチェビシェフ）を用いた半陰的オペレータ分割法は、将来的なパフォーマンス向上の強力な可能性を秘めているが、現在のOpenCLのサポート制限により移植性が制限される。
• 3次元シミュレーションは、GPUアクセcelレーションによりさらに大きな恩恵を受けると予想され、かつては計算不能とされた問題が現在は計算可能になる。