Here we present my numerical code precession, which implements our multi-timescale way to look at spinning black-hole binaries. The paper has a detailed description of the various functions as well as lots of examples.

Update: typos in Eq. (36-37) have been fixed in v3 on the arXiv.

Davide Gerosa, Michael Kesden.
Physical Review D 93 (2016) 124066.
arXiv:1605.01067 [astro-ph.HE].
Open-source code: homepage, repository, documentation.