[논문 리뷰] Recurrence relations and fast algorithms

Mark Tygert|arXiv (Cornell University)|2006. 09. 14.

Numerical Methods and Algorithms참고 문헌 12인용 수 7

한 줄 요약

이 논문은 선형 재귀 관계를 만족하는 특수 함수—예를 들어 가중 치수 수직 다항식과 베셀 함수—의 변환을 빠르게 계산하는 알고리즘을 제시한다. 이러한 재귀 관계로부터 유도된 삼중대각 자기수반 행렬의 스펙트럼 분해를 활용함으로써, 계산 복잡도를 $ O(n \ln n (\ln(1/\varepsilon))^{3}) $로 줄여 기존의 표준적 $ O(n^2) $ 방법보다 훨씬 빠르면서도 높은 정밀도를 유지한다.

ABSTRACT

We construct fast algorithms for evaluating transforms associated with families of functions which satisfy recurrence relations. These include algorithms both for computing the coefficients in linear combinations of the functions, given the values of these linear combinations at certain points, and, vice versa, for evaluating such linear combinations at those points, given the coefficients in the linear combinations; such procedures are also known as analysis and synthesis of series of certain special functions. The algorithms of the present paper are efficient in the sense that their computational costs are proportional to n (ln n) (ln(1/epsilon))^3, where n is the amount of input and output data, and epsilon is the precision of computations. Stated somewhat more precisely, we find a positive real number C such that, for any positive integer n > 10, the algorithms require at most C n (ln n) (ln(1/epsilon))^3 floating-point operations and words of memory to evaluate at n appropriately chosen points any linear combination of n special functions, given the coefficients in the linear combination, where epsilon is the precision of computations.

연구 동기 및 목표

• 재귀 관계를 만족하는 특수 함수의 선형 조합을 평가하기 위한 효율적인 알고리즘을 개발하기 위해.
• 정규직교 다항식과 베셀 함수를 포함한 급수의 빠른 분석 및 합성 가능하게 하기 위해.
• 고정밀도 변환에 대해 계산 비용을 $ O(n^2) $ 에서 $ O(n \ln n (\ln(1/\varepsilon))^3) $ 으로 감소시키기 위해.
• 표준 FFT를 넘어서 보다 넓은 범주에 속하는 특수 함수에까지 빠른 스펙트럼 방법의 적용 가능성을 확장하기 위해.
• 관련 레전드르 함수를 통한 구면 조화함수 및 관련 변환에 대한 실용적 구현을 지원하기 위해.

제안 방법

• 재귀 관계의 해가 삼중대각 자기수반 행렬의 고유벡터와 대응된다는 사실을 활용하기 위해.
• 고유값과 고유벡터 행렬 $ U $ 와 $ U^T $ 를 효율적으로 계산하기 위해 분할 정복 스펙트럼 알고리즘을 사용하기 위해.
• 고유벡터 행렬 $ U $ 와 $ U^T $ 와의 행렬-벡터 곱셈을 빠르게 수행하기 위한 빠른 알고리즘을 적용하기 위해, 이는 순환 이분법과 역반복 기반의 방법을 포함한다.
• 계수 공간과 함수 값 간의 변환을 위해 대각 스케일링 행렬 $ S $ 와 $ S^{-1} $ 를 활용하기 위해.
• 재귀 기반 함수 표현과 구적 노드, 크리스토펠 수를 결합하여 정확한 함수 평가를 가능하게 하기 위해.
• 가중 함수 $ w(x) = (1-x)^m(1+x)^m $ 를 갖는 정규직교 다항식과 다양한 차수의 베셀 함수에 대해 이 방법을 적용하기 위해.

실험 결과

연구 질문

• RQ1 선형 재귀 관계를 만족하는 특수 함수의 변환을 위한 빠른 알고리즘을 구축할 수 있는가?
• RQ2 이러한 알고리즘의 계산 복잡도는 무엇이며, 입력 크기 $ n $ 과 정밀도 $ \varepsilon $ 에 따라 유계로 표현될 수 있는가?
• RQ3 삼중대각 행렬의 스펙트럼 분해를 어떻게 활용하여 정규직교 다항식의 변환을 가속화할 수 있는가?
• RQ4 동일한 프레임워크를 베셀 함수와 구면 조화함수로 확장할 수 있는가?
• RQ5 크리스토펠 수와 구적 노드는 선형 조합의 빠른 평가를 가능하게 하는 데 어떤 역할을 하는가?

주요 결과

• 제안된 알고리즘은 $ n \geq 10 $ 이고 $ \varepsilon \leq 1/10 $ 인 경우, 최대 $ C n \ln n (\ln(1/\varepsilon))^3 $ 개의 부동소수점 연산과 메모리 단위를 요구하며, 여기서 $ C $ 는 전역 상수이다.
• 이 방법을 통해 $ n $ 개의 특수 함수를 $ n $ 개의 점에서 선형 조합으로 평가하는 데 $ O(n \ln n (\ln(1/\varepsilon))^3) $ 번의 연산으로 수행할 수 있다.
• 함수 값에서 계수를 복원하는 역문제 역시 고유벡터 행렬의 수반 행렬을 사용하여 동일한 복잡도로 해결할 수 있다.
• 크리스토펠 수는 고유벡터 행렬의 첫 번째 열을 통해 동일한 프레임워크를 통해 빠르게 계산할 수 있다.
• 베셀 함수의 영점은 스펙트럼 알고리즘의 고유값 계산 단계를 통해 효율적으로 계산할 수 있다.
• 적절한 가중 함수를 갖는 관련 레전드르 함수에 이 프레임워크를 적용함으로써, 구면 조화함수에 대한 빠른 변환을 지원할 수 있다.