[論文レビュー] General Complex Polynomial Root Solver and Its Further Optimization for Binary Microlenses

J. Skowron, Andrew Gould|arXiv (Cornell University)|Mar 5, 2012

Astronomy and Astrophysical Research参考文献 2被引用数 24

ひとこと要約

本稿では、判別式に基づいてラウージェル法、ニュートン法、および新規の中間手法を動的に選択する新しい複素多項式根の求解法を提示する。商用ソルバー（ZROOTSなど）と比較して1.6–3倍の高速化を達成した。本手法は、二重重力レンズの問題に最適化されており、問題固有のヒューリスティクスと失敗に強い数値的安定性を用いて、根の特定処理を高度化することで計算時間を短縮する。

ABSTRACT

We present a new algorithm to solve polynomial equations, and publish its code, which is 1.6-3 times faster than the ZROOTS subroutine that is commercially available from Numerical Recipes, depending on application. The largest improvement, when compared to naive solvers, comes from a fail-safe procedure that permits us to skip the majority of the calculations in the great majority of cases, without risking catastrophic failure in the few cases that these are actually required. Second, we identify a discriminant that enables a rational choice between Laguerre's Method and Newton's Method (or a new intermediate method) on a case-by-case basis. We briefly review the history of root solving and demonstrate that "Newton's Method" was discovered neither by Newton (1671) nor by Raphson (1690), but only by Simpson (1740). Some of the arguments leading to this conclusion were first given by the British historian of science Nick Kollerstrom in 1992, but these do not appear to have penetrated the astronomical community. Finally, we argue that Numerical Recipes should voluntarily surrender its copyright protection for non-profit applications, despite the fact that, in this particular case, such protection was the major stimulant for developing our improved algorithm.

研究の動機と目的

• 天文的応用、特に二重重力レンズを対象として、より高速で数値的に安定した複素多項式根の求解法を開発すること。
• 重力レンズの光曲線をモデリングする際、根の特定処理が処理時間の最大80％を占める計算コストを低減すること。
• ラウージェル法、ニュートン法、および新規の中間根の特定法の間を知的に選択することで、数値的安定性と効率性を向上させること。
• 二重レンズ方程式に現れる5次複素多項式に特化して、求解法を最適化すること。
• Numerical Recipesなどの既存の商用ソルバーにおけるアクセス性と透明性の欠如に対応し、許可の緩いライセンスでオープンソースコードを公開すること。

提案手法

• 大多数のケースで計算を飛ばしてもエッジケースで失敗しない、フェイルセーフな手順を採用している。
• 各根ごとに最適な根の特定法（ラウージェル法、中間法、またはニュートン法）を選択するための判別式に基づく意思決定ルールを導入している。
• 二段階アプローチを採用：まず、各根が特定された後、逐次除法を用いて根を安定的に特定する。次に、全多項式を用いて根を精密化する。
• 収束の挙動に応じて、ラウージェル法、新規の中間法（SG）、およびニュートン法の間を動的に遷移する新しい動的法（CMPLX_LAGUERRE2NEWTON）を実装している。
• 二重重力レンズに特化した最適化を適用しており、例えば根の分離度に基づくソーティングや、2次および3次部分問題に対して解析的ソルバーを用いるなどしている。
• 収束閾値と反復回数の上限に基づく停止基準を採用し、速度と正確性のバランスを取っている。

実験結果

リサーチクエスチョン

• RQ1 数値的安定性を損なわずに、一般の複素多項式根の求解法を著しく高速化することは可能か？
• RQ2 ラウージェル法、ニュートン法、および中間法の間で最適な性能を得るために、信頼性の高い判別式は何か？
• RQ3 二重レンズ方程式の構造をどのように活用することで、一般の多項式求解法を上回る根の特定最適化が可能か？
• RQ4 レンズ方程式求解法における数値精度の実用的限界は何か？そして、これらは天文学者の間で既知であろうか？
• RQ5 Numerical Recipesのような商用数値ライブラリが、非営利の学術的利用を目的として著作権を自主的に放棄すべきであろうか？

主な発見

• 新規ソルバーは、商業的に利用可能なZROOTSサブルーチンと比較して、アプリケーションに応じて1.6–3倍の高速化を達成した。
• 判別式に基づく動的メソッド選択により、非効率な根の特定戦略を回避し、性能の向上が顕著に見られた。
• フェイルセーフな手順により、大多数のケースで計算を飛ばすことが可能となり、信頼性を損なわず高速化に寄与した。
• 特定の5次複素多項式の構造を活用することで、二重重力レンズの文脈において、一般の多項式求解法を上回る顕著な性能向上が達成された。
• 研究者たちは、レンズ求解法における数値精度の限界を同定したが、これらは天文学者の間では未知のものであるように思われるが、数値解析の分野では既知かもしれない。
• オープンソースコードは許可の緩いライセンスで公開されており、学術的利用と透明性を促進している。著者らは、数値ソフトウェアにおいても同様の著作権放棄を推奨している。