[论文解读] 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

一句话总结

本文比較了高階單步求解器——外推法、譜延遲校正法與嵌入式龍格-庫塔對—在初值問題中的表現，並將其視為固定階次的龍格-庫塔方法。研究發現，雖然在串行環境下DOP8（八階）最有效率，但在嚴格容差條件下，平行化的外推法反而超越DOP8，使N=400的N體問題運算速度最快可達兩倍。

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.

研究动机与目标

• 比較四至十二階之間高階單步求解器（外推法、譜延遲校正法與嵌入式龍格-庫塔對）在效率與準確性上的表現。
• 利用理論指標分析穩定性與準確性，並透過數值測試加以驗證。
• 評估這些方法在串行與平行環境下的實務效能。
• 確定在何種條件下，特別是平行化後，高階求解器能超越既有的DOP853等求解器。

提出的方法

• 作者將外推法與延遲校正法重新表述為固定階次的龍格-庫塔格式，以建立統一的比較框架。
• 利用適用於龍格-庫塔方法的理論指標分析穩定性與準確性。
• 針對N=400的N體問題進行數值實驗，比較不同容差水平下的效能表現。
• 將平行化的ODEX外推法程式碼與串行的DOP853程式碼進行基準測試，以評估可擴展性與效率。
• 比較四至十二階之間的串行與平行實作版本。
• 將理論上的穩定性與準確性特性與實際觀察到的數值表現進行關聯。

实验结果

研究问题

• RQ1 在串行與平行環境下，哪種高階單步求解法對初值問題最有效率？
• RQ2 外推法、延遲校正法與嵌入式龍格-庫塔對在四至十二階之間的穩定性與準確性特性如何比較？
• RQ3 平行化的外推法是否能超越高度優化的串行DOP853程式碼？
• RQ4 高階求解器在嚴格容差下相較於寬鬆容差下的效能優勢為何？
• RQ5 方法的選擇如何影響大規模N體模擬中的計算成本？

主要发现

• 在串行模式下，八階的Prince與Dormand對（DOP8）是在測試方法中效率最高的。
• 平行化的外推法可超越串行的DOP853程式碼，在嚴格容差下達成最多兩倍的加速效果。
• 對於N=400的N體問題，實驗性外推法程式碼在寬鬆容差下與DOP853速度相當，而在嚴格容差下則超越DOP853。
• 其他高階方法在平行化實作後，其效率亦超越DOP8，顯示最佳選擇會因架構而異。
• 理論上的穩定性與準確性指標與實際觀察到的數值表現高度一致，驗證了比較框架的有效性。
• 結果顯示，特別是外推法，平行高階求解器在高精度模擬中具有極大潛力。