Randomized vector iterative linear solvers of high precision for large dense system

Abstract

Abstract In this paper we suggest randomized linear solvers with a focus on refinement issue to achieve a high precision while maintaining all the advantages of the Monte Carlo method for solving systems of large dimension with dense matrices. It is shown that each iterative refinement step reduces the error by one order of magnitude. The crucial point of the suggested method is, in contrast to the standard Monte Carlo method, that the randomized vector algorithm computes the entire solution column at once, rather than a single component. This makes it possible to efficiently construct the iterative refinement method. We apply the developed method for solving a system of elasticity equations.