[論文レビュー] Exactification of Stirling's Approximation for the Logarithm of the Gamma Function

Victor Kowalenko|arXiv (Cornell University)|Apr 10, 2014

Iterative Methods for Nonlinear Equations参考文献 22被引用数 3

ひとこと要約

本稿では、発散する漸近級数のボレルおよびメリン＝バーンズ（MB）正則化を用いて、ガンマ関数の対数のスターリング近似を厳密に修正する手法を提示する。両正則化手法は、切り捨てられた漸近展開と対数項を組み合わせた際に同一で、極めて高い精度の値をもたらし、数学的ソフトウェアと同等の精度に達する。MB正則化は、コンactな積分形を採用し、無限和の指数積分を回避するため、ボレル和よりも計算上優れていることが示された。

ABSTRACT

Exactification is the process of obtaining exact values of a function from its complete asymptotic expansion. Here Stirling's approximation for the logarithm of the gamma function or $\ln Γ(z)$ is derived completely whereby it is composed of the standard leading terms and an asymptotic series that is generally truncated. Nevertheless, to obtain values of $\ln Γ(z)$, the remainder must undergo regularization. Two regularization techniques are then applied: Borel summation and Mellin-Barnes (MB) regularization. The Borel-summed remainder possesses an infinite convergent sum of exponential integrals and discontinuous logarithmic terms across Stokes sectors and lines, while the MB-regularized remainders possess one MB integral, which is better to compute, and similar logarithmic terms. The MB integrals are valid over overlapping domains of convergence. Hence, two MB-regularized asymptotic forms can be used to evaluate $\ln Γ(z)$. Despite having to truncate the Borel-summed remainder, it is found that all the remainders combined with (1) the truncated asymptotic series, (2) the leading terms of Stirling's approximation and (3) their logarithmic terms yield identical values of $\ln Γ(z)$. In the few cases where the accuracy falls away, it is mostly due to a very high value for the truncation parameter, which results in the cancellation of redundant decimal places. Where possible, all the asymptotic forms yield the same values as the LogGamma routine in Mathematica.

研究の動機と目的

• 完全な漸近展開（すべての順位より下位の項および超越的級数を含む）から、ガンマ関数の対数の正確な値を求める手法を開発すること。
• 漸近展開における発散する剰余項の問題に対処し、それらを和が取れる形に正則化する手法を導入することにより、意味のある値を付与すること。
• ガンマ関数の対数を評価する際の、ボレル和とメリン＝バーンズ正則化の計算効率および数値的安定性を比較すること。
• 正則化が、lnΓ(z) のような多価関数を主分岐以外の領域にまで拡張して評価可能にする構造を保持できることを示すこと。
• ピカールの漸近論における伝統的信条に挑戦し、完全な漸近展開の微分が正当であり、正確な結果をもたらすことを示すこと。

提案手法

• 本稿では、漸近展開の発散する剰余項をボレル和により正則化し、発散する無限和の指数積分と、ストークス現象に起因する不連続な対数項の形で表現する。
• メリン＝バーンズ（MB）正則化を導入し、剰余項を1つのMB積分として表現する。これにより収束領域の重なりが可能となり、無限和を回避できる。
• コーシーの留数定理を用いて、発散級数をMB積分で表現し、発散する弧路積分を除外しながら、有限で正則化された値を保持する。
• ボレル和された剰余項は非常に大きな値（例：10^5）で切り捨てられ、高い精度を得る。一方、MB正則化形は、計算の簡潔さを実現するため、無限和をリーマン・ゼータ関数に置き換える。
• ストークス不連続性に起因する対数項（SD_M^{SL}(z)、SD_M^{SS}(z)、および SD_MB(M,z)）を明示的に含め、lnΓ(z) の多価性を捉え、主分岐 −π < arg z ≤ π を超えて評価可能にする。
• 切り捨てられた漸近級数、主要なスターリング項、および正則化された剰余項を組み合わせることで、高精度ソフトウェアの出力と一致する値が得られる。

実験結果

リサーチクエスチョン

• RQ1 スターリングの漸近展開における lnΓ(z) の発散する剰余項は、正則化によって正確な値に修正可能か？
• RQ2 ガンマ関数の対数を評価する際、ボレル和とメリン＝バーンズ正則化の精度および計算効率はどのように比較できるか？
• RQ3 正則化プロセスは、特に異なるリーマン面を越えて、lnΓ(z) の多価構造を保持できるか？
• RQ4 すべての順位より下位の項を含む完全な漸近展開は、適切に正則化された場合に、正確な値に還元可能か？
• RQ5 発散するが、正則化が施された完全な漸近展開に対して微分を行うことは許容可能か、かつ正確な結果をもたらすか？

主な発見

• ボレル和された剰余項とMB正則化された剰余項は、切り捨てられた漸近級数および主要なスターリング項と組み合わせた場合、高精度な数学的ソフトウェアの出力と同一の値をもたらす。
• MB正則化された漸近形態は、計算上優れており、無限和の指数積分を1つのMB積分とリーマン・ゼータ関数に置き換えることで、より高速かつ安定した評価が可能になる。
• ストークス不連続性に起因する対数項は、lnΓ(z) の多価性を捉えるために不可欠であり、主分岐 −π < arg z ≤ π を超えて評価を可能にする。
• ボレル和された剰余項は、MB正則化と同等の精度を得るためには非常に大きな値（例：10^5）で切り捨てられる必要があり、その結果、著しく遅く、非効率になる。
• 正則化プロセスは、無限大を有限部分（ハダマードの有限部に類似）として除去することで、発散級数を制御し、漸近解析におけるこのような手法の有効性を裏付けた。
• 本研究は、完全な漸近展開の微分が許容可能であり、正確な結果をもたらすことを示した。これは、ピカールの漸近論における長年の仮定に挑戦するものである。