[論文レビュー] Computation of the incomplete gamma function for negative values of the argument

Amparo Gil, Diego Ruiz‐Antolín+2 more|arXiv (Cornell University)|Aug 1, 2016

Control Systems and Identification被引用数 1

ひとこと要約

本稿では、実数パラメータ a および負の実数引数 z に対して正規化された不完全ガンマ関数 γ*(a,z) を計算する Fortran 90 モジュール IncgamNEG を提示する。アルゴリズムは、級数展開、ポアンカレ型展開、一様漸近展開、再帰関係を組み合わせることで、パrameter領域 (a,z) ∈ [−500,500] × [−500,0] において約 10⁻¹³ の相対精度を達成する。

ABSTRACT

An algorithm for computing the incomplete gamma function $\gamma^*(a,z)$ for real values of the parameter $a$ and negative real values of the argument $z$ is presented. The algorithm combines the use of series expansions, Poincar\'e-type expansions, uniform asymptotic expansions and recurrence relations, depending on the parameter region. A relative accuracy $\sim 10^{-13}$ in the parameter region $(a,z) \in [-500,\,500] imes [-500,\,0)$ can be obtained when computing the function $\gamma^*(a,z)$ with the Fortran 90 module {\bf IncgamNEG} implementing the algorithm.

研究の動機と目的

• 引数 z が負のときの不完全ガンマ関数 γ*(a,z) の堅牢な計算手法を開発すること。
• 負の実数 z の領域における γ*(a,z) の計算における数値的不安定性と、効率的なアルゴリズムの欠如に対処すること。
• 特に (a,z) ∈ [−500,500] × [−500,0) の広いパrameter範囲にわたり、高い相対精度を保証すること。
• パrameter領域に応じて異なる数学的技法を適応的に切り替える統合アルゴリズムを実装すること。

提案手法

• 小規模な |z| および中程度の |a| の領域では級数展開を用いる。
• 中間的パrameter領域では、収束性を向上させるためにポアンカレ型展開を適用する。
• 大規模な |a| および |z| の領域では、一様漸近展開を用いて精度を維持する。
• パrameterの遷移を処理し、数値的安定性を向上させるために再帰関係を統合する。
• a および z の値に応じて、適切な計算手法を動的に選択する。
• 実装は、移植性および再利用性を考慮し、IncgamNEG という名前の Fortran 90 モジュールにカプセル化されている。

実験結果

リサーチクエスチョン

• RQ1 負の実数 z に対して不完全ガンマ関数 γ*(a,z) をどのように高精度に計算できるか？
• RQ2 どの組み合わせの漸近展開と級数法が、広いパrameter範囲にわたり高い精度を保証するか？
• RQ3 再帰関係をどのように効果的に統合することで、混合パrameter領域における数値的安定性を維持できるか？
• RQ4 このようなハイブリッドアルゴリズムが、領域 (a,z) ∈ [−500,500] × [−500,0) において達成可能な相対精度は何か？
• RQ5 一様なアルゴリズムが、a および z の異なる領域における γ*(a,z) の多様な数学的挙動を効率的に処理できるか？

主な発見

• アルゴリズムはパrameter領域 (a,z) ∈ [−500,500] × [−500,0) 全体で約 10⁻¹³ の相対精度を達成する。
• 級数展開、ポアンカレ型展開、一様漸近展開を組み合わせたハイブリッド手法により、さまざまなパrameter領域でも堅牢性が保証される。
• 再帰関係の使用により、特に重要な遷移領域近くでの数値的安定性が向上する。
• Fortran 90 モジュール IncgamNEG は、一貫したパフォーマンスを発揮する形でアルゴリズムを効果的に実装している。
• 本手法は、不完全ガンマ関数における負の引数に起因する計算上の課題を効果的に処理している。