Analysis of the convergence of the 1/t and Wang-Landau algorithms in the calculation of multidimensional integrals
Abstract
In this Brief Report, the convergence of the 1/t and Wang-Landau algorithms in the calculation of multidimensional numerical integrals is analyzed. Both simulation methods are applied to a wide variety of integrals without restrictions in one, two, and higher dimensions. The efficiency and accuracy of both algorithms are determined by the dynamical behavior of the errors between the exact and the calculated values of the integral. It is observed that the time dependence of the error calculated with the 1/t algorithm varies as N-1/2 [with N the number of Monte Carlo (MC) trials], in quantitative agreement with the simple sampling Monte Carlo method. In contrast, the error calculated with the Wang-Landau algorithm saturates in time, evidencing the nonconvergence of this method. The sources of error for both methods are also determined.
- Publication:
-
Physical Review E
- Pub Date:
- December 2008
- DOI:
- 10.1103/PhysRevE.78.067701
- arXiv:
- arXiv:0806.0268
- Bibcode:
- 2008PhRvE..78f7701B
- Keywords:
-
- 02.60.Jh;
- 05.10.Ln;
- 02.70.Uu;
- Numerical differentiation and integration;
- Monte Carlo methods;
- Applications of Monte Carlo methods;
- Condensed Matter - Statistical Mechanics
- E-Print:
- 8 pages, 5 figures