Hi, I saw the tutorial about the Keppler problem.
https://github.com/SciML/SciMLTutorials.jl/blob/d5ece1515327f5931618e48c3a06eec42461e5ce/notebook/models/05-kepler_problem.ipynb

I'm tyring to study a Hamiltonian where $N~2500$ rotors interact in a periodic 2d grid,

I've implemented the Verlet method using the following work as reference
A Study of XY Model by Spin-dynamic simulation
Unfortunately, my code is really slow.

I’m trying to understand if it’s possible to attack this kind of problem using your package. It is?

I could'nt find any point in the documentation which could help me.