Bond-order potentials (BOPs) provide a local and physically transparent description of the interatomic interaction. Here we describe the efficient implementation of analytic BOPs in the BOPfox program and library. We discuss the integration of the underlying non-magnetic, collinear-magnetic and noncollinear-magnetic tight-binding models that are evaluated by the analytic BOPs. We summarize the flow of an analytic BOP calculation including the determination of self-returning paths for computing the moments, the self-consistency cycle, the estimation of the band-width from the recursion coefficients, and the termination of the BOP expansion. We discuss the implementation of the calculations of forces, stresses and magnetic torques with analytic BOPs. We show the scaling of analytic BOP calculations with the number of atoms and moments, present options for speeding up the calculations and outline different concepts of parallelisation. In the appendix we compile the implemented equations of the analytic BOP methodology and comments on the implementation. This description should be relevant for other implementations and further developments of analytic bond-order potentials.arXiv:1803.07491v2 [physics.comp-ph]