Distributed and multi-core computation of 2-loop integrals

For an automatic computation of Feynman loop integrals in the physical region we rely on an extrapolation technique where the integrals of the sequence are obtained with iterated/repeated adaptive methods from the QUADPACK 1D quadrature package. The integration rule evaluations in the outer level, c...

Full description

Saved in:
Bibliographic Details
Published in:Journal of physics. Conference series 2014-06, Vol.523 (1), p.12052-12
Main Authors: de Doncker, E, Yuasa, F
Format: Article
Language:English
Subjects:
Algorithms
Approximation
Computation
Extrapolation
Integral transforms
Integrals
Microprocessors
Packages
Physics
Quadratures
Singularities
Citations: Items that this one cites
Items that cite this one
Online Access:Get full text
Tags: Add Tag
No Tags, Be the first to tag this record!
Description
Summary:For an automatic computation of Feynman loop integrals in the physical region we rely on an extrapolation technique where the integrals of the sequence are obtained with iterated/repeated adaptive methods from the QUADPACK 1D quadrature package. The integration rule evaluations in the outer level, corresponding to independent inner integral approximations, are assigned to threads dynamically via the OpenMP runtime in the parallel implementation. Furthermore, multi-level (nested) parallelism enables an efficient utilization of hyperthreading or larger numbers of cores. For a class of loop integrals in the unphysical region, which do not suffer from singularities in the interior of the integration domain, we find that the distributed adaptive integration methods in the multivariate PARINT package are highly efficient and accurate. We apply these techniques without resorting to integral transformations and report on the capabilities of the algorithms and the parallel performance for a test set including various types of two-loop integrals.
ISSN:1742-6588
1742-6596
DOI:10.1088/1742-6596/523/1/012052