A Monte Carlo Approach to Sparse Approximate Inverse Matrix Computations
Author(s) -
Janko Straßburg,
Vassil Alexandrov
Publication year - 2013
Publication title -
procedia computer science
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 0.334
H-Index - 76
ISSN - 1877-0509
DOI - 10.1016/j.procs.2013.05.402
Subject(s) - monte carlo method , computer science , markov chain monte carlo , hybrid monte carlo , quasi monte carlo method , mathematical optimization , inverse , algorithm , monte carlo integration , matrix exponential , matrix norm , mathematics , matrix (chemical analysis) , monte carlo molecular modeling , statistical physics , statistics , geometry , materials science , composite material , mathematical analysis , eigenvalues and eigenvectors , physics , quantum mechanics , differential equation
In this paper we present a stochastic SPAI pre-conditioner. In contrast to the standard deterministic SPAI pre-conditioners that use the Frobenius norm, we present a Monte Carlo pre-conditioner that relies on the use of Markov Chain Monte Carlo methods to compute a rough matrix inverse (MI). Monte Carlo methods quantify the uncertainties by enabling us to estimate the non-zero elements of the inverse matrix with a given precision and certain probability. The advantage of this approach is that we use sparse Monte Carlo matrix inversion whose complexity is linear to the size of the matrix. The behaviour of the proposed algorithm is studied, its performance is measured and compared with the standard deterministic SPAI approach, as well as the optimized and parallel MSPAI version. An analysis of the results is also presented
Accelerating Research
Robert Robinson Avenue,
Oxford Science Park, Oxford
OX4 4GP, United Kingdom
Address
John Eccles HouseRobert Robinson Avenue,
Oxford Science Park, Oxford
OX4 4GP, United Kingdom