gala
gala copied to clipboard
Dumping some orbital data in long integrations?
I am finding long integrations (with the integrate_orbit method of the Hamiltonian class) end up eating a lot of memory. Is there a way to dump orbital data from early timesteps to a file, in order to free up some ram? Thank you for your help. (And apologies if this is not the appropriate place to ask this question).
That's a reasonable ask! What kind of API would you find useful? If there was a way to dump information at a regular interval (e.g., every 1000 steps or something) would that be sufficient?
I was thinking something along the following lines
orbit = H.integrate_orbit(w0, dt=step_size, n_steps=num_steps, dump_interval=dump_interval, Integrator=gi.DOPRI853Integrator)
where dump_interval would specify the number of steps between dumps to file. Saving to file at a regular interval is perfectly fine.