phasor_bank()

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

noexcept

Sum of N sinusoidal partials, each with its own slow amplitude swell.

coefs is the array phasor_bank_state builds. Every call advances each partial's phase and swell phase by its increment and returns sum(scale * sin(phase) * (0.5 + 0.5 * sin(swell_phase))). The window is not read. An array too short for its stated count gives zero. Increments must stay below one revolution per sample.

Definition at line 238 of file Kernels.cpp.

{
    if (coefs.empty())
        return 0.0;
 
    const auto n = static_cast<size_t>(std::max(coefs.front(), 0.0));
    if (n == 0 || coefs.size() < 1 + 5 * n)
        return 0.0;
 
    const auto increments = coefs.subspan(1, n);
    const auto scales = coefs.subspan(1 + n, n);
    const auto am_increments = coefs.subspan(1 + 2 * n, n);
    const auto phases = coefs.subspan(1 + 3 * n, n);
    const auto am_phases = coefs.subspan(1 + 4 * n, n);
 
    double sum = 0.0;
    for (size_t i = 0; i < n; ++i) {
        phases[i] += increments[i];
        wrap_revolution(phases[i]);
 
        const double swell = 0.5 + 0.5 * std::sin(am_phases[i]);
 
        am_phases[i] += am_increments[i];
        wrap_revolution(am_phases[i]);
 
        sum += scales[i] * swell * std::sin(phases[i]);
    }
    return sum;
}

References sum().

Here is the call graph for this function: