[論文レビュー] Efficient adaptivity for simulating cardiac electrophysiology with spectral deferred correction methods

Fatemeh Chegini, Thomas Steinke+1 more|arXiv (Cornell University)|Nov 13, 2023

Cardiac electrophysiology and arrhythmias被引用数 4

ひとこと要約

本稿では、心電図シミュレーションにおけるスペクトル遅延補正（SDC）法に代数的自己適応を提案し、後続のSDCスイープにおいて自由度（dofs）の部分集合を段階的に選択することで計算コストを低減する。解の特徴が局所的であることに着目し、事後誤差推定を活用することで、メッシュベースの適応と比較してわずかなオーバーヘッドで、モノドメインおよびEMIモデルにおいて2–4倍の高速化を達成した。

ABSTRACT

The locality of solution features in cardiac electrophysiology simulations calls for adaptive methods. Due to the overhead incurred by established mesh refinement and coarsening, however, such approaches failed in accelerating the computations. Here we investigate a different route to spatial adaptivity that is based on nested subset selection for algebraic degrees of freedom in spectral deferred correction methods. This combination of algebraic adaptivity and iterative solvers for higher order collocation time stepping realizes a multirate integration with minimal overhead. This leads to moderate but significant speedups in both monodomain and cell-by-cell models of cardiac excitation, as demonstrated at four numerical examples.

研究の動機と目的

• 高分解能の心電図シミュレーションにおける高い計算コストに対処すること、特にモノドメインおよびEMIモデルにおいて。
• メッシュの細分化・粗化および行列構築に伴う高いオーバーヘッドにより、従来のメッシュベースの空間的適応が非効率であるという問題を克服すること。
• メッシュの変更なしに線形方程式系のレベルで動作する、低オーバーヘッドの代数的適応戦略を開発すること。
• SDC反復中に活性化フロントが存在する領域に自由度を動的に制限することで、効率的なマルチレート統合を可能にすること。
• 代数的適応が、精度と許容誤差要件を維持したままシミュレーション時間を顕著に短縮できることを実証すること。

提案手法

• モノドメインおよびEMIモデルを解くために、高次の精度を達成するための、implicit-explicit（IMEX）時間積分を用いたスペクトル遅延補正（SDC）法を適用する。
• SDCの線形収束に基づく事後誤差推定を用い、後続のスイープで最も重要な自由度（dofs）のみを特定・保持する。
• 段階的部分集合選択を実装する：最初のSDCスイープではすべてのdofsを使用するが、以降のスイープでは補正寄与が顕著なdofsに制限することで、システムサイズを低減する。
• SDC反復と適応的dofs選択をインターリーブすることで、解の活動が低い領域での計算負荷を最小限に抑えつつ、精度を維持する。
• 補正への影響が無視できるdofsを除外するためのドロップ許容誤差基準を用い、手法の効率性と安定性を保証する。
• 2次元および3次元のモノドメインおよびEMIモデルにこの手法を適用し、壁時計時間と精度の観点から、適応的および非適応的SDCベースラインを比較する。

実験結果

リサーチクエスチョン

• RQ1 SDC補正の局所化に基づく代数的適応は、顕著なオーバーヘッドを伴わず、心電図シミュレーションにおける計算コストを低減できるか？
• RQ2 後続のSDCスイープにおけるdofsの部分集合選択は、システムサイズおよびシミュレーション時間の低減にどの程度効果的か？
• RQ3 非適応的SDCと比較して、モノドメインおよびEMIモデルにおけるこの手法で達成できるスピードアップはどの程度か？
• RQ4 特に3次元シミュレーションにおいて、問題のスケールや活性化フロントのダイナミクスに応じて、この手法の性能はどのように変化するか？
• RQ5 この手法の効率性は、SDCスイープの数やドロップ許容誤差の選択にどの程度依存するか？

主な発見

• 提案された代数的適応により、2次元モノドメインモデルで2.13倍、3次元モノドメインモデルで3.29倍のスピードアップが達成された。EMIモデルではさらに高い増加率を示し、それぞれ3.336倍および4.344倍のスピードアップを記録した。
• 3次元モノドメインの例では、最初のスイープで646,166個の活性dofsが、3番目のスイープで425個にまで減少し、補正が活性化フロントに強く局在していることが示された。
• 許容誤差（TOL = 1e-4）の範囲内で精度を維持しながら、後続のSDCスイープでdofsを最大99.9％まで削減できた。
• 誤差推定およびdofs選択に起因するオーバーヘッドは無視できるほど小さく、システムサイズの低減が直接的にシミュレーション時間の短縮に反映された。
• 領域が大きく、活性化フロントの割合が小さい場合、補正の空間的サポートがより局所的になるため、スピードアップが増大する傾向を示した。
• 本手法は、SDCスイープ数が中程度（3–10）の際に最も効果的であり、今後の改善には高次のコロケーションやヒューリスティックな最初のスイープにおけるdofs選択が含まれる可能性がある。