Analytic treatment of complete and incomplete geodesics in Taub-NUT space-times

We present the complete set of analytical solutions of the geodesic equation in Taub-NUT space-times in terms of the Weierstrass elliptic function. We systematically study the underlying polynomials and characterize the motion of test particles by its zeros. Since the presence of the "Misner string" in the Taub-NUT metric has led to different interpretations, we consider these in terms of the geodesics of the space-time. In particular, we address the geodesic incompleteness at the horizons discussed by Misner and Taub, and the analytic extension of Miller, Kruskal and Godfrey, and compare with the Reissner-Nordstr\"om space-time.