[論文レビュー] Runge-Kutta Theory and Constraint Programming

Julien Alexandre dit Sandretto|arXiv (Cornell University)|Apr 13, 2018

Numerical Methods and Algorithms参考文献 16被引用数 3

ひとこと要約

本稿では、区間解析を用いた制約プログラミングの新規アプローチを提案し、安定性、シンプレクティック性、および次数の正確性を保証するRunge-Kutta係数を計算する。Butcherテーブルの条件を制約充足問題として定式化し、グローバルコスト関数を用いた最適化により、従来の浮動小数点法よりも、解の包含範囲がよりきつい検証済みシミュレーションで優れた性能を示す区間係数Runge-Kuttaスキームを生成する。

ABSTRACT

There exist many Runge-Kutta methods (explicit or implicit), more or less adapted to specific problems. Some of them have interesting properties, such as stability for stiff problems or symplectic capability for problems with energy conservation. Defining a new method suitable to a given problem has become a challenge. The size, the complexity and the order do not stop growing. This informal challenge to implement the best method is interesting but an important unsolved problem persists. Indeed, the coefficients of Runge-Kutta methods are harder and harder to compute, and the result is often expressed in floating-point numbers, which may lead to erroneous integration schemes. Here, we propose to use interval analysis tools to compute Runge-Kutta coefficients. In particular, we use a solver based on guaranteed constraint programming. Moreover, with a global optimization process and a well chosen cost function, we propose a way to define some novel optimal Runge-Kutta methods.

研究の動機と目的

• 浮動小数点演算を用いたRunge-Kutta係数の正確な計算が難しくなる問題に対処し、数値積分における誤差の発生を回避する。
• 浮動小数点係数を区間係数に置き換えることで、次数、安定性、シンプレクティック性といった重要な性質を維持する。
• 制約プログラミングとグローバル最適化を用いて、新たなRunge-Kutta法を体系的かつ検証可能な方法で発見する手法を開発する。
• 区間係数スキームが保証された解の包含範囲を持つ検証済みシミュレーション環境での有効性を示す。
• 理論的性質を厳密な計算のもとで保持する区間係数を有する最適なRunge-Kutta法を生成するフレームワークを提供する。

提案手法

• Runge-Kutta法の条件（次数条件、安定性、シンプレクティック性）を区間変数における非線形方程式および不等式の系として定式化する。
• 保証付きの制約プログラミング（CP）ソルバを用いて、すべてのButcherテーブル条件を満たす区間ボックスを特定し、すべての解が有効であり、包含的に正しいことを保証する。
• 局所切り捨て誤差（LTE）を次に高い次数の方法に対して最小化するように設計されたコスト関数を用いたグローバル最適化手順を適用し、最適なスキームを優遇する。
• 構築および検証プロセスの全過程で区間演算を用いることで、厳密な誤差境界と性質の保証を維持する。
• 得られたスキームをDynIbex検証済みシミュレーションツールに実装し、ベンチマーク問題における実験的評価を実施する。
• RadauIIA、Gauss-Legendre、RK4といった古典的手法と比較することで、解の包含範囲と計算効率の両面から結果を検証する。

実験結果

リサーチクエスチョン

• RQ1 区間係数を用いたRunge-Kutta法が、浮動小数点算術では失われがちな次数、安定性、シンプレクティック性といった本質的性質を保持できるか。
• RQ2 区間係数を用いたButcherテーブル条件から生じる高次元的かつ未定義の高い方程式系に対して、制約プログラミングが効果的に適用可能か。
• RQ3 次に高い次数の方法に対する局所切り捨て誤差（LTE）を最適化することで、精度と効率の両面で競争力のあるRunge-Kuttaスキームが得られるか。
• RQ4 小さな係数の区間ですら、区間係数スキームが浮動小数点法のそれよりも解の包含範囲を狭くできるか。
• RQ5 パrameterの不確実性や剛性のあるダイナミクスの存在下で、区間ベースの手法が収縮性とロバストネスをどの程度向上させるか。

主な発見

• 直径が1×10⁻¹⁰程度の区間係数を有するS3O4法（3次スキーム）は、5次精度のRadauIIA法と同等の性能を示したが、計算がはるかに高速で、高い効率性を示した。
• S3O5法（3次スキーム）は、6次精度のGauss-Legendre法よりも5倍速く、最終的な解の包含範囲も顕著に狭かった。
• ERK33法（明示的3次スキーム）は、RK4と同等の時間消費量であったが、解の包含範囲はRK4の1.2倍にとどまり、強力な競争力が示された。
• 直径が1×10⁻¹⁰程度の区間係数でさえ、古典的手法の正確な係数を用いた場合よりも、最終的な解の包含範囲が狭くなることが多く、収縮性の向上が示唆された。
• 制約プログラミングアプローチにより、係数が区間として表現されても、次数および安定性の性質を保証する新たなRunge-Kuttaスキームが成功裏に生成された。
• 最適化プロセスによりLTEが効果的に最小化され、高次精度の方法に近い性能を示すスキームが得られ、コスト関数の有効性が検証された。