Numerical evaluation of harmonic polylogarithms