Fast logarithmic Fourier-Laplace transform of nonintegrable functions
arXiv:1812.09575 · doi:10.1103/PhysRevE.100.053302
Abstract
We present an efficient and very flexible numerical fast Fourier-Laplace transform, that extends the logarithmic Fourier transform (LFT) introduced by Haines and Jones [Geophys. J. Int. 92(1):171 (1988)] for functions varying over many scales to nonintegrable functions. In particular, these include cases of the asymptotic form and with arbitrary real . Furthermore, we prove that the numerical transform converges exponentially fast in the number of data points, provided that the function is analytic in a cone with a finite opening angle around the real axis and satisfies as with a positive constant , which is the case for the class of functions with power-law tails. Based on these properties we derive ideal transformation parameters and discuss how the logarithmic Fourier transform can be applied to convolutions. The ability of the logarithmic Fourier transform to perform these operations on multiscale (non-integrable) functions with power-law tails with exponentially small errors makes it the method of choice for many physical applications, which we demonstrate on typical examples. These include benchmarks against known analytical results inaccessible to other numerical methods, as well as physical models near criticality.
arXiv comment14 pages, 8 figures
References in corpus (11)
- Feshbach Resonances in Ultracold Gases
- Cold atoms in cavity-generated dynamical optical potentials
- Energetics of a strongly correlated Fermi gas
- The `Higgs' Amplitude Mode at the Two-Dimensional Superfluid-Mott Insulator Transition
- Mode-Coupling Theory (MCT) Lecture Notes
- Feynman diagrams versus Fermi-gas Feynman emulator
- Precise determination of the structure factor and contact in a unitary Fermi gas
- Temperature dependence of the contact in a unitary Fermi gas
- Contact and Momentum Distribution of the Unitary Fermi Gas
- Universal phase diagram and scaling functions of imbalanced Fermi gases
- Diagrammatic Monte Carlo algorithm for the resonant Fermi gas