[论文解读] Optimal Barycentric Gegenbauer Quadrature

Kareem T. Elgindy|arXiv (Cornell University)|Mar 3, 2016

Matrix Theory and Algorithms参考文献 29被引用 1

一句话总结

本文提出了一种新颖、稳定且高效的Gegenbauer伪谱方法，基于最优重心Gegenbauer积分公式，利用稳定的重心表示法和Gegenbauer-Gauss点的显式重心权重。该方法实现了指数收敛，同时计算成本显著降低，生成的代数系统条件良好，可由标准求解器求解，在效率上优于以往工作，同时保持了高精度。

ABSTRACT

The work reported in this article presents a high-order, stable, and efficient Gegenbauer pseudospectral method to solve numerically a wide variety of mathematical models. The proposed numerical scheme exploits the stability and the wellconditioning of the numerical integration operators to produce well-conditioned systems of algebraic equations, which can be solved easily using standard algebraic system solvers. The core of the work lies in the derivation of novel and stable optimal Gegenbauer quadratures based on the stable barycentric representation of Lagrange interpolating polynomials and the explicit barycentric weights for the Gegenbauer-Gauss (GG) points. A rigorous error and convergence analysis of the proposed quadratures is presented along with a detailed set of pseudocodes for the established computational algorithms. The proposed numerical scheme leads to a reduction in the computational cost and time complexity required for computing the numerical quadrature while sharing the same exponential order of accuracy achieved by [Elgindy and Smith-Miles (2013b)]. The bulk of the work includes three numerical test examples to assess the efficiency and accuracy of the numerical scheme. The present method provides a strong addition to the arsenal of numerical pseudospectral methods, and can be extended to solve a wide range of problems arising in numerous applications.

研究动机与目标

• 开发一种用于求解各类数学模型的高阶、稳定且高效的数值方法，采用伪谱技术。
• 通过引入条件良好的积分算子，解决传统谱方法中常见的不稳定性和病态问题。
• 基于拉格朗日多项式的重心表示法，推导出新型、稳定的最优Gegenbauer积分公式。
• 在保持现有方法指数精度的同时，降低计算成本和时间复杂度。
• 为科学计算和工程应用中的广泛问题提供一个稳健且可扩展的求解框架。

提出的方法

• 采用拉格朗日插值多项式的重心表示法，以增强数值稳定性。
• 为Gegenbauer-Gauss积分点推导出显式的重心权重，以确保稳定性和高效性。
• 设计数值积分算子，使其生成的代数系统条件良好，适合标准求解器求解。
• 开展严格的误差分析与收敛性分析，以验证积分公式的理论性能。
• 提供所有关键计算算法的伪代码，以确保可复现性及实际实现。
• 该方法基于Gegenbauer多项式理论，并利用其谱性质实现高阶精度。

实验结果

研究问题

• RQ1 如何利用重心表示法稳定高阶谱方法中的Gegenbauer积分？
• RQ2 在Gegenbauer-Gauss点中，何种积分权重选择可同时确保稳定性和高阶精度？
• RQ3 是否可在不牺牲指数收敛性的前提下降低Gegenbauer积分的计算成本？
• RQ4 所生成代数系统的条件性与现有谱方法相比如何？
• RQ5 所提出方法在求解基准问题时的实际性能如何？

主要发现

• 所提方法实现了指数收敛，其精度与Elgindy和Smith-Miles（2013b）的前期工作相当。
• 采用稳定重心权重和表示法，生成了条件良好的代数系统，可被标准求解器轻松求解。
• 与传统方法相比，计算成本和时间复杂度显著降低，效率得到提升。
• 在全部三个数值测试示例中，该方法均保持了高精度和稳定性，展现出良好的鲁棒性。
• 提供的伪代码确保了可复现性，并促进了在科学计算工作流中的实现。
• 该框架可扩展应用于应用数学和工程领域的广泛问题。