nchoksi4456

Thanks Hanno. I think my question is a little different though: I'm interested in the runtime of REBOUND's hybrid integrator MERCURIUS (WHFast + IAS15) vs. using just IAS15 alone (both...

Okay, thanks Hanno!

Thanks Hanno. Before sharing any working example scripts with you, I am trying to turn off safe_mode as you suggested. But I occasionally have ejections from the system so I...

Thanks. So my setup looks like this: ``` sim = rebound.Simulation() sim.integrator = "mercurius"; sim.G = cgsG; sim.collision = 'direct'; sim.collision_resolve = 'merge' sim.collision_resolve_keep_sorted = 1; sim.track_energy_offset = 1 #...

Hi all, Sorry to revive this old thread. I had a couple other questions about MERCURIUS. In a merger between two bodies m1 and m2, how does REBOUND/MERCURIUS assign dcrit...

Hi Hanno, thanks for the quick response! Just to follow-up, is recalculating dcrit supposed to be computationally expensive? I had experimented with just recalculating dcrit EVERY timestep, regardless of whether...

Gotcha, thanks. Mergers is one concern, which your suggestion to edit the collision routine addresses. My other concern was just scattering: if a planet's semi-major axis changes substantially, then its...

Thanks Hanno and David, this has been super informative! Right now I'm thinking my best bet might be some combination of IAS15 at early times, when the system is most...

Hi Hanno, I made the following update to the `reb_collision_resolve_merge` function in REBOUND to only update dcrit when there is a merger. Would you mind confirming that I'm not accidentally...

Thanks Hanno and Robert for your suggestions! The simulations are short, just a few hundred orbits, so I used the heartbeat method In Python with a fixed WHFast timestep of...