CHAPLIN - Complex Harmonic Polylogarithms in Fortran
  We present a new Fortran library to evaluate all harmonic polylogarithms up to weight four numerically for any complex argument. The algorithm is based on a reduction of harmonic polylogarithms up to weight four to a minimal set of basis functions that are computed numerically using series expansions allowing for fast and reliable numerical results. 
Stephan Buehler, Claude Duhr
hep-ph