[논문 리뷰] Fast Computation of Partial Fourier Transforms

Lexing Ying, Sergey Fomel|ArXiv.org|2008. 02. 12.

Scientific Research and Discoveries참고 문헌 11인용 수 4

한 줄 요약

이 논문은 1차원 및 2차원에서 부분 푸리에 변환을 계산하기 위한 두 가지 빠른 알고리즘을 제시한다. 다중 척도 분해를 통해 합계 영역를 분할하고, 각 성분에 대해 빠른 변환(1D에서는 분수형 FFT, 2D에서는 버터플라이 알고리즘)을 적용한다. 1D 알고리즘은 정확하며 $O(Nackepsilon^2 N)$의 복잡도를 가지며, 2D 버전은 근사적이지만 정확하며 $O(N^2ackepsilon^2 N)$의 복잡도를 가지며, 자유도 수준에 대해 거의 선형 스케일링을 달성한다.

ABSTRACT

We introduce two efficient algorithms for computing the partial Fourier transforms in one and two dimensions. Our study is motivated by the wave extrapolation procedure in reflection seismology. In both algorithms, the main idea is to decompose the summation domain of into simpler components in a multiscale way. Existing fast algorithms are then applied to each component to obtain optimal complexity. The algorithm in 1D is exact and takes $O(N\log^2 N)$ steps. Our solution in 2D is an approximate but accurate algorithm that takes $O(N^2 \log^2 N)$ steps. In both cases, the complexities are almost linear in terms of the degree of freedom. We provide numerical results on several test examples.

연구 동기 및 목표

• 1D 및 2D에서 합계 제약 조건이 있는 부분 푸리에 변환의 계산 병목 현상을 해결하기 위해, 특히 지구물리학적 이미징에서의 파동 전파에 응용한다.
• 자유도 수에 대해 선형에 가까운 복잡도를 가지는 알고리즘을 개발하여 직접 계산의 $O(N^2)$ 및 $O(N^4)$ 비용을 극복한다.
• 변수이고 매끄러운 절단 영역 $c_x = N\lceil c^0(x/N)\rceil$로 정의된 제약 조건 하에서 푸리에 변환을 효율적으로 계산할 수 있도록 한다.
• 반사 지구물리학에서 전파되는 파동 모드만 필요한 상황에서 파동장 전파에 대한 확장 가능한 솔루션을 제공한다.

제안 방법

• 1D 또는 2D에서 합계 영역 $D = \{(x,k) \mid k < c_x\}$ 또는 $|k| < c_x$를 반복적인 다중 척도 분할을 통해 이진 상자로 분해한다.
• 1D에서는 각 상자에 분수형 푸리에 변환을 적용하여 정확한 $O(N\backepsilon^2 N)$ 복잡도를 달성한다.
• 2D에서는 희소 푸리에 변환에 적합한 버터플라이 알고리즘을 사용하여 비정규적인 합계 영역을 효율적으로 처리한다.
• 각 상자 내의 커널 행렬을 대각형 및 진동 성분으로 분해하여 빠른 행렬-벡터 곱셈을 가능하게 한다.
• 경계 근처의 상자는 세분화하고 내부의 큰 상자들이 주도하도록 계층적 분해를 구성하여 총 작업량을 최소화한다.
• 모든 상자들의 합집합이 경계 곡선에서 멀어질수록 기하급수적으로 증가하는 상자 크기를 갖는 합계 영역 $D$를 정확히 커버하도록 보장한다.

실험 결과

연구 질문

• RQ1 변수 합계 한계를 가진 부분 푸리에 변환은 1D에서 $O(N^2)$보다 더 효율적으로 계산될 수 있는가?
• RQ2 2D에서는 고차원 비정규 영역로 인해 유사한 접근 방식이 비록 복잡하지만 거의 선형 복잡도를 달성할 수 있는가?
• RQ3 합계 영역의 다중 척도 분해가 기존의 빠른 변환 알고리즘 사용을 가능하게 하는가?
• RQ4 2D에서 부분 변환을 근사함에 있어 정확성과 효율성 사이의 상충 관계는 어떠한가?
• RQ5 제안된 방법은 지구물리학적 이미징에서 파동 전파 워크플로우에 통합되어 계산 비용을 줄일 수 있는가?

주요 결과

• 1D 알고리즘은 $O(N\backepsilon^2 N)$ 연산으로 정확한 계산을 달성하며, 표준 FFT보다 뿐만 아니라 $O(\backepsilon N)$ 배 더 비싸다.
• 2D 알고리즘은 근사적이지만 정확하며 $O(N^2\backepsilon^2 N)$ 복잡도를 가지며, 자유도 수에 대해 거의 선형으로 스케일링된다.
• 수치 결과에서는 $N=2048$일 때 직접 평가 대비 100배 이상, 전체 FFT 대비 1000배 이상의 속도 향상을 보였다.
• 실제 지구물리학적 속도 모델을 가진 2D의 경우, 상대 오차가 $10^{-8}$ 이하로 유지되면서도 직접 방법 대비 1000배 이상의 속도 향상을 달성했다.
• 2D 알고리즘의 상수 계수는 버터플라이 절차로 인해 높아 보이며, 구현 측면에서 최적화 여지가 있음을 시사한다.
• 다양한 매끄러운 절단 함수에 대해 안정적이며, $N$ 증가에 따라 잘 스케일링됨을 보였으며, $N$이 2048까지 확장된 경우에도 확인되었다.