I get different simulations when I use the same seed, but different number of MPI nodes.
Is this a bug, or a well known feature?
I found this while investigating issue #25 on quasars with Dec=NaN, since I got different number of bad Dec when using MPI and not using it.