[論文レビュー] NumGfun: a Package for Numerical and Analytic Computation with D-finite Functions

Marc Mezzarobba|arXiv (Cornell University)|Feb 16, 2010

Numerical Methods and Algorithms被引用数 5

ひとこと要約

この論文では、多項式係数をもつ線形微分方程式の解であるD-有限関数に対する高速かつ証明可能な数値的・解析的計算を実装する、Maple用パッケージNumGfunを提示する。チャドノフスキー＆チャドノフスキーのビットバースト法を含む高度なアルゴリズムを用い、組合せ的数列や特殊関数の高精度評価を、特異点付近でも可能にし、ソフト線形時間計算量と厳密な誤差境界を達成する。

ABSTRACT

This article describes the implementation in the software package NumGfun of classical algorithms that operate on solutions of linear differential equations or recurrence relations with polynomial coefficients, including what seems to be the first general implementation of the fast high-precision numerical evaluation algorithms of Chudnovsky & Chudnovsky. In some cases, our descriptions contain improvements over existing algorithms. We also provide references to relevant ideas not currently used in NumGfun.

研究の動機と目的

• D-有限関数評価の理論的高速アルゴリズムと、コンピュータ代数システムにおける実用的ソフトウェア実装との間のギャップを埋めること。
• D-有限関数の一般用途で、特異点（特に不規則特異点を含む）や高精度での証明可能な数値的評価フレームワークを提供すること。
• チャドノフスキー＆チャドノフスキーのビットバーストアルゴリズムを完全な誤差制御とともに実装・最適化し、ソフト線形時間評価を可能にすること。
• 特に組合せ論および特殊関数において、D-有限関数の記号的計算と数値的計算の両方をサポートすること。
• 形式的正しさの保証を得た再利用可能かつ拡張可能なソフトウェアパッケージをMapleに統合すること。

提案手法

• パッケージは、D-有限関数の高速数値評価にビットバースト法を用い、nビットの精度に対してO(n log^{O(1)} n)のビット複雑度を達成する。
• 微分方程式および初期条件から導かれる厳密な誤差境界を統合し、証明可能な結果を保証する。
• 特異点付近の評価には、数値的解析接続とメジャントリクセ方程式を用いた遷移行列の境界を適用する。
• 基本解の切断されたテイラー級数による多項式近似を用い、誤差制御された展開次数を採用する。
• diffeqtoproc関数は、低次の多項式を事前計算し誤差を有界にすることで、ディスク上での証明可能な評価手続きを生成する。
• 2段階戦略を採用する：基準点から必要な展開次数を推定し、必要に応じて精度を向上させる。

実験結果

リサーチクエスチョン

• RQ1 チャドノフスキー＆チャドノフスキーのビットバーストアルゴリズムは、完全な誤差制御を伴う一般用途のコンピュータ代数システムに実装可能か？
• RQ2 不規則特異点を含む特異点におけるD-有限関数の証明可能な高精度評価は、どのように達成できるか？
• RQ3 繰り返し評価に適したディスク上での多項式近似の性能と精度のトレードオフは何か？
• RQ4 D-有限関数の数値的評価を、組合せ的数列および特殊関数の両方に対して効率的かつ信頼性高く実現できるか？
• RQ5 任意の多項式係数をもつ線形微分方程式によって定義される関数に対しても、アルゴリズムをどのように適合させられるか？

主な発見

• 10万番目のモーツキン数は4.7秒で計算可能であり、大きなインデックスでは単純な再帰的展開を上回る性能を示す。
• z = -0.99における二重特異的ヘン関数を1000桁の精度で22秒で評価し、特異点付近での高精度評価能力を実証した。
• M₁₀⁶の評価には1分を要し、約477,112桁の数が得られ、スケーラビリティを確認した。
• diffeqtoproc関数は、繰り返し高価な評価を回避する証明可能な評価手続きを生成し、全プロットの合計時間は約9秒にまで短縮された。
• 形式的境界（例：bound_diffeq）と数値的評価（例：evaldiffeq）の両方を、厳密な誤差制御でサポートする。
• 高精度評価においてソフト線形時間計算量を達成し、チャドノフスキー＆チャドノフスキーおよびヴァン・ダー・ホーヴェンの理論的境界と一致した。