[논문 리뷰] Computing Highly Oscillatory Integrals

Yunyun Ma, Yuesheng Xu|arXiv (Cornell University)|2015. 07. 02.

Mathematical functions and polynomials인용 수 3

한 줄 요약

이 논문은 대수적 특이성과 정적 점을 가진 고주파 적분에 대해 복합적인 모멘트 없는 구적법 두 종류를 제안한다. 진동 및 특이성의 행동에 따라 적분 구간을 분할하고, 약한 특이성 부분에는 격자점 기반의 고전적 구적법을, 진동 부분에는 수정된 모멘트 없는 필론 유형 방법을 적용하여, 진동자 역행렬을 계산하지 않으면서도 다항식 또는 지수 수렴을 달성한다. 이는 기존 방법보다 정확도와 효율성 면에서 뛰어나다.

ABSTRACT

We develop two classes of composite moment-free numerical quadratures for computing highly oscillatory integrals having integrable singularities and stationary points. The first class of the quadrature rules has a polynomial order of convergence and the second class has an exponential order of convergence. We first modify the moment-free Filon-type method for the oscillatory integrals without a singularity or a stationary point to accelerate their convergence. We then consider the oscillatory integrals without a singularity or a stationary point and then those with singularities and stationary points. The composite quadrature rules are developed based on partitioning the integration domain according to the wave number and the singularity of the integrand. The integral defined on a subinterval has either a weak singularity without rapid oscillation or oscillation without a singularity. The classical quadrature rules for weakly singular integrals using graded points are employed for the singular integral without rapid oscillation and the modified moment-free Filon-type method is used for the oscillatory integrals without a singularity. Unlike the existing methods, the proposed methods do not have to compute the inverse of the oscillator. Numerical experiments are presented to demonstrate the approximation accuracy and the computational efficiency of the proposed methods. Numerical results show that the proposed methods outperform methods published most recently.

연구 동기 및 목표

• 고주파 수가 높아질수록 기존 방법이 계산 비용이 급격히 증가하는 상황에서 특이성과 정적 점을 가진 고주파 적분을 효율적으로 계산하는 데 도전한다.
• 기존의 필론 유형 및 레빈 유형 방법에서의 계산 병목 현상이 되는 모멘트 또는 진동자의 역행렬 계산을 피하는 구적법을 개발한다.
• 적분자의 국소적 행동에 따라 구적 노드의 수를 적응적으로 조정하여 다항식 수렴 또는 지수 수렴의 더 높은 수렴 차수를 달성한다.
• CFCC 및 FCC 공식과 같은 최신 기법들에 비해 계산 효율성과 근사 정확도 면에서 향상된 성능을 확보한다.
• 양자 화학, 단층 촬영, 유체역학 등 고주파 적분을 포함하는 실제 문제에 적용 가능한 견고한 수치적 프레임워크를 제공한다.

제안 방법

• 진동수와 특이점, 정적 점의 위치에 따라 적분 구간을 부분구간으로 나누어, 각 부분구간이 약한 특이성이 있지만 빠른 진동이 없거나, 진동은 있지만 특이성이 없는 조건을 만족하도록 한다.
• 약한 특이성이 있는 부분구간에는 적분자의 비연속적 행동을 다루기 위해 격자점 기반의 고전적 구적법을 적용한다.
• 빠른 진동이 있지만 특이성이 없는 부분구간에는 모멘트를 계산할 필요가 없는 수정된 모멘트 없는 필론 유형 방법을 사용한다.
• 고정 또는 가변 수의 노드를 사용하는 조각다항식 보간을 통해 다항식 또는 지수 수렴 속도를 달성하는 복합 구적법을 구성한다.
• FCC 및 CFCC 공식에서 요구되는 비선형 진동자 $ g $ 의 역행렬 계산을 피함으로써 계산 비용을 감소시킨다.
• 이론적 수렴 분석은 다양한 진동수와 적분자 유형에서 제안된 방법을 CMFE, FCC, CFCC 공식과 비교한 수치 실험을 통해 뒷받 Circa.

실험 결과

연구 질문

• RQ1 특이성과 정적 점을 동시에 가진 고주파 적분에 대해 다항식 또는 지수 수렴 속도를 달성할 수 있는 복합 구적법을 설계할 수 있는가?
• RQ2 모멘트 없는 필론 유형 방법에서 진동자의 역행렬을 계산하지 않으면서도 높은 정확도를 유지할 수 있는가?
• RQ3 국소적 진동 및 특이성 행동에 기반하여 정확도와 계산 비용을 균형 잡는 최적의 적분 구간 분할 전략은 무엇인가?
• RQ4 제안된 복합 모멘트 없는 구적법이 CFCC 및 FCC와 같은 기존 방법에 비해 정확도와 효율성 면에서 어떻게 비교되는가?
• RQ5 비선형 경우에서 모멘트를 계산하거나 진동자의 역행렬을 계산하지 않으면서도 지수 수렴을 달성할 수 있는가?

주요 결과

• 제안된 복합 모멘트 없는 구적법은 고정 또는 가변 수의 보간 노드를 사용함에 따라 다항식 또는 지수 수렴 속도를 달성한다.
• 수치 결과는 CMFE 공식이 CFCC 공식보다 더 높은 근사 정확도를 보이며, 특히 지수 $ \mu < 0 $ 인 특이성과 차수 $ r > 0 $ 인 정적 점을 가진 적분자에서 뚜렷한 우수성을 보인다.
• CMFE 방법은 FCC 및 CFCC 공식보다 CPU 시간을 크게 절약한다. 예를 들어 $ \kappa = 10^2 $ 인 경우 $ f_1(x) = 1 $ 일 때 3.35e-2 초 대비 2.84e-1 초를 소비함에도 불구하고 더 많은 기능 평가를 수행한다.
• $ f_2(x) = \ln x $ 일 때, CMFE 공식은 88회의 기능 평가와 4.10e-2 초를 소비하지만, CFCC 공식은 73회의 평가와 3.90e-1 초를 소비하여 더 뛰어난 효율성을 보인다.
• CMFE 방법은 FCC 및 CFCC 방법과 달리 $ g^{-1} $ 를 계산할 필요가 없어, 비선형 진동자에 대해 특히 높은 계산 시간 절감 효과를 보인다.
• 모든 시험 케이스에서 CMFE 방법은 $ \kappa $ 가 $ 10^2 $ 에서 $ 10^7 $ 에 이르는 광범위한 진동수 범위에서 일관된 수렴 행동을 유지하며 높은 정확도를 유지한다.