I am not very clear on how your algorithms work, but you will need timesteps much less than a gyroperiod so that the full orbits are resolved. This is true for any particle-based code unless there is some trickery I am not aware of.
I just invented said trickery!
the trickery was a bit trickier than i thought, and i had to convert to double-precision, but it works great!
