Improve performance of the LeapFrog tracking

We saw in the profiling by @JAlameddine that tracking takes a significant fraction of total time. It is very likely that this can be significantly optimized. Maybe a higher-order method can be used. Maybe also some analytical propagation could be implemented (helical motion in magnetic field)? Maybe we are tracking with too high precision?