Gaussian Quadrature: A Powerful Numerical Integration Technique

2025-06-08

This blog post explores Gaussian quadrature, a powerful numerical integration technique, specifically Chebyshev-Gauss quadrature. It approximates definite integrals by evaluating the function at specific nodes and summing the weighted values. Compared to traditional methods, it achieves higher accuracy with fewer nodes, particularly for integrals over the interval [-1,1]. The post explains how to adapt general intervals and function forms to fit the Chebyshev-Gauss quadrature, demonstrating its application and advantages with an example. The technique found application in estimating sea level change rates.

Read more

Solving the Eikonal Equation with Fast Sweeping in JAX

2025-05-11

This blog post explores solving the Eikonal equation, crucial in interface evolution and image processing, using the fast sweeping method implemented in JAX. It begins by explaining level set methods and the Eikonal equation itself. The core of the post details the fast sweeping algorithm, covering grid setup, iterative updates, and the Godunov upwind scheme. NumPy and JAX code implementations are provided, with benchmarks demonstrating JAX's significant speed advantage. The author also discusses attempts at parallelizing the algorithm and the challenges encountered.

Read more