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.
Open-source code: homepage, repository, documentation.