tapped_sum()

MAYAFLUX_API double MayaFlux::Kinesis::Discrete::tapped_sum ( std::span< const double >  window, std::span< const double >  coefs  )

noexcept

Sum of gain * window[lag] over (lag, gain) pairs.

coefs holds pairs: coefs[2j] is a lag in samples, coefs[2j + 1] its gain. Fractional lags interpolate linearly between neighbours. A pair whose lag lies outside the window contributes nothing, so a short window simply reaches fewer taps.

Definition at line 122 of file Kernels.cpp.

{
    if (window.empty())
        return 0.0;
 
    const auto reach = static_cast<double>(window.size() - 1);
    double sum = 0.0;
 
    for (size_t j = 0; j + 1 < coefs.size(); j += 2) {
        const double lag = coefs[j];
        if (!(lag >= 0.0) || lag > reach)
            continue;
 
        const auto index = static_cast<size_t>(lag);
        const double fraction = lag - static_cast<double>(index);
        double value = window[index];
        if (fraction > 0.0 && index + 1 < window.size())
            value += (window[index + 1] - value) * fraction;
 
        sum += coefs[j + 1] * value;
    }
    return sum;
}

References index, sum(), value, and window.

Here is the call graph for this function: