HyperPrecision: a Mathematica package for high-precision numerical evaluation of multivariate hypergeometric functions
Description
In this paper, we present HyperPrecision, a Mathematica package for high-precision numerical evaluation of general Horn-type multivariate hypergeometric functions and their Laurent expansions in a small parameter ϵ. Such functions appear widely in physics and mathematics, with applications ranging from quantum field theory and string theory to number theory and statistics. Their high-precision numerical evaluation, however, remains challenging, since their defining series converge only in restricted domains and analytic continuation beyond these domains is, in general, non-trivial. HyperPrecision addresses this problem by automatically constructing the Pfaffian system of partial differential equations for a given hypergeometric function and restricting it to a one-dimensional contour in the space of variables connecting the starting point to the target point. The resulting ordinary differential equation is then solved by the Frobenius method, with the boundary conditions analytically determined by the defining series. We illustrate the use of the package by evaluating commonly occurring multivariate hypergeometric functions, including the Appell F1, F2, F3, and F4 functions, the Horn G- and H-series, and the Lauricella FA, FB, FC, and FD functions, as well as by considering applications to angular integrals, Feynman integrals, and cosmological and holographic correlators.