2014/07/28 by Dan Mazur, Jeremy Heyl, Mazur, Dan +2 · 1 citation
Computer Science · Mathematics · Physics and Astronomy · #Computational Physics and Python Applications #FOS: Physical sciences #High Energy Physics - Theory (hep-th) #Scientific Research and Discoveries #Tensor decomposition and applications #hep-th
paper · pdf · doi:10.48550/arxiv.1407.7486
17 pages, 12 figures
arxiv created 2014/07/28 · openalex publication_date 2014/07/28 · arxiv updated 2014/07/29 · openalex created_date 2025/10/10 · openalex updated_date 2026/07/28
We give an overview of the worldline numerics technique, and discuss the parallel CUDA implementation of a worldline numerics algorithm. In the worldline numerics technique, we wish to generate an ensemble of representative closed-loop particle trajectories, and use these to compute an approximate average value for Wilson loops. We show how this can be done with a specific emphasis on cylindrically symmetric magnetic fields. The fine-grained, massive parallelism provided by the GPU architecture results in considerable speedup in computing Wilson loop averages. Furthermore, we give a brief overview of uncertainty analysis in the worldline numerics method. There are uncertainties from discretizing each loop, and from using a statistical ensemble of representative loops. The former can be minimized so that the latter dominates. However, determining the statistical uncertainties is complicated by two subtleties. Firstly, the distributions generated by the worldline ensembles are highly non-Gaussian, and so the standard error in the mean is not a good measure of the statistical uncertainty. Secondly, because the same ensemble of worldlines is used to compute the Wilson loops at different values of T and x_ cm, the uncertainties associated with each computed value of the integrand are strongly correlated. We recommend a form of jackknife analysis which deals with both of these problems.