Document Type
Article
Publication Date
2012
DOI
10.1016/j.procs.2012.04.041
Publication Title
Procedia Computer Science
Volume
9
Pages
383-392
Abstract
In this paper, we present an approach of reusing random walks in Monte Carlo methods for linear systems. The fundamental idea is, during the Monte Carlo sampling process, the random walks generated to estimate one unknown element can also be effectively reused to estimate the other unknowns in the solution vector. As a result, when the random walks are reused, a single random walk can contribute samples for estimations of multiple unknowns in the solution simultaneously while ensuring that the samples for the same unknown element are statistically independent. Consequently, the total number of random walk transition steps needed for estimating the overall solution vector is reduced, which improves the performance of the Monte Carlo algorithm. We apply this approach to the Monte Carlo algorithm in two linear algebra applications, including solving a system of linear equations and approximating the inversion of a matrix. Our computational results show that compared to the conventional implementations of Monte Carlo algorithms for linear systems without random walk reusing, our approach can significantly improve the performance of Monte Carlo sampling process by reducing the overall number of transition steps in random walks to obtain the entire solution within desired precision.
Rights
© 2012 The Authors. Open Access under the terms of the Creative Commons Attribution-NonCommercial NoDerivs 3.0 Unported License.
Original Publication Citation
Ji, H., & Li, Y. (2012). Reusing random walks in Monte Carlo methods for linear systems. Procedia Computer Science, 9, 383-392. https://doi.org/10.1016/j.procs.2012.04.041
Repository Citation
Ji, H., & Li, Y. (2012). Reusing random walks in Monte Carlo methods for linear systems. Procedia Computer Science, 9, 383-392. https://doi.org/10.1016/j.procs.2012.04.041
ORCID
0000-0001-8303-1491 (Ji)