[논문 리뷰] Computation of the incomplete gamma function for negative values of the argument

Amparo Gil, Diego Ruiz‐Antolín외 2명|arXiv (Cornell University)|2016. 08. 01.

Control Systems and Identification인용 수 1

한 줄 요약

이 논문은 실수 매개수 a와 음의 실수 인자 z에 대해 정규화된 불완전 감마 함수 γ*(a,z)를 계산하는 Fortran 90 모듈 IncgamNEG를 제시한다. 알고리즘은 시리즈 전개, Poincaré형 전개, 균일 점근 전개 및 재귀 관계를 조합하여 매개수 영역 (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) 영역에서 넓은 매개수 범위에서 높은 상대 정밀도를 확보하는 것.
• 다양한 수학적 기법을 매개수 영역에 따라 자동으로 전환하는 통합 알고리즘을 구현하는 것.

제안 방법

• 작은 |z|와 중간 정도의 |a|에 대해 시리즈 전개를 사용한다.
• 중간 매개수 영역에서는 수렴성을 향상시키기 위해 Poincaré형 전개를 적용한다.
• 큰 |a|와 |z|에 대해서는 정확도를 유지하기 위해 균일 점근 전개를 활용한다.
• 매개수 전환과 수치적 안정성을 향상시키기 위해 재귀 관계를 통합한다.
• a와 z의 값에 따라 적절한 계산 기법을 동적으로 선택하는 방법을 사용한다.
• 포괄성과 재사용성을 확보하기 위해 IncgamNEG라는 이름의 Fortran 90 모듈에 구현한다.

실험 결과

연구 질문

• RQ1 어떻게 하면 음의 실수 값 z에 대해 불완전 감마 함수 γ*(a,z)를 정확하게 계산할 수 있는가?
• RQ2 어떤 점근 전개와 시리즈 방법의 조합이 넓은 매개수 범위에서 높은 정밀도를 보장하는가?
• RQ3 재귀 관계는 어떻게 효과적으로 통합되어 혼합 매개수 영역에서 수치적 안정성을 유지하는가?
• RQ4 이러한 하이브리드 알고리즘은 (a,z) ∈ [−500,500] × [−500,0) 영역에서 달성 가능한 상대 정밀도는 얼마인가?
• RQ5 단일의 통합 알고리즘이 a와 z의 다양한 영역에서 γ*(a,z)의 다양한 수학적 행동을 효율적으로 처리할 수 있는가?

주요 결과

• 알고리즘은 매개수 영역 (a,z) ∈ [−500,500] × [−500,0) 전역에서 약 10⁻¹³의 상대 정밀도를 달성한다.
• 시리즈, Poincaré형, 균일 점근 전개를 조합한 하이브리드 접근 방식이 다양한 매개수 영역에서 강건성을 확보한다.
• 재귀 관계의 사용은 특히 중요한 전이 영역 근처에서 수치적 안정성을 향상시킨다.
• Fortran 90 모듈 IncgamNEG는 일관된 성능을 보이며 알고리즘을 성공적으로 구현한다.
• 이 방법은 불완전 감마 함수에서 음의 인자와 관련된 계산적 과제를 효과적으로 처리한다.