[論文レビュー] New Progress in Univariate Polynomial Root-finding

Victor Y. Pan|arXiv (Cornell University)|May 30, 2018

Advanced Numerical Analysis Techniques被引用数 2

ひとこと要約

本稿では、未定義の多項式根の求め方を、新しい技法を用いて分割とEhrlichの関数的反復法を強化することで高速化し、特にスパースな多項式において顕著な高速化を達成している。また、パスフォロwingニュートン法の反復処理の改善も行われている。これらの改善により、既存のツール（例：MPSolve）と同等またはそれ以上の実用性と競争力を持つようになっている。

ABSTRACT

Univariate polynomial root-finding has been studied for four millennia and is still the subject of intensive research. Hundreds of efficient algorithms for this task have been proposed. Two of them are nearly optimal. The first one was proposed in 1995; it relies on recursive factorization of a polynomial, is quite involved, and has never been implemented. The second one was proposed in 2016, relies on subdivision iterations, was implemented in 2018, and promises to be practically competitive, although user's current choice for univariate polynomial root-finding is the package MPSolve, proposed in 2000, revised in 2014, and based on Ehrlich's functional iterations. By incorporating some old and new techniques we significantly accelerate subdivision and Ehrlich's iterations; as by-product we also accelerate path-following Newton's iterations for root-finding.. Moreover our acceleration of the known subdivision root-finders is dramatic in the case of sparse input polynomials. Some of our techniques promise to be valuable for the design and analysis of other polynomial root-finders as well.

研究の動機と目的

• 4,000年以上にわたり研究されてきた、一変数多項式の根を効率的に求めるという長年の課題に取り組む。
• 特にスパース多項式において、実装の欠如や実用的な効率性の不足といった、既存の根の求め方の限界を克服すること。
• 分割法とEhrlichの関数的反復法といった既知の反復的手法の性能を、新しいアルゴリズム的技法によって向上させること。
• 提案された改善の副産物として、パスフォロウイングニュートン法の反復処理を高速化すること。
• 現在の手法に有効な技術を、将来的な多項式根の求め方のアルゴリズム設計に対しても応用可能な可能性を秘めた技術を開発すること。

提案手法

• スパース多項式において特に顕著な高速化を実現するため、分割に基づく根の求め方を加速するため、古くからの手法と新規手法を統合する。
• 加速技法を適用することで、Ehrlichの関数的反復法の収束速度を向上させること。
• スパース多項式の構造的性質を活用し、分割法における計算コストを顕著に削減すること。
• ある手法（例：Ehrlichの手法）からの加速戦略を他手法（例：パスフォロウイングニュートン法）へと応用すること。
• 数値的安定性を保ちつつ速度を向上させるアルゴリズムの微調整を実施すること。
• 再帰的因数分解と反復的精錬を基盤として用い、新規の加速メカニズムによってこれを強化すること。

実験結果

リサーチクエスチョン

• RQ1 分割に基づく根の求め方を、特にスパース多項式に対してどのように高速化できるか？
• RQ2 Ehrlichの関数的反復法は、収束性と実用的性能を向上させるために、どのような方法で改善できるか？
• RQ3 ある根の求め方のための加速技法を、他の手法（例：パスフォロウイングニュートン法）へ効果的に転用できるか？
• RQ4 多項式のスパarsityと構造的性質は、顕著な性能向上を実現するために果たす役割は何か？
• RQ5 古典的手法と現代的手法のどの組み合わせが、根の求め方のアルゴリズムにおいて最も効果的な高速化をもたらすか？

主な発見

• 分割に基づく根の求め方の高速化により、特にスパースな入力多項式において顕著な性能向上が達成された。
• 新規の加速技法を適用することで、Ehrlichの関数的反復法の速度が著しく向上した。
• 提案された手法の副産物として、パスフォロウイングニュートン法の反復処理も高速化された。
• 提案された技術は、既存のアルゴリズムに有効であるだけでなく、将来的な多項式根の求め方のアルゴリズム設計に対しても有用である可能性を秘めている。
• 改善されたアルゴリズムは、特定の状況下でMPSolveのような現在の業界標準ツールと同等またはそれ以上の性能を発揮するようになった。
• 古典的手法と新規手法の統合により、一変数多項式根の求め方において実用的かつ効率的なフレームワークが構築された。