hypersynchronization.rhs_pairwise_a2a_harmonic#
- hypersynchronization.rhs_pairwise_a2a_harmonic(t, theta, omega, k1, k2)[source]#
Right-hand side of the ODE, all-to-all pairwise coupling with a harmonic.
Coupling function: k1*sin(theta_j - theta_i) + k2*sin(2*theta_j - 2*theta_i). Uses a vectorised implementation that avoids explicit matrix products.
- Parameters:
- tfloat
Time (unused, present for compatibility with integrators).
- thetanumpy.ndarray
Phases of the oscillators.
- omegafloat or numpy.ndarray
Natural frequencies of the oscillators.
- k1float
Coupling strength of the first harmonic.
- k2float
Coupling strength of the second harmonic.
- Returns:
- dthetanumpy.ndarray of shape (N,)
Rate of change of phases.
See also