[論文レビュー] Computing Highly Oscillatory Integrals

Yunyun Ma, Yuesheng Xu|arXiv (Cornell University)|Jul 2, 2015

Mathematical functions and polynomials被引用数 3

ひとこと要約

本稿では、代数的特異性および停留点を有する高周波数積分に対して、2種類の合成モーメントフリー quadrature ルールを提案する。周波数と特異性の挙動に基づいて積分区間を分割し、弱い特異性部分にはグレーデッド quadrature を適用し、振動部分には修正されたモーメントフリー Filon 形式の手法を適用することで、振動子の逆行列を計算せずに多項式的または指数的収束を達成し、既存手法よりも精度と効率に優れている。

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.

研究の動機と目的

• 高周波数において従来の手法が計算コストが高くなる状況において、特異性および停留点を有する高周波数積分を効率的に計算する課題に対処すること。
• 既存の Filon 形式および Levin 形式手法におけるボトル neck であるモーメントや振動子の逆行列の計算を回避する quadrature ルールの開発。
• 被積分関数の局所的挙動に応じて quadrature ノードの数を調整することにより、多項式的および指数的収束の高い収束次数を達成すること。
• CFCC や FCC 公式などの最先端手法と比較して、計算効率および近似精度を向上させること。
• 量子化学、トモグラフィー、流体動力学における振動的積分を含む実世界の問題に適用可能な堅牢な数値フレームワークの提供。

提案手法

• 波数および特異点・停留点の位置に基づいて積分区間を部分に分割し、各部分区間が弱い特異性を伴いながら急激な振動を有しない、または振動を伴いながら特異性を持たないように保証する。
• 弱い特異性を有する部分区間に対しては、被積分関数の滑らかでない挙動に対処するため、グレーデッド点を用いた古典的 quadrature ルールを適用する。
• 急激な振動を伴うが特異性のない部分区間に対しては、振動子のモーメントの計算を回避する修正されたモーメントフリー Filon 形式の手法を用いる。
• 固定または可変のノード数を用いた区分的多項式補間を用いて、多項式的または指数的収束率を達成する合成 quadrature ルールを構築する。
• FCC や CFCC のような手法で必要とされる非線形振動子 $ g $ の逆行列の計算を回避することで、計算コストを低減する。
• 理論的収束解析は、さまざまな波数および被積分関数の種類において、提案手法と CMFE、FCC、CFCC 公式の比較を含む数値実験によって裏付けられる。

実験結果

リサーチクエスチョン

• RQ1 特異性および停留点を有する高周波数積分に対して、多項式的または指数的収束率を達成する合成 quadrature ルールを設計できるか？
• RQ2 モーメントフリー Filon 形式手法において、振動子の逆行列を回避しつつも高い精度を維持する方法は何か？
• RQ3 局所的な振動および特異性の挙動に基づいて、精度と計算コストのバランスを最適化する積分区間の分割戦略は何か？
• RQ4 提案された合成モーメントフリー quadrature ルールは、CFCC や FCC と比較して精度および効率に優れているか？
• RQ5 非線形ケースにおいて、モーメントの計算や振動子の逆行列の計算を伴わずして指数的収束を達成できるか？

主な発見

• 提案された合成モーメントフリー quadrature ルールは、固定ノード数または可変ノード数の使用に応じて、多項式的または指数的収束率を達成する。
• 数値結果から、CMFE 公式は、特に指数 $ \mu < 0 $ の特異性および $ r > 0 $ の次数の停留点を有する被積分関数において、CFCC 公式よりも高い近似精度を達成することが示された。
• CMFE 法は、FCC や CFCC 公式と比較して、CPU 時間を著しく削減している。例えば、$ \kappa = 10^2 $ で $ f_1(x) = 1 $ の場合、3.35e-2 秒 vs. 2.84e-1 秒であり、関数評価回数は多いにもかかわらず。
• $ f_2(x) = \ln x $ の場合、CMFE 公式は 88 回の関数評価と 4.10e-2 秒を要したが、CFCC 公式は 73 回の評価と 3.90e-1 秒を要した。これは、優れた効率性を示している。
• CMFE 法は $ g^{-1} $ の計算を必要とせず、特に非線形振動子において、FCC や CFCC 手法と比較して顕著な計算時間の短縮が達成されている。
• すべてのテストケースにおいて、CMFE 法は $ \kappa $ が $ 10^2 $ から $ 10^7 $ の広い範囲にわたり、一貫した収束挙動を示しながら高い精度を維持した。