  Accurate spin tracking is a valuable tool for understanding spin dynamics in
particle accelerators and can help improve the performance of an accelerator.
In this paper, we present a detailed discussion of the integrators in the spin
tracking code gpuSpinTrack. We have implemented orbital integrators based on
drift-kick, bend-kick, and matrix-kick splits. On top of the orbital
integrators, we have implemented various integrators for the spin motion. These
integrators use quaternions and Romberg quadratures to accelerate both the
computation and the convergence of spin rotations. We evaluate their
performance and accuracy in quantitative detail for individual elements as well
as for the entire RHIC lattice. We exploit the inherently data-parallel nature
of spin tracking to accelerate our algorithms on graphics processing units.
