×

Automatic differentiation for error analysis of Monte Carlo data. (English) Zbl 07683953

Summary: Automatic Differentiation (AD) allows to determine exactly the Taylor series of any function truncated at any order. Here we propose to use AD techniques for Monte Carlo data analysis. We discuss how to estimate errors of a general function of measured observables in different Monte Carlo simulations. Our proposal combines the \(\varGamma\)-method with Automatic differentiation, allowing exact error propagation in arbitrary observables, even those defined via iterative algorithms. The case of special interest where we estimate the error in fit parameters is discussed in detail. We also present a freely available reference implementation of the ideas discussed in this work.

MSC:

74-XX Mechanics of deformable solids
65-XX Numerical analysis

Software:

uncertainties; UNEW

References:

[1] Wolff, U., Comput. Phys. Comm., 156, 143-153 (2004), arXiv:hep-lat/0306017 [hep-lat]
[2] Madras, N.; Sokal, A. D., J. Stat. Phys., 50, 1, 109-186 (1988) · Zbl 1084.82503
[3] Schaefer, S.; Sommer, R.; Virotta, F., Nuclear Phys., B845, 93-119 (2011), arXiv:10095228 [hep-lat] · Zbl 1207.81209
[4] Wikipedia contributors, Automatic differentiation — Wikipedia, the free encyclopedia, 2018 (online). https://en.wikipedia.org/wiki/Automatic_differentiation. (Accessed 22 June 2018). · Zbl 1455.65036
[5] Virotta, F., Critical Slowing Down and Error Analysis of Lattice QCD Simulations (2012), Humboldt-Universität zu Berlin, Mathematisch-Naturwissenschaftliche Fakultät I, (Ph.D. thesis)
[6] Lottini, S.; Sommer, R., (Theory and Practice of Data Analysis (2014))
[7] Simma, Hubert; Sommer, Rainer; Virotta, Francesco, General Error Computation in Lattice Gauge Theory (2012), Internal notes ALPHA collaboration
[8] Lebigot, Eric O., (Uncertainties: A Python Package for Calculations with Uncertainties (2018))
[9] B. De Palma, M. Erba, L. Mantovani, N. Mosco, A Python program for the implementation of the \(\Gamma \)-method for Monte Carlo simulations, arXiv:170302766 [hep-lat].
[10] Alberto, Ramos, : Error Analysis of Monte Carlo Data with Automatic Differentiation (2018)
[11] Linnainmaa, S., BIT Numer. Math., 16, 2, 146-160 (1976) · Zbl 0332.65024
[12] Griewank, A., Evaluating Derivatives: Principles and Techniques of Algorithmic Differentiation (2008), Society for Industrial and Applied Mathematics: Society for Industrial and Applied Mathematics Philadelphia, PA · Zbl 1159.65026
[13] J.A. Fike, J.J. Alonso, The development of hyper-dual numbers for exact second-derivative calculations, in: AIAA paper 2011-886, 49th AIAA Aerospace Sciences Meeting, 2011.
[14] Del Debbio, L.; Manca, G. M.; Vicari, E., Phys. Lett., B594, 315-323 (2004), arXiv:hep-lat/0403001 [hep-lat]
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. In some cases that data have been complemented/enhanced by data from zbMATH Open. This attempts to reflect the references listed in the original paper as accurately as possible without claiming completeness or a perfect matching.