A highly efficient Voigt program for line profile computation

Evaluation of the Voigt function, a convolution of a Lorentzian and a Gaussian profile, is essential in various fields such as spectroscopy, atmospheric science, and astrophysics. Efficient computation of the function is crucial, especially in applications where the function may be called for an enormous number of times. In this paper, we present a highly efficient novel algorithm and its Fortran90 implementation for the practical evaluation of the Voigt function with accuracy in the order of 1 0 − 6 . The algorithm uses improved fits based on Chebyshev subinterval polynomial approximation for functions in two variables. The algorithm significantly outperforms widely-used competitive algorithms in the literature, in terms of computational speed, making it highly suitable for real-time applications and large-scale data processing tasks. The substantial improvement in efficiency positions the present algorithm and computer code as a valuable tool in relevant scientific domains. The algorithm has been adopted and implemented in the Meudon PDR code at Paris Observatory and is recommended for similar applications and simulation packages. • Highly efficient algorithm for evaluating the Voigt function with an accuracy of 1 0 − 6 . • Chebyshev subinterval polynomial approximation enhances computational performance. • Significant speed improvement over existing algorithms, suitable for real-time applications. • Optimized for large-scale data processing. • Adopted by the Meudon PDR code for its proven superior efficiency.