[论文解读] 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作为黑箱，加速对插值核函数的求和计算，将整体复杂度降低为节点数的线性函数。
• 该方法利用指数核的渐近光滑性，并应用基于斯特林近似和分式插值公式的误差界。
• 通过一维方法的张量化，将该方法扩展至多变量问题，保持线性复杂度与精度。

实验结果

研究问题

• 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$，表明其具有鲁棒性与精度控制能力。
• 数值实验验证了理论预测，显示在各种参数配置下均具有快速收敛性与稳定性。
• 该方法通过将问题视为具有分层核近似的广义非均匀FFT，实现了复单位圆盘中非均匀节点上多项式的快速求值。