[論文レビュー] A theoretical comparison of high order explicit Runge-Kutta, extrapolation, and deferred correction methods

David I. Ketcheson, Umair bin Waheed|arXiv (Cornell University)|May 27, 2013

Numerical methods for differential equations参考文献 20被引用数 2

ひとこと要約

この論文は、初期値問題に対する高次1段階解法—補外法、スペクトル遅延補正法、埋め込みルンゲ・クッタ対—を比較し、それらを固定次数ルンゲ・クッタ法として定式化する。並列化された補外法は、タイトな許容誤差においてDOP8（8次）を上回り、N=400のN体問題では最大2倍速くなることが判明した。

ABSTRACT

We compare the three main types of high-order one-step initial value solvers: extrapolation, spectral deferred correction, and embedded Runge--Kutta pairs. We consider orders four through twelve, including both serial and parallel implementations. We cast extrapolation and deferred correction methods as fixed-order Runge--Kutta methods, providing a natural framework for the comparison. The stability and accuracy properties of the methods are analyzed by theoretical measures, and these are compared with the results of numerical tests. In serial, the 8th-order pair of Prince and Dormand (DOP8) is most efficient. But other high order methods can be more efficient than DOP8 when implemented in parallel. This is demonstrated by comparing a parallelized version of the well-known ODEX code with the (serial) DOP853 code. For an $N$-body problem with $N=400$, the experimental extrapolation code is as fast as the tuned Runge--Kutta pair at loose tolerances, and is up to two times as fast at tight tolerances.

研究の動機と目的

• 4次から12次までの範囲で、高次1段階解法—補外法、スペクトル遅延補正法、埋め込みルンゲ・クッタ対—の効率性と精度を比較すること。
• 理論的指標を用いて安定性と精度を分析し、数値実験でその妥当性を検証すること。
• これらの手法のシリアルおよび並列実装の実用的性能を評価すること。
• 特に並列化された手法が、DOP853のような標準的な解法を上回る条件を特定すること。

提案手法

• 著者たちは、補外法と遅延補正法を固定次数ルンゲ・クッタスキームに再定式化することで、統一的な比較フレームワークを構築した。
• ルンゲ・クッタ法に適用可能な理論的指標を用いて、安定性と精度を分析した。
• 許容誤差レベルを変化させたN=400のN体問題を用いて、性能を比較する数値実験を実施した。
• O-DEX補外コードの並列化版を、シリアルDOP853コードと比較してベンチマークし、スケーラビリティと効率性を評価した。
• 4次から12次までの次数において、シリアルおよび並列実装の両方を含めた比較を実施した。
• 理論的安定性および精度特性と、観察された数値的性能の相関を検証した。

実験結果

リサーチクエスチョン

• RQ1 シリアルおよび並列環境下で、初期値問題に対するどの高次1段階法が最も効率的か？
• RQ2 4次から12次までの範囲で、補外法、遅延補正法、埋め込みルンゲ・クッタ対の安定性および精度特性はどのように比較されるか？
• RQ3 並列化された補外法は、高度に最適化されたシリアルDOP853コードを上回ることができるか？
• RQ4 タイトな許容誤差と緩い許容誤差の間で、高次法の性能優位性はどのように変化するか？
• RQ5 手法の選択が、大規模なN体シミュレーションにおける計算コストにどのように影響するか？

主な発見

• シリアルモードでは、8次精度のプリンスとドゥルマンド対（DOP8）が、テストされた手法の中で最も効率的であった。
• 並列化された補外法は、シリアルDOP853コードを上回り、タイトな許容誤差において最大2倍の高速化が達成された。
• N=400のN体問題において、実験的補外コードは緩い許容誤差ではDOP853と同等の速度を示し、タイトな許容誤差ではそれを上回った。
• 他の高次法も、並列化された実装ではDOP8を上回る効率性を示し、アーキテクチャに応じて最適な選択が変化することを示した。
• 理論的安定性および精度指標と、観察された数値的性能の間に良好な整合性が確認され、比較フレームワークの妥当性が裏付けられた。
• 結果から、特に補外法を用いた並列高次解法が、高精度シミュレーションにおいて有望であることが示唆された。