kuramoto_state()

MAYAFLUX_API std::vector< double > MayaFlux::Kinesis::Discrete::kuramoto_state	(	std::span< const double >	omegas,
		double	coupling,
		std::span< const double >	phases = {}
	)

Builds the array kuramoto reads and writes.

Parameters

omegas	Natural increment of each oscillator in radians per sample; sets the count
coupling	Strength of the pull toward the group
phases	Starting phases; missing entries are 0

Layout: [N, coupling, omegas, phases].

Definition at line 486 of file Kernels.cpp.

{
    const size_t n = omegas.size();
    std::vector<double> state(2 + 2 * n, 0.0);
    state.front() = static_cast<double>(n);
    state.at(1) = coupling;
 
    for (size_t i = 0; i < n; ++i) {
        state.at(2 + i) = omegas[i];
        state.at(2 + n + i) = i < phases.size() ? phases[i] : 0.0;
    }
    return state;
}