Related papers: FIESTA 2: parallelizeable multiloop numerical calc…
This article is the second of a series of three presenting an alternative method to compute the one-loop scalar integrals. It extends the results of the first article to general complex masses. Let us remind the main features enjoyed by…
There are three main types of numerical computations for the Bessel function of the second kind: series expansion, continued fraction, and asymptotic expansion. In addition, they are combined in the appropriate domain for each. However,…
With an integrated software package {\tt GRACE}, it is possible to generate Feynman diagrams, calculate the total cross section and generate physics events automatically. We outline the hybrid method of parallel computation of the…
In this paper, we give a detailed account of the algorithm outlined in [1] for Feynman integral reduction and $\varepsilon$-factorised differential equations. The algorithm consists of two steps. In the first step, we use a new geometric…
We present a new method for numerically computing generic multi-loop Feynman integrals. The method relies on an iterative application of Feynman's trick for combining two propagators. Each application of Feynman's trick introduces a…
We present pySecDec, a new version of the program SecDec, which performs the factorisation of dimensionally regulated poles in parametric integrals, and the subsequent numerical evaluation of the finite coefficients. The algebraic part of…
In this paper the C++ version of FIRE is presented - a powerful program performing Feynman integral reduction to master integrals. All previous versions used only Wolfram Mathematica, the current version mostly uses Wolfram Mathematica as a…
We calculate 3-loop master integrals for heavy quark correlators and the 3-loop QCD corrections to the $\rho$-parameter. They obey non-factorizing differential equations of second order with more than three singularities, which cannot be…
An improved method is presented for the numerical evaluation of multi-loop integrals in dimensional regularization. The technique is based on Mellin-Barnes representations, which have been used earlier to develop algorithms for the…
We report on the calculation of multi-loop Feynman integrals for single-scale problems by means of difference equations in Mellin space. The solution to these difference equations in terms of harmonic sums can be constructed algorithmically…
We present the main improvements and new features in version $\texttt{2.0}$ of the open-source $\texttt{C++}$ library $\texttt{FireFly}$ for the interpolation of rational functions. This includes algorithmic improvements, e.g. a hybrid…
We analyze the FEAST method for computing selected eigenvalues and eigenvectors of large sparse matrix pencils. After establishing the close connection between FEAST and the well-known Rayleigh-Ritz method, we identify several critical…
I present an algorithm based on sector decomposition and Mellin-Barnes techniques to power expand Feynman integrals. The coefficients of this expansion are given in terms of finite integrals that can be calculated numerically. I show in an…
In a recent article [1], the FEAST algorithm has been presented as a general purpose eigenvalue solver which is ideally suited for addressing the numerical challenges in electronic structure calculations. Here, FEAST is presented beyond the…
We present GraSSP, a novel approach to perform automated parallelization relying on recent advances in formal verification and synthesis. GraSSP augments an existing sequential program with an additional functionality to decompose data…
The evaluation of multi-loop Feynman integrals is one of the main challenges in the computation of precise theoretical predictions for the cross sections measured at the LHC. In recent years, the method of differential equations has proven…
We present a method for rewriting dimensionally regulated Feynman parameter integrals in the Minkowski regime as a sum of real, positive integrands multiplied by complex prefactors. This representation eliminates the need for contour…
Present and future high-precision tests of the Standard Model and beyond for the fundamental constituents and interactions in Nature are demanding complex perturbative calculations involving multi-leg and multi-loop Feynman diagrams.…
In this paper we study the calculation of multiloop Feynman integrals that cannot be expressed in terms of multiple polylogarithms. We show in detail how certain types of two- and three-point functions at two loops, which appear in the…
Numerical evaluations of Feynman integrals often proceed via a deformation of the integration contour into the complex plane. While valid contours are easy to construct, the numerical precision for a multi-loop integral can depend…