Feynman integral reduction: balanced reconstruction of sparse rational functions and implementation on supercomputers in a co-design approach

Integration-by-parts (IBP) reduction is one of the essential steps in evaluating Feynman integrals. A modern approach to IBP reduction uses modular arithmetic evaluations at the specific numerical values of parameters with subsequent reconstruction of the analytic rational coefficients. Due to the large number of sample points needed, problems at the frontier of science require an application of supercomputers. In this article, we present a rational function reconstruction method that fully takes advantage of sparsity, combining the balanced reconstruction method and the Zippel method. Additionally, to improve the efficiency of the finite-field IBP reduction runs, at each run several numerical probes are computed simultaneously, which allows to decrease the resource overhead. We describe what performance issues one encounters on the way to an efficient implementation on supercomputers, and how one should co-design the algorithm and the supercomputer infrastructure. We present characteristic examples of IBP-reduction in the case of massless two-loop fourand five-point Feynman diagrams using a development version of FIRE and give illustrative examples mimicking the reduction of coefficients appearing in scattering amplitudes for post-Minkowskian gravitational binary dynamics.