高斯求积:高效数值积分的利器

2025-06-08

本文介绍了一种高效的数值积分技术——高斯求积,特别是切比雪夫-高斯求积法。它通过在特殊节点上评估函数值并进行加权求和来逼近定积分。与传统方法相比,高斯求积法只需较少的节点就能达到更高的精度,尤其适用于[-1,1]区间上的积分。文章还讲解了如何将一般区间和函数形式转化为切比雪夫-高斯求积法适用的形式,并通过实例演示了其应用和优势,最终在海平面变化率估计中得到应用。

阅读更多

JAX 中的快速扫描法求解 Eikonal 方程

2025-05-11

本文介绍了使用 JAX 实现快速扫描法求解 Eikonal 方程,该方程在界面演化和图像处理等领域应用广泛。文章首先解释了水平集方法和 Eikonal 方程的概念,然后详细阐述了快速扫描法的算法原理,包括网格设置、迭代更新以及 Godunov 上风差分格式的使用。最后,作者提供了基于 NumPy 和 JAX 的代码实现,并进行了性能比较,结果显示 JAX 版本显著提升了计算效率。此外,文章还探讨了并行化快速扫描法的尝试及遇到的挑战。

阅读更多