hannorein / rebound

💫 An open-source multi-purpose N-body code.
https://rebound.readthedocs.io/
GNU General Public License v3.0
820 stars 217 forks source link

selfgravity_disc_mpi Issue #710

Closed ruaSulaiman closed 9 months ago

ruaSulaiman commented 11 months ago

Hello,

I have noticed a weird thing while testing the code in the examples provided in rebound (c code), specifically the selfgravity_disc_mpi example. The code is using 4 MPI nodes. I tried to change the number of nodes involved to 1, then 2, then went back to 4. The final number of particles left in the simulation is different after integrating the system for the same amount of time (around 1000 yrs). Does this make sense?

hannorein commented 11 months ago

The example uses random initial conditions, so you will alway get a slightly different answer.