“…In this paper we concentrate exclusively on harmonic polylogarithms (HPL's) up to weight four, which since their introduction have found many applications in computations up to two-loop order in the perturbative expansion, e.g., [7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29]. In order to confront the theoretical next-to-next-to-leading order (NNLO) predictions to experiment, it is mandatory to be able to evaluate HPL's numerically in a fast and accurate way.…”