We describe a novel N-body code designed for simulations of the central regions of galaxies containing massive black holes. The code incorporates Mikkola’s “algorithmic” chain regularization scheme including post-Newtonian ...