“…Its modular design paves the way for future methodological extensions, including detection of Hopf bifurcations [107], limit cycles and chaos [108], higher-order Krylov-Bogoliubov averaging method [109], as well as interfaces with existing dedicated libraries to treat nonlinear spatially-extended or quantum systems [110,111]. Its usage can assist a breadth of fields, where nonlinear harmonically-driven systems appear, such as modal analysis in structural dynamics [10,33], electric circuits [11,34,35], nonlinear optics [12,[36][37][38][39][40][41][42], optomechanics [22,[43][44][45], micro-and nanomechanics [18,46,[51][52][53][54][55][56][57], oscillator networks [58][59][60][61][62][63], Ising machines [64][65][66][67][68][69], and many-body light-matter systems [70]…”