  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.
