  We present FastGPL, a C++ library for the fast evaluation of generalized
polylogarithms which appear in many multi-loop Feynman integrals. We implement
the iterative algorithm proposed by Vollinga and Weinzierl in a two-step
approach, i.e., we generate concise expressions using an external program and
hard-code them into the numeric library. This allows efficient and accurate
numeric evaluations of generalized polylogarithms suitable for Monte Carlo
integration and event generation. Floating-point arithmetics are carefully
taken care of to avoid loss of accuracy. As an application and demonstration,
we calculate the two-loop corrections for Higgs boson production in the vector
boson fusion channel at electron-positron colliders. FastGPL is expected to be
useful for event generators at the next-to-next-to-leading order accuracy.
