[論文レビュー] Inexact spectral deferred corrections using single-cycle multigrid

Robert Speck, Daniel Ruprecht+3 more|arXiv (Cornell University)|Jan 30, 2014

Advanced Numerical Methods in Computational Mathematics参考文献 11被引用数 3

ひとこと要約

本稿では、剛性初期値問題に対して、大規模な連立線形方程式を効率的に解くために、1サイクルマルチグリッドを活用する非正確なスペクトル遅延補正（SDC）法を提案する。SDCの反復補正フレームワークとマルチグリッド予め条件付けを組み合わせることで、計算コストを低減しつつ高次の精度を達成し、標準的なSDCに比べて収束速度を顕著に向上させる。

ABSTRACT

Implicit integration methods based on collocation are attractive for a num-ber of reasons, e.g. their ideal (for Gauss-Legendre nodes) or near ideal(Gauss-Radau or Gauss-Lobatto nodes) order and stability properties. How-ever, straightforward application of a collocation formula with M nodes toan initial value problem with dimension d requires the solution of one largeMd×Md system of nonlinear equations.Spectral deferred correction (SDC) methods, introduced by Dutt et al.[2000],arean attractiveapproachforiterativelycomputing the solution to thecollocation problem using a low-order method (like implicit or IMEX Euler)as a building block. Instead of solving one huge system of size Md × Md,SDC iteratively solves M smaller d×d systems to approximate the solutionof the full system (see also the discussion in Huang et al. [2006]). It has beenshown e.g. by Xia et al. [2007] that each iteration/sweep of SDC raises theorder by one, so that SDC with k iterations and a first-order base method isof order k, up to the order of the underlying collocation formula. Therefore,to achieve formal order p, SDC requires p/2 nodes and p iterations and thusp

研究の動機と目的

• 剛性初期値問題に対する陰的コロケーション法における大規模かつ結合された方程式系を解く際の高い計算コストを低減すること。
• 収束に必要な反復回数を削減することで、スペクトル遅延補正（SDC）法の効率を向上させること。
• マルチグリッド技術をSDCに統合し、各SDCスイープにおける線形方程式の予め条件付けを効果的に行い、収束を加速すること。
• 多次元問題における大規模な方程式系の解法に伴う計算負担を軽減しながらも、高次の精度を維持すること。
• 非正確SDCに1サイクルマルチグリッドを適用した場合、最小限の追加オーバーヘッドで最適な収束率を達成できることを示すこと。

提案手法

• 各M個のコロケーションノードで、低次の基本法（例：陰的エイラー法）を用いて反復的に解を補正する、スペクトル遅延補正（SDC）を基本フレームワークとして用いる。
• 各SDCスイープで生じる線形方程式を解くために、直接法や反復法の代わりに、マルチグリッド予め条件付けアプローチを採用する1サイクルマルチグリッド法を適用する。
• 収束を加速するために、粗いグリッドの階層を用い、各SDCスイープあたり必要な反復回数を削減する。
• マルチグリッド解に対して制御された誤差を許容する非正確な解法を導入し、精度と計算コストのバランスを取る。
• 補正ステップが各スイープにわたり一貫性を保つようにすることで、元のコロケーション式の精度の順序を維持する。
• SDCの反復的性質とマルチグリッドの高速収束性を組み合わせることで、計算作業を低減しつつ高次の時間積分を達成する。

実験結果

リサーチクエスチョン

• RQ1 1サイクルマルチグリッドは、SDCにおける線形方程式を効果的に予め条件付けし、精度を損なわせることなく収束を加速できるか？
• RQ2 マルチグリッドを用いた非正確SDC法は、標準的なSDCに比べて収束速度と計算コストの点でどのように異なるか？
• RQ3 マルチグリッドサイクル内で非正確な解法を許容することで、SDC法全体の次数と安定性にどのような影響を与えるか？
• RQ4 SDCとマルチグリッドの組み合わせは、反復回数を顕著に削減しつつも、高次の精度を維持できるか？
• RQ5 提案手法は、問題次元とコロケーションノード数の増加に対して、効率的にスケーリングできるか？

主な発見

• 非正確SDC法に1サイクルマルチグリッドを適用した場合、正確なSDCと同等の収束速度を達成するが、反復ごとの計算コストは著しく低減された。
• マルチグリッド予め条件付けの導入により、特に条件数が高い剛性問題において、各SDCスイープに必要な反復回数が削減された。
• 線形方程式を非正確に解いても、元のコロケーション式の形式的精度の順序を維持したままである。
• 1サイクルマルチグリッド手法は、精度と効率の間で良好なトレードオフを提供し、正確なSDCに比べて収束を高速化しつつも、次数の損失がない。
• SDC反復ごとの計算作業は、マルチグリッド収束速度に比例して低減され、大規模問題において顕著な節約効果が得られた。
• 数値実験により、非正確SDCにマルチグリッドを適用した手法は、さまざまな剛性ODEおよび問題サイズの拡大に対しても、安定性と精度を維持していることが確認された。