Header

UZH-Logo

Maintenance Infos

A GPU compatible quasi-Monte Carlo integrator interfaced to pySecDec


Borowka, S; Heinrich, G; Jahn, S; Jones, S P; Kerner, M; Schlenk, J (2018). A GPU compatible quasi-Monte Carlo integrator interfaced to pySecDec. CERN 246, University of Zurich.

Abstract

The purely numerical evaluation of multi-loop integrals and amplitudes can be a viable alternative to analytic approaches, in particular in the presence of several mass scales, provided sufficient accuracy can be achieved in an acceptable amount of time. For many multi-loop integrals, the fraction of time required to perform the numerical integration is significant and it is therefore beneficial to have efficient and well-implemented numerical integration methods. With this goal in mind, we present a new stand-alone integrator based on the use of (quasi-Monte Carlo) rank-1 shifted lattice rules. For integrals with high variance we also implement a variance reduction algorithm based on fitting a smooth function to the inverse cumulative distribution function of the integrand dimension-by-dimension. Additionally, the new integrator is interfaced to pySecDec to allow the straightforward evaluation of multi-loop integrals and dimensionally regulated parameter integrals. In order to make use of recent advances in parallel computing hardware, our integrator can be used both on CPUs and CUDA compatible GPUs where available.

Abstract

The purely numerical evaluation of multi-loop integrals and amplitudes can be a viable alternative to analytic approaches, in particular in the presence of several mass scales, provided sufficient accuracy can be achieved in an acceptable amount of time. For many multi-loop integrals, the fraction of time required to perform the numerical integration is significant and it is therefore beneficial to have efficient and well-implemented numerical integration methods. With this goal in mind, we present a new stand-alone integrator based on the use of (quasi-Monte Carlo) rank-1 shifted lattice rules. For integrals with high variance we also implement a variance reduction algorithm based on fitting a smooth function to the inverse cumulative distribution function of the integrand dimension-by-dimension. Additionally, the new integrator is interfaced to pySecDec to allow the straightforward evaluation of multi-loop integrals and dimensionally regulated parameter integrals. In order to make use of recent advances in parallel computing hardware, our integrator can be used both on CPUs and CUDA compatible GPUs where available.

Statistics

Downloads

15 downloads since deposited on 18 Jan 2019
15 downloads since 12 months
Detailed statistics

Additional indexing

Item Type:Working Paper
Communities & Collections:07 Faculty of Science > Physics Institute
Dewey Decimal Classification:530 Physics
Language:English
Date:28 November 2018
Deposited On:18 Jan 2019 07:41
Last Modified:18 Jan 2019 07:47
Series Name:CERN
OA Status:Green

Download

Download PDF  'A GPU compatible quasi-Monte Carlo integrator interfaced to pySecDec'.
Preview
Content: Accepted Version
Filetype: PDF
Size: 850kB