A fixed-point iteration method for the number Pi with arbitrary odd order of convergence based on the sine function

논문 개요

이 논문은 사인 (sine) 함수를 기반으로 하여, 임의의 홀수 차수 (arbitrary odd order) 로 수렴하는 고정점 반복법 (fixed-point iteration method) 을 제시하여 원주율 ($\pi$) 을 고정밀도로 계산하는 새로운 알고리즘을 제안합니다. 저자는 이 방법이 초기값이 $\pi$에 충분히 가까울 때, 매우 빠른 수렴 속도로 $\pi$를 계산할 수 있음을 수학적으로 증명하고 Maple 을 이용한 수치 실험을 통해 그 효율성을 입증했습니다.

1. 문제 제기 (Problem)

• 고정밀도 $\pi$ 계산의 필요성: 과학, 공학 및 암호학 등 다양한 분야에서 원주율 $\pi$의 고정밀도 값이 필요하지만, 기존 알고리즘들은 특정 수렴 속도 제한이나 계산 복잡도 문제를 가질 수 있습니다.
• 수렴 속도의 한계: 일반적인 고정점 반복법 (예: $x_{n+1} = x_n + \sin(x_n)$) 은 3 차 (cubic) 수렴을 보이지만, 더 높은 차수의 수렴을 달성하기 위한 체계적인 방법론이 필요했습니다.

2. 방법론 (Methodology)

가. 핵심 아이디어: arcsin 급수의 역이용

논문은 $\arcsin(t)$의 테일러 급수 전개식을 출발점으로 삼습니다.
$$\arcsin(t) = \sum_{k=1}^{\infty} \left( \prod_{\ell=1}^{k-1} \frac{2\ell-1}{2\ell} \right) \frac{t^{2k-1}}{2k-1}$$
이 급수의 앞 $P$항만을 취하고, $t = \sin(x)$로 치환하여 유한 합을 구성합니다. 이를 통해 $\pi$를 근사하는 고정점 함수 $S(x)$를 정의합니다.

나. 고정점 함수 정의

임의의 자연수 $P \in \mathbb{N}$에 대해 다음과 같은 함수를 정의합니다:
$$S(x) = x + \sum_{k=1}^{P} \left( \prod_{\ell=1}^{k-1} \frac{2\ell-1}{2\ell} \right) \frac{\sin(x)^{2k-1}}{2k-1}$$
이때 반복식은 다음과 같습니다:
$$x_{n+1} = S(x_n)$$

다. 수렴성 증명 (Theorem)

저자는 다음 세 가지 사실을 증명합니다:

1. 고정점 성질: $S(\pi) = \pi$가 성립합니다. ($\sin(\pi)=0$이므로 모든 항이 소거됨)
2. 미분계수 성질: $\pi$에서의 $k$계 도함수 $S^{(k)}(\pi)$는 $1 \le k \le 2P$일 때 $0$이며,$k=2P+1$일 때만 0 이 아닌 값을 가집니다.
   $$S^{(k)}(\pi) = \begin{cases} 0 & 1 \le k \le 2P \\ \left( \prod_{\ell=1}^{P} (2\ell-1) \right)^2 & k = 2P+1 \end{cases}$$
3. 수렴 차수: Banach 고정점 정리에 따라 초기값 $x_0$가 $\pi$에 충분히 가까우면, 이 반복식은 정확히 **$(2P+1)$차 (odd order)**로 $\pi$에 수렴합니다.

라. 최적화 전략 (실제 계산 시)

• 계수 사전 계산: 반복식 내의 계수들은 상수이므로 미리 계산하여 벡터로 저장합니다.
• 점진적 정밀도 향상 (Self-correcting): 초기 반복에서는 낮은 정밀도로 계산하고, 수렴 차수 $(2P+1)$에 따라 필요한 유효숫자만 증가시키며 계산합니다. 이는 고정밀도 사인 함수 계산에 소요되는 시간을 크게 줄여줍니다.
• 초기값 설정: 이미 18 자리 정도 정확한 $\pi$ 근사값을 초기값으로 사용하여 수렴을 가속화합니다.

3. 주요 결과 (Results)

가. 수치 실험 1: 3 차 수렴 ($P=1$)

• 함수: $S(x) = x + \sin(x)$
• 초기값: $x_0 = 3$
• 결과: 7 번의 반복으로 900 자리 이상의 $\pi$를 얻었습니다.
• 수렴 확인: 각 반복 단계마다 오차의 지수 (negative exponent) 가 약 3 배씩 증가하여 3 차 수렴을 확인했습니다.

나. 수치 실험 2: 100 만 자리 계산 ($P=4$, 9 차 수렴)

• 설정: $P=4$로 설정하여 **9 차 수렴 (nonic convergence)**을 달성했습니다.
• 초기값: 18 자리 정밀도의 $\pi$ 근사값 사용.
• 최적화: 6 단계의 반복을 수행하며, 각 단계마다 계산 정밀도 (Digits) 를 $18 \times 9^k$에 맞춰 동적으로 증가시켰습니다.
• 결과: 6 번의 반복으로 **약 100 만 자리 (정확히는 $10^{59070}$ 수준의 정밀도)**의 $\pi$ 값을 얻었습니다.
• 성능: AMD Ryzen 5 5600G CPU 와 Maple 2025.2 를 사용하여 효율적으로 계산되었습니다.

4. 주요 기여 (Key Contributions)

1. 임의의 홀수 차수 수렴 달성: 파라미터 $P$를 조절함으로써 원하는 임의의 홀수 차수 ($3, 5, 7, \dots, 2P+1$) 의 수렴 속도를 가진 $\pi$ 계산 알고리즘을 체계적으로 제시했습니다.
2. 수학적 엄밀성: 사인 함수의 급수 전개와 고계 도함수 분석을 통해 수렴 차수가 정확히 $(2P+1)$임을 엄밀하게 증명했습니다.
3. 계산 효율성 최적화: 고정밀도 사인 함수 계산의 비효율성을 해결하기 위해 '점진적 정밀도 향상' 기법을 적용하여, 100 만 자리 이상의 초고정밀도 계산을 매우 적은 반복 횟수 (6 회) 로 달성할 수 있음을 보였습니다.

5. 의의 및 결론 (Significance)

이 논문은 $\pi$ 계산 분야에서 기존에 알려지지 않았던 사인 함수 기반의 고정점 반복법을 제안하여, 매우 높은 수렴 속도를 가진 새로운 계산 패러다임을 제시했습니다. 특히, $P$ 값을 조절하여 계산 자원과 요구되는 정밀도 사이의 균형을 맞출 수 있는 유연성을 제공하며, 초고정밀도 계산 (High-precision computation) 분야에서 효율적인 알고리즘으로 활용될 수 있는 잠재력을 보여줍니다. 이 방법은 수학적 아름다움 (arcsin 급수의 역이용) 과 실용적인 계산 효율성을 동시에 갖춘 사례로 평가됩니다.