[論文レビュー] Two-dimensional generalization of the Muller root-finding algorithm and its applications

Plamen Fiziev, Denitsa Staicova|arXiv (Cornell University)|May 28, 2010

Model Reduction and Neural Networks参考文献 18被引用数 7

ひとこと要約

本稿では、複素変数を含む2つの非線形超越方程式の連立系を解くための、ミュラー・ルート・ファインディング・アルゴリズムの2次元一般化を提案する。特に、共形ヘン関数を含む問題に対して特に有効である。導関数の評価に失敗する、または収束が遅いような困難なケース—例えばシュワルツシルト黑洞やカー黑洞の準正規モード計算—において、ニュートン法やブロイデン法よりも高速な収束を達成する。

ABSTRACT

We propose a new algorithm for solving a system of two nonlinear transcendental equations with two complex variables based on the Muller algorithm. The two-dimensional Muller algorithm is tested on systems of different type and is found to work comparably to Newton's method and Broyden's method in many cases. The new algorithm is particularly useful in systems featuring the Heun functions whose complexity may make the already known algorithms not efficient enough or not working at all. In those specific cases, the new algorithm gives distinctly better results than the other two methods. As an example for its application in physics, the new algorithm was used to find the quasi-normal modes (QNM) of Schwarzschild black hole described by the Regge-Wheeler equation. The numerical results obtained by our method are compared with the already published QNM frequencies and are found to coincide to a great extent with them. Also discussed are the QNM of the Kerr black hole, described by the Teukolsky Master equation.

研究の動機と目的

• 特殊関数（ヘン関数など）を含む複素非線形系を解く際のニュートン法やブロイデン法の限界を克服する。
• 共形ヘン関数を含む系において、導関数の評価が困難で収束が悪いという計算上の課題に対処する。
• 微分を必要とせず、高速に収束するアルゴリズムを、理論物理学のスペクトル問題（例えばブラックホールの準正規モード）に適した形で開発する。
• 標準的な反復法が失敗する、または収束が遅すぎるような系に対して、信頼性の高い数値的代替手法を提供する。
• レッジ＝ウェーラー方程式およびテウコルスキー方程式に記述されるシュワルツシルトおよびカー黑洞の準正規モードに対する、実際の物理問題へのアルゴリズムの有効性を実証する。

提案手法

• 1次元ミュラー法（元来は2次セカント法で、超線形収束（約1.84）を示す）を2つの複素変数へ一般化する。
• 3つの初期点に基づく2次元の補間を用いて根を近似し、明示的な導関数計算を回避する。
• 1次元ミュラー部分計算を繰り返し適用する2次元フレームワーク内で統合し、根の推定値を改善する。
• 1次元部分計算における最大反復回数 $ P $ を調整することで収束をチューニングし、全体の性能と収束速度に影響を与える。
• 2つの複素方程式の連立系にアルゴリズムを適用し、収束の信頼性を高めるために関数の順序付けや線形結合の工夫に特に注意を払う。
• ヘン関数を含む系を用いた数値的テストで本手法を実装し、ニュートン法やブロイデン法と結果を比較する。

実験結果

リサーチクエスチョン

• RQ1 微分を必要とせず、2次元の根の探査アルゴリズムが、複雑で多項式でない関数を含む系において、ニュートン法やブロイデン法を上回る性能を示せるか？
• RQ2 導関数の評価が困難な共形ヘン関数を含む系において、2次元ミュラー法はどのように性能を発揮するか？
• RQ3 初期推定値と関数の順序付けが、新しいアルゴリズムにおける収束行動および根の選択に与える影響は何か？
• RQ4 1次元ミュラー部分計算における反復回数 $ P $ のチューニングが、2次元手法の収束速度と精度に顕著な影響を与えるか？
• RQ5 本手法は、レッジ＝ウェーラー方程式およびテウコルスキー方程式で記述されるシュワルツシルトおよびカー黑洞の準正規モード周波数を信頼性高く計算できるか？

主な発見

• 共形ヘン関数を含む系において、導関数の評価に失敗する、または不正確であるようなケースにおいて、2次元ミュラー法はニュートン法やブロイデン法よりも高速に収束する。
• シュワルツシルト黑洞の準正規モードに関して、本手法で得られた周波数は、既に発表済みの結果とよく一致しており、その正確性が裏付けられる。
• 一部のケースでは、1次元部分計算における最大反復回数 $ P $ を増加させることで、合計計算時間を短縮できた。これは $ P $ に対して非単調な依存関係があることを示唆している。
• アルゴリズムの性能は関数の順序付けと初期条件に敏感である。関数の再順序付けや方程式の線形結合を用いることで、発散や望ましくない根への収束を防げる。
• ヘン関数を含む系において、ブロイデン法は一般にロバストであるが、本手法は収束速度においてブロイデン法を上回る。
• 実数解や単純な非線形性を含む系では、2次元ミュラー法はニュートン法やブロイデン法と同等の性能を示すが、特殊関数を含む複雑で特異なケースにおいて顕著な優位性を示す。