Function utilities::lowpass_fft_deconv

Function Documentation

Eigen::VectorXd utilities::lowpass_fft_deconv(const Eigen::VectorXd &input, const Eigen::VectorXd &kernel, const std::string &flag, double stab = 0.01)

stably solve Fredholm Integral of the first kind by fft based deconvolution with a gaussian low pass filter

Return
Parameters
  • input:
  • kernel:
  • stab: Stabilisation parameter as a proportion of the max magnitude of the kernel. Fixes the lowpass total cutoff point. For example, default 0.01 creates a low pass gaussian filter, between 0 and the point where the FFT of the kernel reaches 1% of its max (usually DC) magnitude. Set <= 0 to circumvent low pass filtering.
  • flag: