[論文レビュー] Convergence Acceleration Techniques

U. D. Jentschura, Sergey Aksenov+3 more|ArXiv.org|Feb 1, 2002

Embedded Systems Design Techniques被引用数 3

ひとこと要約

本稿では、非交代級数を交代級数に変換するヴァン・ヴァインガーデン変換を最初に適用し、その後デルタ変換を適用してさらに収束を加速する、組み合わせ非線形凝縮変換（CNCT）を導入する。この手法により、収束が遅いまたは発散する級数の収束が高速化され、200桁の精度を数時間で達成可能となり、10^205項が必要な困難なゼータ級数の和に対しては約84,000項で実現可能である。

ABSTRACT

This work describes numerical methods that are useful in many areas: examples include statistical modelling (bioinformatics, computational biology), theoretical physics, and even pure mathematics. The methods are primarily useful for the acceleration of slowly convergent and the summation of divergent series that are ubiquitous in relevant applications. The computing time is reduced in many cases by orders of magnitude.

研究の動機と目的

• 物理学、生物学、数学に共通する収束が遅いまたは発散する級数を効果的に加速する堅牢な数値的手法の開発。
• 非交代級数に対する標準的な収束加速法における数値的不安定性の克服。
• 最小限の計算コストで特殊関数（多対数関数、ゼータ関数、レルヒトランスェンデントなど）の高精度評価を可能にする。
• 量子物理学やバイオインフォマティクスを含む科学計算分野に応用可能な汎用的ツールの提供。
• 極めて高い精度を要する級数評価において、計算作業量を桁単位で削減できるかの実証。

提案手法

• 非交代級数を交代級数に変換するため、ヴァン・ヴァインガーデン変換を適用し、切り捨て誤差の影響を低減する。
• 変換後の級数にデルタ変換を適用して、さらなる収束加速を実現する。
• 組み合わせた変換は、部分和 $ s_n $ を変換値 $ \\_n}^{(0)}(1,\mathbf{S}_n) $ に写像する「組み合わせ非線形凝縮変換（CNCT）」と呼ばれる。
• 部分和の末尾桁に隠された情報を用いて、非線形補外により切り捨て誤差をモデル化・除去する。
• C言語、Fortran、Mathematicaでの実装が行われ、サンプルコードはダウンロード可能である。
• 特に $ |z| \approx 1 $ の場合、例えば多対数関数やレルヒトランスェンデントにおいて、標準的手法が失敗するが、本手法は特に有効である。

実験結果

リサーチクエスチョン

• RQ1 収束が遅い非交代級数に対して、高い数値的安定性を保ちながら収束加速技術を一般化できるか。
• RQ2 CNCTは、多対数関数やゼータ関数といった特殊関数の評価において、計算コストをどの程度削減できるか。
• RQ3 84,000項しか使用しない状況で、10^205項が必要な級数をどの程度の精度で計算できるか。
• RQ4 CNCTは、量子電磁力学計算において10^-15以上の精度を達成できるか。
• RQ5 $ \sum_{k=1}^\infty \left( H_k \right)^2 k^{-2} $ のような級数に対して、直接和を取るのと比較してCNCTがどの程度の性能向上を示すか。

主な発見

• CNCTは、12回の変換順序で $ 10^{-1} \mathrm{Li}_3(0.99999) $ の計算において15桁の相対精度を達成した。これは直接的に約1億項を和算するのと同等の精度である。
• $ \sum_{k=1}^\infty \left( H_k \right)^2 k^{-2} $ の級数に対して、CNCTは84,000項で200桁の精度を数時間で達成し、計算コストを約200桁単位で削減した。
• 水素原子の4P状態におけるベーテ対数の計算において、29桁の精度を達成した。これは、以前の最高精度結果よりも9桁も高い精度であった。
• 標準的手法が収束が遅いため失敗する $ |z| \approx 1 $ の $ \Phi(z,s,v) $ に対しても、CNCTは収束を達成した。
• 量子電磁力学の束縛状態計算において、変換が収束を安定化・加速させ、相対精度10^-14レベルの測定を可能にした。
• 数学、物理学、計算生物学の分野における特殊関数（ゼータ関数、多対数関数、レルヒトランスェンデント関数など）の評価において、本手法は堅牢性と高い精度を示した。