A two-domain MATLAB implementation for efficient computation of the Voigt/complex error function

TL;DR

This paper introduces a two-domain MATLAB algorithm for fast and accurate computation of the Voigt/complex error function, optimized for radiative transfer applications and variable input parameters.

Contribution

It presents a novel two-domain scheme that adapts the interpolation grid based on input parameters, improving efficiency over existing methods.

Findings

Meets accuracy requirements for HITRAN-based radiative transfer

Provides rapid computation at small x ranges

Efficient MATLAB implementation for practical use

Abstract

In this work we develop a new algorithm for the efficient computation of the Voigt/complex error function. In particular, in this approach we propose a two-domain scheme where the number of the interpolation grid-points is dependent on the input parameter $y$. The error analysis we performed shows that the MATLAB implementation meets the requirements for radiative transfer applications involving the HITRAN molecular spectroscopic database. The run-time test shows that this MATLAB implementation provides rapid computation, especially at smaller ranges of the parameter $x$.