Hanno Rein

Note that you don't need to run it for very long or need to do any plots. Just compare any coordinate of any particle at any given time. The coordinates...

Is the semi-major axis so large that the timestep is larger than 100yr?

Keep us posted! I'm curious to find out what the issue is. Just to make sure: are you using the latest version of REBOUND? Did you make any modifications?

Hm. Strange. Just to make sure: Are the timesteps always shorter than the output interval? If not, you might end up giving `sim.integrate()` a time in the past (and it...

Maybe this is obvious: The if statement doesn't solve the issue of integrating backwards. It just reports it. Are you saying in the new tests, the if statement never triggers...

Oops. I missed that. Sorry!

I still have a feeling this is somehow related to the issue. What happens if you run it with max_time = 6e4 yr?

Weird indeed. I'll put my thinking cap on.

I've integrated two simulations for some ~7200 yrs, one with max_time=3e4*yr_cgs and 1.5e4*yr_cgs. Both give exactly the same answer down to the last digit. According to your plot, one should...

It can sometimes be tricky to tell which version of rebound and reboundx python is using. Especially if you've compiled reboundX from scratch. The shared library files get installed in...