A single-domain implementation of the Voigt/complex error function by vectorized interpolation

Abstract

In this work we show how to perform a rapid computation of the Voigt/complex error over a single domain by vectorized interpolation. This approach enables us to cover the entire set of the parameters $x,y\in \mathbb{R}$ required for the HITRAN-based spectroscopic applications. The computational test reveals that within domains $x\in \left[0,15\right]\cap y\in \left[10^{-8},15\right]$ and $x\in \left[0,50000\right]\cap y\ge 10^{-8}$ our algorithmic implementation is faster in computation by factors of about $8$ and $3$, respectively, as compared to the fastest known C/C++ code for the Voigt/complex error function. A rapid MATLAB code is presented.

arXiv comment22 pages

References in corpus (5)

• HELIOS-K: An Ultrafast, Open-source Opacity Calculator for Radiative Transfer
• A rational approximation of the Dawson's integral for efficient computation of the complex error function
• Sampling by incomplete cosine expansion of the sinc function: application to the Voigt/complex error function
• The Voigt and complex error function: Humlíček's rational approximation generalized
• A sampling-based approximation of the complex error function and its implementation without poles

Cited by in corpus (2)

• A two-domain MATLAB implementation for efficient computation of the Voigt/complex error function
• Effect of the instrument slit function on upwelling radiance from a wavelength dependent surface reflectance