2024/09/27 by Yigit Gunsur Elmacioglu, Elmacioglu, Yigit Gunsur, Rory Conlin +7
Computer Science · Environmental Science · #Computational Physics and Python Applications #FOS: Computer and information sciences #Performance (cs.PF) #Soil Moisture and Remote Sensing
paper · pdf · doi:10.48550/arxiv.2409.19156
openalex publication_date 2024/09/27 · openalex created_date 2025/10/10 · openalex updated_date 2026/07/30
Zernike polynomials serve as an orthogonal basis on the unit disc, and have proven to be effective in optics simulations, astrophysics, and more recently in plasma simulations. Unlike Bessel functions, Zernike polynomials are inherently finite and smooth at the disc center (r=0), ensuring continuous differentiability along the axis. This property makes them particularly suitable for simulations, requiring no additional handling at the origin. We developed ZERNIPAX, an open-source Python package capable of utilizing CPU/GPUs, leveraging Google's JAX package and available on GitHub as well as the Python software repository PyPI. Our implementation of the recursion relation between Jacobi polynomials significantly improves computation time compared to alternative methods by use of parallel computing while still performing more accurately for high-mode numbers.