Abstract
In this paper we analyze a quasi-Monte Carlo method for solving systems of linear algebraic equations. It is well known that the convergence of Monte Carlo methods for numerical integration can often be improved by replacing pseudorandom numbers with more uniformly distributed numbers known as quasirandom numbers. Here the convergence of a Monte Carlo method for solving systems of linear algebraic equations is studied when quasirandom sequences are used. An error bound is established and numerical experiments with large sparse matrices are performed using Soboĺ, Halton and Faure sequences. The results indicate that an improvement in both the magnitude of the error and the convergence rate are achieved.
Supported by the Ministry of Education and Science of Bulgaria under Grant # MM 902/99 and by Center of Excellence BIS-21 grant ICA1-2000-70016
Access this chapter
Tax calculation will be finalised at checkout
Purchases are for personal use only
Preview
Unable to display preview. Download preview PDF.
References
P. Bratley, B. L. Fox, and H. Niederreiter. Implementation and tests of low discrepancy point sets, ACM Trans. on Modeling and Comp. Simul., 2, 195–213, 1992.
R. E. Caflisch. Monte Carlo and quasi-Monte Carlo methods, Acta Numerica, 7, 1–49, 1998.
I. Dimov and A. Karaivanova. Iterative Monte Carlo algorithms for linear algebra problems, Lecture Notes in Computer Science, 1196, Springer, 66–77, 1996.
I. Dimov, V. Alexandrov, and A. Karaivanova. Resolvent Monte Carlo Methods for linear algebra problems, in The second IMACS Seminar on Monte Carlo Methods, Mathematics and Computers in Simulation, 55(1–3), 25–35, 2001.
H. Faure. Discrépance de suites associées á un système de numération (en dimension s), Acta Arithmetica, XLI, 337–351, 1992.
J. H. Halton. Sequential Monte Carlo techniques for the solution of linear systems, TR 92-033, University of North Carolina at Chapel Hill, Department of Computer Science, 1992.
J. H. Halton. On the efficiency of certain quasi-random sequences of points in evaluating multi-dimensional integrals, Numer. Math., 2, 84–90, 1960.
J. M. Hammersley and D. C. Handscomb. Monte Carlo Methods, Chapman and Hall, London, New York, 1964.
M. Mascagni and A. Karaivanova. Are quasirandom numbers good for anything besides integration?, in Proc. of Advances in Reactor Physics and Mathematics and Computation into the Next Millennium (PHYSOR2000), 2000, (to appear).
M. Mascagni and A. Karaivanova. Matrix computations using quasirandom sequences, in Numerical Analysis and its Applications, Lecture Notes in Computer Science, 1988, Springer, 552–559, 2001.
H. Niederreiter. Random Number Generation and Quasi-Monte Carlo Methods, SIAM, Philadelphia, 1992.
I. M. Soboĺ. Monte Carlo Numerical Methods, Nauka, Moscow, 1973, (in Russian).
I. M. Soboĺ. The distribution of points in a cube and approximate evaluation of integrals, Zh. Vychisl. Mat. Mat. Fiz., 7, 784–802, 1967, (in Russian).
Author information
Authors and Affiliations
Editor information
Editors and Affiliations
Rights and permissions
Copyright information
© 2001 Springer-Verlag Berlin Heidelberg
About this paper
Cite this paper
Karaivanova, A., Georgieva, R. (2001). Solving Systems of Linear Algebraic Equations Using Quasirandom Numbers. In: Margenov, S., Waśniewski, J., Yalamov, P. (eds) Large-Scale Scientific Computing. LSSC 2001. Lecture Notes in Computer Science, vol 2179. Springer, Berlin, Heidelberg. https://doi.org/10.1007/3-540-45346-6_16
Download citation
DOI: https://doi.org/10.1007/3-540-45346-6_16
Published:
Publisher Name: Springer, Berlin, Heidelberg
Print ISBN: 978-3-540-43043-8
Online ISBN: 978-3-540-45346-8
eBook Packages: Springer Book Archive