[논문 리뷰] The High Precision Numerical Calculation of Stieltjes Constants. Simple and Fast Algorithm

K. Maślanka, Andrzej Koleżyński|arXiv (Cornell University)|2022. 06. 22.

Analytic Number Theory Research인용 수 1

한 줄 요약

이 논문은 리만 제타 함수의 초함수 유사 전개와 등간격 실수 변수에서의 고정밀도 ζ(1 + jε) 평가를 활용하여, 초고정밀도—최대 80,000자리 정밀도까지—로 스틸티jes 상수를 계산하는 매우 효율적이고 수학적으로 단순한 알고리즘을 제시한다. 이 방법의 성공은 이전에 초고정밀도 수준에서 오류를 유발했던 PARI/GP의 제타 함수 계산에서 발생하는 핵심 버그를 해결한 데 기반한다.

ABSTRACT

We present a simple but efficient method of calculating Stieltjes constants at a very high level of precision, up to about 80000 significant digits. This method is based on the hypergeometric-like expansion for the Riemann zeta function presented by one of the authors in 1997 \cite{Maslanka 1}. The crucial ingredient in this method is a sequence of high-precision numerical values of the Riemann zeta function computed in equally spaced real arguments, i.e. $ζ(1+\varepsilon),ζ(1+2\varepsilon),ζ(1+3\varepsilon),...$ where $\varepsilon$ is some real parameter. (Practical choice of $\varepsilon$ is described in the main text.) Such values of zeta may be readily obtained using the PARI/GP program, which is especially suitable for this.

연구 동기 및 목표

• 초고정밀도에서 스틸티제스 상수를 수치적으로 효율적이고 수학적으로 단순한 방법으로 계산하는 데 목적이 있다.
• 고정밀도를 확보하기 위해 천문학적 합계가 필요한 고전적 정의의 극도로 느린 수렴 문제를 극복하는 데 목적이 있다.
• 이전에는 계산 및 수치적 불안정성으로 인해 접근이 어려웠던 80,000자리 정밀도까지의 스틸티제스 상수를 신뢰성 있게 계산할 수 있도록 하는 데 목적이 있다.
• 초고정밀도에서 결과를 손상시켰던 PARI/GP의 제타 함수 구현에 존재하는 핵심 버그를 식별하고 수정하는 데 목적이 있다. 이는 알고리즘의 신뢰성에 직접적인 영향을 미쳤다.
• 미래 연구에서 더 정확한 스틸티제스 상수의 점근적 전개를 유도하기 위한 견고한 기초를 제공하는 데 목적이 있다.

제안 방법

• 이 방법은 저자 중 한 명이 1997년에 개발한 리만 제타 함수의 초함수 유사 전개를 핵심 분석적 프레임워크로 사용한다.
• PARI/GP 컴퓨터 대수 시스템을 활용하여 j = 1, 2, ... 에 대해 실수 변수 ζ(1 + jε)에서의 고정밀도 제타 함수 값 계산을 수행한다.
• 알고리즘은 유한 차분 기법과 급수 역행렬 기법을 통해 스틸티제스 상수를 추출하기 위해 정규화된 제타 값 시퀀스 ζ(1 + jε)에 의존한다.
• 제타 함수 평가의 정밀도가 핵심적이다. 저자들은 정확성을 확보하기 위해 PARI/GP를 커스터마이징한 설정을 사용했으며, 여러 컴파일러, 운영 체제, 하드웨어 플랫폼을 테스트했다.
• PARI/GP의 여러 환경과 버전에서 광범위한 테스트를 수행하여, 오일러-매클라린 알고리즘 내에서 베르누이 수 계산 시 발생하는 반올림 오류를 규명하고 분리했다.
• 베르누이 수 캐싱 중 이중 반올림이 발생하는 것으로 밝혀진 PARI/GP 버그의 수정은 74,000자리 및 80,000자리 정밀도에서 신뢰할 수 있는 결과를 도출하기 위해 필수적이었다.

실험 결과

연구 질문

• RQ1 간단하고 효율적인 알고리즘을 사용하여 스틸티제스 상수를 80,000자리 정밀도까지 안정적으로 계산할 수 있는가?
• RQ2 고정밀도 스틸티제스 상수 계산에서 주요 수치적 장애물은 무엇이며, 이를 어떻게 극복할 수 있는가?
• RQ3 특히 ζ(1 + jε)에 대해 PARI/GP의 제타 함수 구현이 초고정밀도 수준에서도 정확한가?
• RQ4 초고정밀도 제타 함수 평가에서 체계적인 오류는 무엇을 원인으로 하고 있으며, 이를 어떻게 추적하고 수정할 수 있는가?
• RQ5 수정된 제타 함수 데이터는 향후 더 정확한 스틸티제스 상수의 점근적 전개를 도출하는 데 기여할 수 있는가?

주요 결과

• 알고리즘이 성공적으로 스틸티제스 상수를 80,000자리 정밀도까지 계산하여, 이 종류의 상수에 대해 이전에 없었던 초고정밀도를 달성했다.
• 이 방법의 효율성과 단순성은 등간격 실수 변수에서의 정규화된 제타 함수 평가와 초함수 유사 전개의 사용에 기반한다.
• 베르누이 수 캐싱 중 이중 반올림이 발생하는 것으로 밝혀진 PARI/GP의 제타 함수 계산에서의 핵심 버그가 식별되고 수정되었으며, 이는 이전에 74,000자리 및 80,000자리 정밀도에서 결과를 손상시켰다.
• 오류는 초고정밀도에서만 나타났으며, 낮은 정밀도 계산에서는 탐지되지 않아, 초고정밀 산술에서 철저한 테스트의 필요성을 강조한다.
• 수정된 PARI/GP 구현은 신뢰할 수 있는 제타 함수 데이터를 회복시켰으며, 이는 알고리즘의 성공에 필수적이었다.
• 결과는 알고리즘의 견고성을 검증하며, 향후 연구에서 더 정확한 스틸티제스 상수의 점근적 전개를 도출할 수 있는 길을 열었다.