FIESTA 3: Cluster-parallelizable multiloop numerical calculations in physical regions

The goal of this paper is to present a new major release of the program FIESTA (Feynman Integral Evaluation by a Sector decomposiTion Approach). This version presents features like cluster-parallelization, new asymptotic expansion algorithms, calculations in physical regions, new sector-decomposition strategies, as well as multiple speed, memory, and stability improvements. Program title: FIESTA 3 Catalogue identifier: AECP_v3_0 Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AECP_v3_0.html Program obtainable from: CPC Program Library, Queen’s University, Belfast, N. Ireland Licensing provisions: GNU General Public License, version 2 No. of lines in distributed program, including test data, etc.: 109631 No. of bytes in distributed program, including test data, etc.: 6777017 Distribution format: tar.gz Programming language: Wolfram Mathematica 7.0 or higher, C++. Computer: From a desktop PC to a supercomputer. Operating system: Unix, Linux, Mac OS X. Has the code been vectorized or parallelized?: Yes. Number of processors used: from 1 processor up to loading a supercomputer (tests were performed up to 1024 cores) RAM: Depends on the complexity of the problem Catalogue identifier of previous version: AECP_v2_0 Journal reference of previous version: Comput. Phys. Comm. 182 (2011) 790 Classification: 4.4, 4.12, 5, 6.5. External routines: Wolfram Mathematica [1], KyotoCabinet [2], Cuba [3], QHull [4] Does the new version supercede the previous version?: Yes. Some obsolete options were removed as being superseded by better approaches. Nature of problem: The sector decomposition approach to evaluating Feynman integrals falls apart into the sector decomposition itself, where one has to minimize the number of sectors; the pole resolution and epsilon expansion; and the numerical integration of the resulting expression. Moreover, in cases where the integrand is complex, one has to perform a contour deformation. Solution method: The program has a number of sector decomposition strategies. One of the most important features is the ability to perform a contour deformation, as well as the so-called pre-resolution in the case of integrals at the threshold. Everything except the integration is performed in Wolfram Mathematica [1] (required version is 7.0 or higher). This part of the calculation is parallelizable with the use of shared memory. The database is stored on hard disk with the use of the KyotoCabinet [2] database engine. The integration part of the algorith