We present the FORTRAN-code HPOLY.f for the numerical calculation of harmonic
polylogarithms up to w = 8 at an absolute accuracy of ∼4.9⋅10−15
or better. Using algebraic and argument relations the numerical representation
can be limited to the range x∈[0,2−1]. We provide replacement
files to map all harmonic polylogarithms to a basis and the usual range of
arguments x∈]−∞,+∞[ to the above interval analytically. We also
briefly comment on a numerical implementation of real valued cyclotomic
harmonic polylogarithms.Comment: 19 pages LATEX, 3 Figures, ancillary dat