[논문 리뷰] Rigorous high-precision computation of the Hurwitz zeta function and its derivatives

Fredrik Johansson|arXiv (Cornell University)|2013. 09. 11.

Analytic Number Theory Research인용 수 4

한 줄 요약

이 논문은 증명 가능한 정확한 오차 한계를 갖는 엄밀한 고정밀 알고리즘을 제시하며, 오일러-매클로린 공식을 사용하여 허위츠 제타 함수 ζ(s,a)와 그 도함수를 계산한다. 이 방법은 스티엘티제 상수, 키퍼-리 계수, 리만 제타 함수의 첫 번째 비자명한 영점에 대한 기록적인 계산을 가능하게 하며, γ₁₀₀₀₀₀에 대해 최대 10,860자리의 정확한 숫자를 확보하고 오차 제어를 검증한다.

ABSTRACT

We study the use of the Euler-Maclaurin formula to numerically evaluate the Hurwitz zeta function $ζ(s,a)$ for $s, a \in \mathbb{C}$, along with an arbitrary number of derivatives with respect to $s$, to arbitrary precision with rigorous error bounds. Techniques that lead to a fast implementation are discussed. We present new record computations of Stieltjes constants, Keiper-Li coefficients and the first nontrivial zero of the Riemann zeta function, obtained using an open source implementation of the algorithms described in this paper.

연구 동기 및 목표

• 임의의 복소수 s와 a에 대해, 증명 가능한 오차 한계를 갖는 수치적으로 엄밀한 고정밀도로 허위츠 제타 함수 ζ(s,a)와 그 도함수를 평가하는 방법을 개발한다.
• 기존의 오일러-매클로린 접근법을 일반화하여 임의의 복소수 s와 a, 그리고 s에 대한 임의의 수의 도함수를 다룰 수 있도록 확장한다.
• 스티엘티제 상수나 키퍼-리 계수와 같은 특수 상수를 이전에 보다 높은 정확도로 계산하고 정밀도를 검증할 수 있도록 한다.
• 고정밀 산술 및 대규모 계산에 적합한 오픈소스이자 효율적인 구현을 제공한다.
• 고정밀 평가에서 치명적인 취약성으로 이어질 수 있는 중간 계산의 불안정성과 치명적인 취소 현상을 극복한다.

제안 방법

• 허위츠 제타 함수에 오일러-매클로린 합공식을 적용하고, 적분 잔여항을 통한 해석적 계속을 수행한다.
• 주요 급수와 주기적 베르누이 다항식의 경계를 사용하여 잔여 적분 R에 대한 엄밀한 오차 경계를 유도한다.
• s를 형식적 거듭제곱 변수로 간주함으로써 s에 대한 도함수를 동시에 평가할 수 있도록 한다.
• s=1에서의 극을 처리하기 위해 발산하는 항 I를 형식적 거듭제곱 급수의 복합을 통해 해석적으로 계속한다.
• 꼬리 합에 대해 이진 분할을 사용하고, 빠른 다항식 산술을 통해 계산 비용을 감소시켜 효율적인 구현을 한다.
• 알고리즘은 고정밀 산술에 최적화되어 있으며, 각 항에 대해 정밀도를 동적으로 조정하여 메모리와 계산량을 줄인다.

실험 결과

연구 질문

• RQ1 오일러-매클로린 공식이 임의의 복소수 s와 a, 그리고 s에 대한 임의의 수의 도함수에 대해 허위츠 제타 함수와 그 도함수에 대해 증명 가능한 정확한 오차 한계를 제공하도록 어떻게 적응시킬 수 있는가?
• RQ2 고정밀도, 단일 점 평가(예: 수백에서 수천 자리)에 대해 ζ(s,a)의 고정밀 평가를 어떻게 효율적이고 확장 가능하게 만들 수 있는가?
• RQ3 스티엘티제 상수와 키퍼-리 계수의 기록적 계산에 적용했을 때 이 방법의 성능과 정확도는 어떠한가?
• RQ4 고정밀 산술에서의 속도, 메모리 사용량, 신뢰성 측면에서 기존 알고리즘과 비교해 볼 때 이 방법은 어떠한가?
• RQ5 이 방법을 사용하여 히ュ리스틱 추정치, 예를 들어 스티엘티제 상수에 대한 키프네슬-코피 공식의 점근적 근사치를 검증하거나 기각할 수 있는가?

주요 결과

• 이 방법은 n ≤ 100,000인 모든 스티엘티제 상수 γₙ을 최소 10,860자리의 정확도로 계산하였으며, γ₁₀₀₀₀₀은 1.991927306312541095658…×10⁸³⁴³²로 계산되었다.
• n ≤ 100,000인 모든 γₙ의 계산은 16스레드 시스템에서 26시간이 소요되었고, 피크 메모리는 80 GiB를 기록하여 대규모 표의 실현 가능성을 입증했다.
• 키프네슬-코피 근사 공식은 n ≤ 100,000 범위에서 n=137를 제외한 모든 경우에서 부호가 정확하다는 것이 검증되었으며, 엄밀한 경계를 통한 정확도를 확인했다.
• 이 구현은 n ≤ 1000인 모든 γₙ을 1000자리 정밀도로 11초 이내에 계산하며, 마테마티카 9.0보다 거의 4개 자리 수준 빠르게 작동한다.
• 이전 연구들이 히ュ리스틱 추정에 의존한 것과 달리, 이 방법은 검증 가능한 오차 경계를 제공하여 고정밀 산술에서의 정확성을 보장한다.