deconvolve()

std::vector< double > MayaFlux::Kinesis::Discrete::deconvolve	(	std::span< const double >	src,
		std::span< const double >	ir,
		double	regularization = 1e-6
	)

Frequency-domain deconvolution with Tikhonov regularisation.

Parameters

src	Input signal
ir	Impulse response to invert
regularization	Regularisation factor for numerical stability

Returns

Deconvolved output, same length as src

Definition at line 124 of file Convolution.cpp.

{
    return apply_convolution(src, ir,
        [regularization](const auto& s, const auto& k, auto& r) {
            P::transform(P::par_unseq, s.begin(), s.end(), k.begin(), r.begin(),
                [regularization](const std::complex<double>& sig, const std::complex<double>& ker) {
                    const double mag_sq = std::norm(ker);
                    if (mag_sq < regularization)
                        return std::complex<double> {};
                    return sig * std::conj(ker) / (mag_sq + regularization);
                });
        });
}

References apply_convolution().

Here is the call graph for this function: