[論文レビュー] Optimal Runge--Kutta Stability Regions

David I. Ketcheson, Aron Ahmadia|arXiv (Cornell University)|Jan 14, 2012

Numerical methods for differential equations参考文献 21被引用数 3

ひとこと要約

本稿では、スペクトルが既知の初期値ODEおよびPDEに対して、最小偏差実行可能性問題を解くことで、最適に安定したルンゲ＝クッタ法を計算する凸最適化に基づくアルゴリズムを提案する。この手法は、初期値問題の安定ステップサイズを最大化し、特定の条件下でグローバル収束を保証するとともに、それらの条件を満たさない場合でも頑健な性能を示す。

ABSTRACT

We consider the problem of finding optimally stable polynomial approximations to the exponential for application to one-step integration of initial value ordinary and partial differential equations. The objective is to find the largest stable step size and corresponding method for a given problem when the spectrum of the initial value problem is known. The problem is expressed in terms of a general least deviation feasibility problem. Its solution is obtained by a new fast, accurate, and robust algorithm based on convex optimization techniques. Global convergence of the algorithm is proven in the case that the order of approximation is one and in the case that the spectrum encloses a starlike region. Examples demonstrate the effectiveness of the proposed algorithm even when these conditions are not satisfied.

研究の動機と目的

• スペクトルが既知の初期値問題に対する1ステップ積分において、最大の安定ステップサイズを特定すること。
• 指数関数への多項式近似のための一般化された最小偏差実行可能性問題として問題を定式化すること。
• 凸最適化技術に基づく、高速かつ高精度で頑健なアルゴリズムの開発。
• 近似次数が1またはスペクトルがスターリーク型領域を形成する場合に、アルゴリズムのグローバル収束を保証すること。
• 理論的収束条件を越えて実用的適用性を拡張するため、数値的検証を通じてその有効性を示すこと。

提案手法

• 多項式近似が指数関数からどれほど離れているかを最小化する最小偏差実行可能性問題として問題を定式化する。
• 効率的かつ高精度に実行可能性問題を解くために、凸最適化フレームワークを採用する。
• 与えられたスペクトルと近似次数に対して、安定ステップサイズを最大化するようにアルゴリズムを設計する。
• 近似次数が1またはスペクトルがスターリーク型領域を形成する特殊ケースにおいて、グローバル収束を証明する。
• スペクトルから導かれる多項式安定領域を最適化プロセスの指針として用いる。
• 数値実験により、収束条件を満たさない場合でも頑健な性能を確認する。

実験結果

リサーチクエスチョン

• RQ1 スペクトルが既知の初期値問題に対して、ルンゲ＝クッタ法を用いて到達可能な最大の安定ステップサイズは何か？
• RQ2 1ステップ積分における最大安定性を実現するため、指数関数への多項式近似をどのように最適化できるか？
• RQ3 どのような状況で最適化アルゴリズムのグローバル収束を保証できるか？
• RQ4 スペクトル領域がスターリーク型でない、または近似次数が1を超過する場合、提案手法はどのように動作するか？
• RQ5 一般のスペクトル構成下でも、アルゴリズムは精度と頑健性を維持できるか？

主な発見

• 提案アルゴリズムは、近似次数が1の場合にグローバル収束を達成する。
• スペクトルがスターリーク型領域を囲む場合にも、グローバル収束が確立される。
• 理論的収束条件を満たさない場合でさえ、最適に安定したステップサイズを効果的に計算可能である。
• 数値例により、実用的状況下でのアルゴリズムの頑健性と有効性が確認された。
• 凸最適化フレームワークにより、ルンゲ＝クッタ法の安定領域を高速かつ高精度に計算可能である。
• スペクトル情報に基づいて最大安定ステップサイズを体系的に特定する手法を提供する。