Reducing round-off errors in rigid body dynamics

In several recent publications, numerical integrators based on Jacobi elliptic functions are proposed for solving the equations of motion of the rigid body. Although this approach yields theoretically the exact solution, a standard implementation shows an unexpected linear propagation of round-off errors. We explain how deterministic error contribution can be avoided, so that round-off behaves like a random walk.