[논문 리뷰] Numerical Computations Concerning the GRH

David J. Platt|arXiv (Cornell University)|2013. 05. 14.

Analytic Number Theory Research참고 문헌 11인용 수 5

한 줄 요약

이 논문은 Dirichlet L-함수의 효율적이고 엄밀한 계산을 위한 두 가지 새로운 알고리즘을 제시하며, 모듈러스 $ q \leq 400,000 $인 모든 원시 Dirichlet L-함수에 대해 임의의 높이까지 일반화된 리만 가설(GRH)의 검증을 가능하게 한다. 고성능 컴퓨팅 클러스터와 FFT 및 간격 산술을 포함한 고급 수치 기법을 사용한 계산을 통해, GRH는 290억 개 이상의 L-함수에 대해 확인되었으며, $ q \leq 2,000,000 $인 모든 원시 특성에 대해 중심점에서의 영이 아님을 입증하였다.

ABSTRACT

We describe two new algorithms for the efficient and rigorous computation of Dirichlet L-functions and their use to verify the Generalised Riemann Hypothesis for all such L-functions associated with primitive characters of modulus q<=400,000. For even q, we check to height t_0=max(1e8/q,7.5e7/q+200) and for odd q to height t_0=max(1e8/q,3.75e7/q+200).

연구 동기 및 목표

• 모듈러스 $ q \leq 400,000 $인 모든 원시 Dirichlet L-함수에 대해 일반화된 리만 가설(GRH)을 엄밀히 검증하는 것.
• 이전 결과에서 $ q \leq 432 $ 및 낮은 높이에만 도달한 바를 초월하여 GRH의 계산적 검증을 확장하는 것.
• 효율적으로 Dirichlet L-함수와 그 영들을 고정밀도로 계산할 수 있는 두 가지 새로운 수치 알고리즘을 개발하고 구현하는 것.
• $ q \leq 2,000,000 $인 모든 원시 특성에 대해 $ L_{\chi}(1/2) \neq 0 $임을 확인하여 삼항 골드바흐 문제와 같은 해석수론 응용을 지원하는 것.
• 간격 산술과 고정밀도 오일러-매클라린 합산을 사용하여 함수방정식에서 발생하는 수치적 불안정성과 부호 모호성을 엄밀한 경계로 해결하는 것.

제안 방법

• 임의의 길이의 이산 푸리에 변환을 효율적으로 계산하기 위해 수정된 블루스타인 FFT 알고리즘을 사용하여, 근사 함수방정식을 통한 L-함수의 평가를 가능하게 하였다.
• 푸아송 합공식을 적용하여 정수에 대한 합을 푸리에 변환에 대한 합과 연결함으로써, $ L_{\chi}(1/2 + it) $의 이산점에서의 평가를 위한 샘플링 기반 방법을 구축하였다.
• 감마 인자와 지수 항으로 정규화된 완비된 $ \Lambda_{\chi}(t) $-함수를 사용하여 실수 성질과 대칭성을 확보함으로써, 부호 변화 탐지에 기여하였다.
• 두 단계 알고리즘을 구현: $ q \leq 10,000 $에 적합한 '작은 $ q $' 알고리즘과 고조도 모듈러스에 적합한 '큰 $ q $' 알고리즘으로, 메모리 및 캐시 효율성에 최적화되었으며 병렬 처리가 가능하도록 설계되었다.
• 간격 산술을 이중 정밀도 및 100비트 MPFI 정밀도로 적용하여 $ \Lambda_{\chi}(1/2 + it) $의 모호한 부호 변화를 해결하고, 영의 엄밀한 탐지 보장을 하였다.
• 영의 수를 추정하기 위해 터닝의 방법을 적용하고 실제 영 수와 교차 검증함으로써 격차를 수정하였으며, 적응형 샘플링과 고정밀 재계산을 통해 이를 수행하였다.

실험 결과

연구 질문

• RQ1 모듈러스 $ q \leq 400,000 $인 모든 원시 Dirichlet L-함수에 대해 GRH를 엄밀히 검증할 수 있으며, 임계대에서 어느 정도의 높이까지 검증 가능한가?
• RQ2 특히 큰 $ q $에 대해 고정밀도로 효율적이고 정확하게 Dirichlet L-함수와 그 영들을 계산할 수 있는 수치 알고리즘은 무엇인가?
• RQ3 간격 산술을 사용하여 함수방정식의 수치적 불안정성과 부호 모호성을 엄밀한 경계로 해결할 수 있는가?
• RQ4 $ q \leq 2,000,000 $인 모든 원시 Dirichlet L-함수에 대해 중심값 $ L_{\chi}(1/2) \neq 0 $임이 성립하는가, 그리고 이는 계산적으로 확인될 수 있는가?
• RQ5 현대 컴퓨팅 인프라와 알고리즘 발전을 통해 이전 작업의 범위를 얼마나 초월하여 GRH 검증의 범위를 확장할 수 있는가?

주요 결과

• 모듈러스 $ q \leq 400,000 $인 모든 원시 Dirichlet L-함수에 대해 GRH가 성립하며, 짝수 $ q $에 대해서는 높이 $ t_0 = \max\left(\frac{10^8}{q}, \frac{7.5 \cdot 10^7}{q} + 200\right) $까지, 홀수 $ q $에 대해서는 $ t_0 = \max\left(\frac{10^8}{q}, \frac{3.75 \cdot 10^7}{q} + 200\right) $까지 검증되었다.
• 모듈러스 $ q \leq 2,000,000 $인 모든 $ 739,151,526,102 $개의 원시 특성에 대해 $ L_{\chi}(1/2) $가 영이 아님을 확인하였으며, 이 중 20개만이 100비트 MPFI 정밀도를 사용해 해결이 필요하였다.
• 모듈러스 $ q \leq 400,000 $인 모든 경우에 대해 약 $ 3.8 \cdot 10^{13} $개의 영이 임계선 상에 고립되었으며, 선상 외부에 영이 발견되지 않았다.
• 작은 $ q $와 큰 $ q $에 최적화된 두 가지 새로운 알고리즘은 이전의 엄밀한 계산 범위보다 6개의 지수 차수 더 큰 범위에서 GRH의 검증을 가능하게 하였다.
• 간격 산술과 고정밀도 오일러-매클라린 합산은 $ \Lambda_{\chi}(1/2 + it) $의 부호 모호성을 성공적으로 해결하여, 수치적 불안정성으로 인한 영 누락을 방지하였다.
• 전체 계산은 영국 및 프랑스의 고성능 컴퓨팅 클러스터를 통해 약 40만 코어 시간을 소비하였으며, 두 알고리즘이 다중 코어 시스템에서 효율적으로 병렬 처리되었다.