  Drell-Yan lepton pair production processes are extremely important for
Standard Model (SM) precision tests and for beyond the SM searches at hadron
colliders. Fast and accurate predictions are essential to enable the best use
of the precision measurements of these processes; they are used for parton
density fits, for the extraction of fundamental parameters of the SM, and for
the estimation of background processes in searches. This paper describes a new
numerical program, DYTurbo, for the calculation of the QCD transverse-momentum
resummation of Drell-Yan cross sections up to next-to-next-to-leading
logarithmic accuracy combined with the fixed-order results at
next-to-next-to-leading order ($\mathcal{O}(\alpha_{\mathrm{S}}^2)$), including
the full kinematical dependence of the decaying lepton pair with the
corresponding spin correlations and the finite-width effects. The DYTurbo
program is an improved reimplementation of the DYqT, DYqT and DYNNLO programs,
which provides fast and numerically precise predictions through the
factorisation of the cross section into production and decay variables, and the
usage of quadrature rules based on interpolating functions for the integration
over kinematic variables.
