[論文レビュー] Peer Methods for the Solution of Large-Scale Differential Matrix Equations

Peter Benner, Norman Lang|arXiv (Cornell University)|Jul 23, 2018

Matrix Theory and Algorithms参考文献 30被引用数 4

ひとこと要約

本稿では、大規模な微分行列方程式、特に微分リッカティ方程式（DRE）を解くための暗黙的およびローゼンブロック型ペア法を提案する。低ランクの対称不定因子分解と再定式化されたジャコビアンフリー手法を活用することで、計算コストを低減しつつ高い精度を維持し、大規模問題において従来のBDFおよびローゼンブロックスキームよりも効率的で安定した性能を発揮する。

ABSTRACT

We consider the application of implicit and linearly implicit (Rosenbrock-type) peer methods to matrix-valued ordinary differential equations. In particular the differential Riccati equation (DRE) is investigated. For the Rosenbrock-type schemes, a reformulation capable of avoiding a number of Jacobian applications is developed that, in the autonomous case, reduces the computational complexity of the algorithms. Dealing with large-scale problems, an efficient implementation based on low-rank symmetric indefinite factorizations is presented. The performance of both peer approaches up to order 4 is compared to existing implicit time integration schemes for matrix-valued differential equations.

研究の動機と目的

• 最適制御およびモデル低減に現れる大規模な微分リッカティ方程式（DRE）を効率的に解く課題に対処する。
• 行列値ODEの古典的暗黙時間積分スキームが直面する高い計算コストと記憶容量要件を克服する。
• DRE解の数値的ランク構造を活用する低ランク定式化を開発し、計算複雑度を低減する。
• ジャコビアン評価回数を最小限に抑えるために、ローゼンブロック型ペア法の再定式化を導入し、計算効率を向上させる。
• 大規模問題において、特にローゼンブロック型変種を含むペア法が、古典的BDFおよび暗黙スキームに比べて優れた性能を示すことを実証する。

提案手法

• 同一の精度および安定性特性を持つペア段階を有する一般線形法の一種であるペア法を行列値ODEに適用する。
• 微分リッカティ方程式（DRE）$ \dot{X}(t) = A(t)^T X(t) + X(t)A(t) - X(t)S(t)X(t) + W(t) $ の形に適応した、暗黙的および線形的暗黙的（ローゼンブロック型）ペアスキームを構築する。
• 段階評価を再利用することで、重複するジャコビアン適用を回避するローゼンブロック型ペア法の再定式化を開発し、計算オーバーヘッドを低減する。
• 各時間ステップで生じる代数的リッカティ方程式（AREs）を効率的に解くために、低ランク対称不定因子分解（LRSIF）を実装する。DRE解の低ランク構造を活用する。
• 解を $ X = ZZ^T $（$ Z \in \mathbb{R}^{n \times k} $、$ k \ll n $）という低ランクコレスキー型因子分解で表現し、記憶および計算コストを削減する。
• 大規模な設定でも完全な行列構築を避けるために、ペアスキームとKrylov部分空間法を統合する。

実験結果

リサーチクエスチョン

• RQ1 ペア法は、特にDREを対象として、効率的かつ安定的に大規模な微分行列方程式を解くために効果的に適応可能か？
• RQ2 提案されたローゼンブロック型ペア法の再定式化は、定常状態においてジャコビアン評価回数と計算コストをどの程度低減するか？
• RQ3 低ランク対称不定因子分解（LRSIF）は、大規模DREの時間積分におけるスケーラビリティをどの程度向上させるか？
• RQ4 大規模DREに対して、ペア法は古典的BDFおよびローゼンブロックスキームに比べて、精度および計算効率で優れているか？
• RQ5 DRE解の低ランク構造は、ペア統合スキームと組み合わせることで、顕著な性能向上をもたらすか？

主な発見

• 特に再定式化されたバージョン（mRosPeer）を含むローゼンブロック型ペア法は、大規模DREに対して、古典的BDFおよびローゼンブロックスキームを上回る計算効率を示した。
• 800ステップの対流拡散LTV問題において、mRosPeer(1)スキームは11.16秒で相対誤差2.09e-02を達成し、類似の精度を保ちながらBDF(1)を速度面で上回った。
• n = 2025の大きな規模のLTI対流拡散問題において、mRosPeer(2)スキームは543.07秒で実行されたが、BDF(3)（870.76 s）およびBDF(4)（813.07 s）よりも著しく速く、高い精度を維持した。
• Peer(4)およびRosPeer(4)は、段階数および次数の増加に伴い計算時間が長くなった（それぞれ2652.84 sおよび913.28 s）、これは次数と効率のトレードオフを示している。
• 再定式化されたローゼンブロック型ペアスキーム（mRosPeer）は、高次スキームにおいても、元のバージョン（RosPeer）を一貫して上回る計算時間性能を示した。
• 低ランク対称不定因子分解（LRSIF）により、各時間ステップにおけるAREの効率的解法が可能になり、大規模DREの解法が、完全な行列の保存や演算を伴わずに行えるようになった。