[논문 리뷰] Fast and highly accurate computation of Chebyshev expansion coefficients of analytic functions

Haiyong Wang, Daan Huybrechs|arXiv (Cornell University)|2014. 04. 09.

Iterative Methods for Nonlinear Equations참고 문헌 19인용 수 4

한 줄 요약

이 논문은 복소평면에서의 경로적분을 사용하여 해석적 함수의 체비셰프 전개 계수를 계산하는 새로운 방법을 제시한다. 이 방법은 기계 정밀도 상대오차를 달성한다. 각 계수에 대해 최적의 경로 반지름을 최적화함으로써, 매우 작은 계수일지라도 높은 정밀도를 유지하며, 기계 정밀도에 가까운 고차 도함수 및 도함수의 근을 정확도를 잃지 않고 계산할 수 있다.

ABSTRACT

Chebyshev expansion coefficients can be computed efficiently by using the FFT, and for smooth functions the resulting approximation is close to optimal, with computations that are numerically stable. Given sufficiently accurate function samples, the Chebyshev expansion coefficients can be computed to machine precision accuracy. However, the accuracy is only with respect to absolute error, and this implies that very small expansion coefficients typically have very large relative error. Upon differentiating a Chebyshev expansion, this relative error in the small coefficients is magnified and accuracy may be lost, especially after repeated differentiation. At first sight, this seems unavoidable. Yet, in this paper, we focus on an alternative computation of Chebyshev expansion coefficients using contour integrals in the complex plane. The main result is that the coefficients can be computed with machine precision relative error, rather than absolute error. This implies that even very small coefficients can be computed with full floating point accuracy, even when they are themselves much smaller than machine precision. As a result, no accuracy is lost after differentiating the expansion, and even the 100th derivative of an analytic function can be computed with near machine precision accuracy using standard floating point arithmetic. In some cases, the contour integrals can be evaluated using the FFT, making the approach both highly accurate and fast.

연구 동기 및 목표

• 표준 FFT 기반 체비셰프 계수 계산에서 작은 계수 및 그 도함수의 정밀도 손실 문제를 해결하기 위해.
• 계수 계산에서 절대오차가 아니라 기계 정밀도 상대오차를 달성하는 방법을 개발하기 위해.
• 표준 단정밀도 부동소수점 산술을 사용하여 해석적 함수의 고차 도함수 계산을 정확히 수행할 수 있도록 하기 위해.
• 스펙트럼 방법의 안정성과 정밀도를 부드러운 함수의 고차 도함수 및 근 찾기 문제로 확장하기 위해.
• 각 계수에 대해 최적의 경로 반지름을 사용하여 체비셰프 계수를 계산하는 프레임워크를 제공함으로써, 강건성과 정밀도를 보장하기 위해.

제안 방법

• 이 방법은 복소평면에서의 경로적분을 사용하여 체비셰프 전개 계수를 계산하며, 수치적 적분에는 사다리꼴 법칙을 적용한다.
• 각 계수 $ a_n $에 대해 상대오차를 최소화하기 위해 최적의 반지름 $ \rho^*(n) $을 선택함으로써, 매우 작은 계수일지라도 높은 정밀도를 확보한다.
• 경로적분 공식은 동일한 반지름을 모든 계수에 사용할 경우 FFT를 활용할 수 있도록 하여 계산 효율성을 유지한다.
• 각 계수에 대해 최적의 반지름을 사용할 경우, 도함수 계산 후 상대오차가 증폭되는 표준 FFT 기반 접근법의 문제를 피할 수 있다.
• 이 방법은 복소해석학과 카우치 적분 공식을 활용하며, 테일러 계수를 위한 보른레만의 접근법과 유사하지만, 체비셰프 전개에 적응시켰다.
• 도함수의 계수는 도함수 전개의 계수에 대한 재귀관계를 사용하여 체비셰프 전개를 미분함으로써 계산된다.

실험 결과

연구 질문

• RQ1 해석적 함수의 체비셰프 전개 계수를 절대오차가 아니라 기계 정밀도 상대오차로 계산할 수 있는가?
• RQ2 각 계수에 대해 최적의 경로 반지름을 사용할 경우, 고차 도함수 계산의 정밀도가 상당히 향상되는가?
• RQ3 표준 부동소수점 산술을 사용하여 100차 도함수의 정확도를 높게 유지할 수 있는가?
• RQ4 최적의 반지름을 사용한 계수 계산을 통해 고차 도함수의 근을 정확히 계산할 수 있는가?
• RQ5 각 계수에 대한 최적의 반지름을 효율적이고 자동으로 결정할 수 있는가?

주요 결과

• 제안된 방법은 기계 에psilon보다 훨씬 작은 계수일지라도 기계 정밀도 상대오차로 체비셰프 계수를 계산한다.
• 고차 도함수, 특히 100차 도함수까지도 표준 부동소수점 산술을 사용하여 기계 정밀도에 가까운 정밀도로 계산할 수 있다.
• 각 계수 $ a_n $에 대해 최적의 반지름 $ \rho^*(n) $을 사용함으로써 상대오차가 작게 유지되어 도함수 계산 중 오차 증폭을 방지한다.
• 수치 실험 결과, 각 계수에 대해 개별 최적 반지름을 사용할 경우, 모든 계수에 동일한 반지름을 사용하는 것보다 도함수의 근 계산에서 훨씬 낮은 오차를 기록한다.
• 모든 계수에 동일한 반지름을 사용할 경우 계산 효율성이 유지되어 FFT 기반 가속이 가능하다.
• 이 방법은 도함수 전개 계수를 정확히 계산함으로써 고차 도함수의 근 찾기(예: 국소 최대점, 변곡점 등)를 정확히 수행할 수 있다.