  The program ftint is introduced which numerically evaluates dimensionally
regulated integrals as they occur in the perturbative approach to the
gradient-flow formalism in quantum field theory. It relies on sector
decomposition in order to determine the coefficients of the individual orders
in $\epsilon=(4-D)/2$, where $D$ is the space-time dimension. For that purpose,
it implements an interface to the public library pySecDec. The current version
works for massive and massless integrals up to three-loop level with vanishing
external momenta, but the underlying method is extendable to more general
cases.
