[論文レビュー] Fast evaluation of real and complex exponential sums

Stefan Kunis, Ines Melzer|arXiv (Cornell University)|Jun 10, 2016

Electromagnetic Scattering and Analysis参考文献 20被引用数 5

ひとこと要約

本稿では、階層行列技術と一般化された高速フーリエ変換（FFT）を組み合わせることで、複素単位円板内の非一様ノードにおける多項式の評価を効率的に行う高速なアルゴリズムを提示する。この手法により、振動的および滑らかなカーネル変換に対して、保証された精度を伴いながら線形計算量が達成され、従来の手法に比べ顕著に改善される。

ABSTRACT

Recently, the butterfly approximation scheme and hierarchical approximations have been proposed for the efficient computation of integral transforms with oscillatory and with asymptotically smooth kernels. Combining both approaches, we propose a certain fast Fourier-Laplace transform, which in particular allows for the fast evaluation of polynomials at nodes in the complex unit disk. All theoretical results are illustrated by numerical experiments.

研究の動機と目的

• 振動的および滑らかなカーネルを有する指数和を高速かつ高精度に評価するアルゴリズムの開発。
• 複素単位円板内の非一様ノードにおける多項式評価の計算ボトルネックの解消。
• 階層行列近似と一般化FFTを組み合わせることで、効率性と精度の向上。
• 指数カーネルに対する既存の離散ラプラス変換手法の一般化と改善。
• 理論的誤差境界を備えた自己完結的で数値的に安定な、高速積分変換のフレームワークの提供。

提案手法

• ドメインの二分木的分割を用い、原点からの距離に基づいて積分または和分のドメインを許容可能な区間に分割する。
• 各部分区間内で、チェビシェフノード上でのバーチカルラグランジュ補間を用いてカーネル $\kappa(y,\xi) = e^{-y\xi}$ を近似する階層行列技術を適用する。
• $y\xi$ が小さい領域ではカーネルを 1 で近似し、 $y\xi$ が大きい領域では 0 で近似し、中間領域では $q$ 点のチェビシェフ補間を用いる。
• 一般化FFTをブラックボックスとして用い、補間されたカーネルの和の評価を高速化し、全体の計算量をノード数に線形に抑える。
• 指数カーネルの漸近的滑らかさを活用し、スターリングの近似およびバーチカル補間式から導かれる誤差境界を適用する。
• 1次元手法のテンソル化により多次元問題へと拡張し、線形複雑性と精度を維持する。

実験結果

リサーチクエスチョン

• RQ1 離散ラプラス変換は、線形複雑性と高精度を保ちつつ、より一般化されたカーネル関数を扱えるように一般化可能か？
• RQ2 階層行列技術と一般化FFTを効果的に組み合わせることで、指数和の評価をどのように高速化できるか？
• RQ3 二分木的分割ドメイン上で指数カーネル $e^{-y\xi}$ を補間する際の厳密な誤差境界はどのように導出可能か？
• RQ4 本手法は非一様ノードにおける複素単位円板上での多項式評価を効率的に可能にするか？その計算コストはいかほどか？
• RQ5 振動的および滑らかな積分カーネルの近似において、精度と計算コストのトレードオフはどのようなものか？

主な発見

• 提案手法は、問題サイズおよび目標精度の対数因子を除き、ノード数に対して線形にスケーリングする計算複雑性を達成する。
• カーネル近似の誤差境界は $\left\|\kappa - \mathcal{I}_{q}^{A\times B}\kappa\right\|_{C(A\times B)} \leq \frac{C\mu^{q}q^{\nu}}{2^{2q-1}}\left(2 + \frac{2}{\pi}\log q\right)(\operatorname{dist}(A,0)\operatorname{dist}(B,0))^{-s}$ であり、$C=1/\sqrt{2\pi}$、$\mu=1$、$s=0$、$-\nu=1/2$ で、指数カーネルに対して成り立つ。
• 指数カーネルに対しては補間誤差が $2^{1-2q}$ で抑えられ、$q = \lceil \frac{1}{2} + \log_4(1/\varepsilon) \rceil$ と設定することで、所望の精度 $\varepsilon$ が達成される。
• 指数和の評価における全体の誤差は $\|\mathbf{f} - \mathbf{\tilde{f}}\|_{\infty} \leq \varepsilon \|\mathbf{\hat{f}}\|_1$ を満たし、ロバスト性と精度制御が確認される。
• 数値実験により理論的予測が確認され、さまざまなパrameter領域で高速収束性と安定性が示された。
• 本手法により、階層的カーネル近似を用いた一般化された非等間隔FFTとしての問題に取り組むことで、複素単位円板内の非一様ノードにおける多項式評価が高速に可能となった。