[論文レビュー] Algebraic criteria for A-stability of peer two-step methods

Bernhard A. Schmitt|arXiv (Cornell University)|Jun 18, 2015

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

ひとこと要約

本稿では、2つの正定値重み行列から導かれるテスト行列の半正定性に基づいて、ピア2段階法のA安定性の新しい代数的基準を提示する。この手法により、有理数演算とコレスキー分解を用いてA安定性を正確に検証可能となり、次数 $ s-1 $ のA安定ピア法の明示的パラメータ化が可能となる。また、$ s=3,4 $ に対して、対角的陰的および並列的なL安定法が構築される。

ABSTRACT

A new criterion for A-stability of peer two-step methods is presented which is verifiable exactly in exact arithmetic by checking semi-definiteness of a certain test matrix. It depends on the existence of two positive definite weight matrices for a given method. Although the initial approach is different using properties of the numerical radius the criterion itself resembles the one from algebraic stability of General Linear Methods. Known numerical algorithms for the computation of the unknown weight matrices suffer from rank deficiencies of the test matrix. For $s$-stage peer methods of order $s-1$ this rank defect is identified with an explicit block diagonal decomposition of the test matrix in trivial and definite blocks. In the design of methods its coefficients are unknown and an explicit parametrization of A-stable peer methods of order $s-1$ is presented with a weight matrix as parameter. This leads to a general existence result for any number of stages. The restrictions for efficient L-stable peer methods like diagonally-implicit and parallel ones are also discussed and such methods with $3$ and $4$ stages are constructed.

研究の動機と目的

• ピア2段階法のA安定性に対する、数値的固有値計算に依存しない検証可能な代数的基準の開発。
• $ z=0 $ においてスペクトル半径が正確に1となるため、標準的な浮動小数点固有値チェックに起因する不安定性を解消する。
• ガウス＝ルジャンドル法における代数的安定性の概念を、行列値安定関数を用いてピア2段階法へ一般化する。
• 実用的な剛性ODEソルバーにおける使用を想定し、効率的なL安定ピア法（対角的陰的および並列型）の構築。
• 任意の段数 $ s $ に対して存在が保証される、次数 $ s-1 $ のA安定ピア法のパラメータ化された族の提供。

提案手法

• 行列値安定関数 $ \mathbf{M}(z) $ に対するスペクトル半径条件 $ \varrho(\mathbf{M}(z)) < 1 $ を用いてA安定性を定式化。$ \Re(z) < 0 $ の範囲で成り立つ。
• サイズ $ (2s) \times (2s) $ のテスト行列 $ \mathbb{M} $ を導入。この行列の半正定性がA安定性を保証する。この行列は、未知の対称正定値重み行列 $ Z $ および $ W $ から導出される。
• 行列の数値半径の性質を用いて、スペクトル半径条件を $ \mathbb{M} $ の半正定性条件に再定式化。これはGLMsにおける代数的安定性に類似している。
• 半正定性計画法の実行可能性を用いて $ Z $ および $ W $ を特定。既知のアルゴリズムにおけるランク不足問題は、次数 $ s-1 $ の $ s $ 段法に対してブロック対角分解を用いて解消。
• 重み行列 $ Z $ を自由パラメータと見なすことで、A安定ピア法のパラメータ化された族を構築。これにより体系的な設計が可能となる。
• 浮動小数点誤差を避けるために、有理数演算とコレスキー分解を用いて結果を検証。特に $ \varrho(\mathbf{M}(0)) = 1 $ の場合に重要となる。

実験結果

リサーチクエスチョン

• RQ1 ピア2段階法のA安定性に対する、正確な算術計算で検証可能な代数的基準を定式化できるか？
• RQ2 A安定性検証において、$ z=0 $ における $ \mathbf{M}(z) $ の固有値計算に起因する数値的不安定性はどのように克服できるか？
• RQ3 正定値重み行列 $ Z $ および $ W $ は、ピア法におけるA安定性の特徴付けにおいて果たす役割は何か？また、これらを信頼性を持って計算する方法は？
• RQ4 任意の $ s $ に対して、次数 $ s-1 $ のA安定ピア法の一般化されたパラメータ化を導出できるか？また、その存在を保証できるか？
• RQ5 対角的陰的または並列スキームのような効率的なL安定ピア法を構築するための構造的制約は何か？また、それらを満たす方法は？

主な発見

• 本稿では、$ (2s) \times (2s) $ のテスト行列 $ \mathbb{M} $ の半正定性に基づく、A安定性の新しい代数的基準を確立した。この基準は、有理数演算とコレスキー分解を用いて正確に検証可能である。
• $ s $ 段法の次数 $ s-1 $ のピア法において、テスト行列はランク欠損を示すが、これは明示的な分解により自明なブロックと正定値ブロックに分解され、信頼性の高い計算が可能となる。
• 重み行列 $ Z $ を自由パラメータと見なすことで、次数 $ s-1 $ のA安定ピア法のパラメータ化された族が構築され、任意の $ s \geq 3 $ に対して存在が保証される。
• $ s=3 $ の場合、並列3段L安定ピア法が構築され、重み行列 $ Z $ および $ W $ の最小固有値はそれぞれ $ \approx 0.122 $ および $ \approx 0.021 $ であり、テスト行列 $ \hat{\mathbb{M}} $ はランク3で非自明な固有値が $[0.12, 40.16]$ の範囲に存在する。
• $ s=4 $ の場合、並列4段L安定ピア法が構築され、$ Z $ および $ W $ の最小固有値はそれぞれ $ \approx 0.109 $ および $ \approx 0.062 $ であり、テスト行列 $ \mathbb{M} $ はランク6で非自明な固有値が $[0.021, 194.1]$ の範囲に存在する。
• 数値実験により、変換（例：$ U_H $ を用いたもの）が行われない場合、小さな負の固有値（例：$ -6 \cdot 10^{-9} $）が生じるため、半正定性は信頼性を持って検証できないことが確認された。これにより、提案された行列変換の必要性が強調された。