GPU Implementation of Zippel Method for Feynman Integral Reconstruction

The Zippel algorithm performs a rational reconstruction of multivariate polynomials and aims specifically at the sparse case. It is applied in different fields of science, lately becoming an important step in Feynman integral reduction in elementary particle physics. In some cases with multiple variables it might become a bottleneck for the whole evaluation so that different optimizations are required. In this paper we describe how we ported the classical Zippel algorithm together with its balanced version for rational functions to graphical processor units and perform its evaluation on several GPUs.