Convert standard leapfrog to Boris method (add support for electric fields)
The current tracking algorithm, a basic leapfrog, can sometimes exceed c
and isn't a proper relativistic tracker. I think the next step is to move to the Boris integrator, described here, which is relativistic, symplectic, phase-space-preserving, and supports both electric and magnetic fields.