Towards Monte Carlo preconditioning approach and hybrid Monte Carlo algorithms for Matrix Computations

Towards Monte Carlo preconditioning approach and hybrid Monte Carlo algorithms for Matrix Computations

0.00 Avg rating0 Votes
Article ID: iaor201530429
Volume: 70
Issue: 11
Start Page Number: 2709
End Page Number: 2718
Publication Date: Dec 2015
Journal: Computers and Mathematics with Applications
Authors: ,
Keywords: matrices, numerical analysis
Abstract:

An enhanced version of a stochastic SParse Approximate Inverse (SPAI) preconditioner for general matrices is presented in this paper. This method is used in contrast to the standard deterministic preconditioners computed by the Modified SParse Approximate Inverse Preconditioner (MSPAI). Thus we present a Monte Carlo preconditioner that relies on the use of Markov Chain Monte Carlo (MCMC) methods to compute a rough approximate matrix inverse first, which can further be optimized by an iterative filter process and a parallel refinement, to enhance the accuracy of the inverse and the preconditioner respectively. The advantage of the proposed approach is that finding the sparse Monte Carlo matrix inversion has a computational complexity linear of the size of the matrix, it is inherently parallel and thus can be obtained very efficiently for large matrices and can be used also as an efficient preconditioner while solving systems of linear algebraic equations. The behaviour of the proposed algorithm is studied and its performance measured, evaluated and compared with MSPAI.

Reviews

Required fields are marked *. Your email address will not be published.