Abstract:The fitting of spectral lines is a common step in the analysis of line observations and simulations. However, the observational noise, the presence of multiple velocity components, and potentially large data sets make it a non-trivial task. We present a new computer program Spectrum Iterative Fitter (SPIF) for the fitting of spectra with Gaussians or with hyperfine line profiles. The aim is to show the computational efficiency of the program and to use it to examine the general accuracy of approximating spectra with simple models. We describe the implementation of the program. To characterise its performance, we examined spectra with isolated Gaussian components or a hyperfine structure, also using synthetic observations from numerical simulations of interstellar clouds. We examined the search for the globally optimal fit and the accuracy to which single-velocity-component and multi-component fits recover true values for parameters such as line areas, velocity dispersion, and optical depth. The program is shown to be fast, with fits of single Gaussian components reaching on graphics processing units speeds approaching one million spectra per second. This also makes it feasible to use Monte Carlo simulations or Markov chain Monte Carlo calculations for the error estimation. However, in the case of hyperfine structure lines, degeneracies affect the parameter estimation and can complicate the derivation of the error estimates. The use of many random initial values makes the fits more robust, both for locating the global $\chi^2$ minimum and for the selection of the optimal number of velocity components.

Fast and Accurate Linear Fitting for an Incompletely Sampled Gaussian Function With a Long Tail [Tips & Tricks]

A Fast, Accurate, and Separable Method for Fitting a Gaussian Function

Curve Fitting and Fairing Using Conic Splines

Fast Gaussian particle filtering algorithm

Fast fitting of spectral lines with Gaussian and hyperfine structure models

Integral Transforms from Finite Data: An Application of Gaussian Process Regression to Fourier Analysis

A Novel Probability Weighting Method To Fit Gaussian Functions

Gauss-Seidel Progressive and Iterative Approximation for Least Squares Fitting

IMPROVED RELIABILITY DATA CURVE FITTING METHOD BY CONSIDERING SAMPLES DISTINCTION

Feasibility study on the least square method for fitting non-Gaussian noise data

Fast Adaptive Fourier Integration for Spectral Densities of Gaussian Processes

Fast simulation of truncated Gaussian distributions

A Gaussian approximation filter for nonlinear systems with compound noises

Theoretical And Experimental Analyses Of Gaussian Filtering For Ultraprecision Circular Profile

A novel nonlinear optimization method for fitting a noisy Gaussian activation function

Gaussian Approximation For Non-stationary Time Series with Optimal Rate and Explicit Construction

Sub-Nyquist sensing of Gaussian pulse streams with unknown shape factor based on information fitting

Improved reliability data curve fitting method by considering samples distinction | udoskonalona metoda dopasowywania krzywych do danych niezawodnosciowych uwzgledniajaca róznice miçdzy probkami

The Linear Template Fit

Research on Interpolation and Data Fitting: Basis and Applications