fourier_qwe¶

empymod.transform.fourier_qwe(fEM, time, freq, ftarg)[source]¶

Fourier Transform using Quadrature-With-Extrapolation.

It follows the QWE methodology [Key12] for the Hankel transform, see hankel_qwe() for more information.

The function is called from one of the modelling routines in empymod.model. Consult these modelling routines for a description of the input and output parameters.

This function is based on get_CSEM1D_TD_QWE.m from the source code distributed with [Key12].

fourier_qwe() checks how steep the decay of the frequency-domain result is, and calls QUAD for the very steep interval, for which QWE is not suited.

Returns

tEMarray: Returns time-domain EM response of fEM for given time.
convbool: If true, QWE/QUAD converged. If not, ftarg might have to be adjusted.