[論文レビュー] A time-parallel algorithm for almost integrable Hamiltonian systems

Hugo Jiménez-Pérez, J. Laskar|arXiv (Cornell University)|Jun 18, 2011

Numerical methods for differential equations参考文献 1被引用数 8

ひとこと要約

本稿では、摂動項の微分情報を用いたピカード反復と作用角座標を活用することで、時間並列化されたほぼ可積分なハミルトニアン系の数値積分を高速化するアルゴリズムを提示する。この手法は、逐次積分器への完全収束を達成し、速度向上の上限が 1/a に制限される。SST97アルゴリズムを上回る収束性を実現するため、高次微分の使用により収束を改善している。

ABSTRACT

We introduce a time-parallel algorithm for solving numerically almost integrable Hamiltonian systems in action-angle coordinates. This algorithm is a refinement of that introduced by Saha, Stadel and Tremaine in 1997 (SST97) for the same type of problems. Our refined algorithm has a better convergence obtained from the use of derivatives of the perturbing term not considered in the original SST97 algorithm. An advantage of this algorithm is its independence of the step-size for the parallelized procedures which can be consider as a particular case of the parareal scheme.

研究の動機と目的

• ほぼ可積分なハミルトニアン系のための時間並列数値アルゴリズムを構築すること。
• SST97法の収束が遅い問題を、摂動項の微分を組み込むことで改善すること。
• 時間刻みサイズに依存しない並列計算を可能とし、大規模かつ長期的なシミュレーションに適すること。
• 逐次積分器への収束を達成しつつ、予測子ステップの計算コストを最小限に抑えること。

提案手法

• アルゴリズムは、時間区間全体にわたり作用角座標上でピカード反復を用いて解を繰り返し改善する。
• 時間領域を重複するウィンドウに分割し、初期条件の並列補正を可能にする。
• 補正ステップで摂動項の微分を用いることで、元のSST97法を上回る収束速度を達成する。
• 粗いソルバと精度の高い補正ステップを組み合わせた、パラレル・タイム・スケジューリングに類似したスキームに基づく。補正ステップでは複数発射法とニュートン反復を用いる。
• 最適なウィンドウサイズ P* は解析的に導出され、P* = √(bT_P / aT_C) として与えられ、与えられた問題における速度向上を最大化する。
• アルゴリズムはSBAB4積分法を用いて実装され、摂動評価コストが高価な問題（T_B ≪ T_A）に対して最適化されている。

実験結果

リサーチクエスチョン

• RQ1 ほぼ可積分なハミルトニアン系において、SST97時間並列アルゴリズムの収束速度をどのように向上させられるか？
• RQ2 摂動項の微分が時間並列ソルバの収束をどのように加速するか？
• RQ3 逐次解への収束を保証しつつ速度向上を最大化する最適なウィンドウサイズ P* は何か？
• RQ4 計算コストと問題次元数の変化に伴うアルゴリズムのスケーリング特性は？特に太陽系のような高次元系においては？
• RQ5 時間刻みサイズに依存しない上限のある速度向上を達成できるか？その上限は何かによって制限されるか？

主な発見

• 提案手法は、逐次積分器への100%の収束を達成し、正確性を保証する。
• 速度向上は 1/a に上限がある。ここで a は粗いソルバのコストに関連するパrameterであり、理論的上限を示している。
• 最適なウィンドウサイズ P* は P* = √(bT_P / aT_C) として解析的に導出され、与えられた問題における速度向上を最大化する。
• 摂動評価コスト T_B が可積分部分のコスト T_A よりも著しく小さい（T_B ≪ T_A）場合に最も効果的である。
• 太陽系のような高次元系では、連結された積分器の数 j が制限される（例：j ≤ 8）。これ以上になると性能が劣化する。
• 予測子ステップ（T_P）と補正ステップ（T_C）の実行時間をバランスさせることで、並列スレッドの無駄な待機時間を削減する。最適効率を得るには T_P ≈ (P - C_Δt)T_C である必要がある。