Matrix computations (3rd ed.)
A new iterative Monte Carlo approach for inverse matrix problem
Journal of Computational and Applied Mathematics
Parallel resolvent Monte Carlo algorithms for linear algebra problems
Mathematics and Computers in Simulation - IMACS sponsored Special issue on the second IMACS seminar on Monte Carlo methods
Parallel Monte Carlo Algorithms for Sparse SLAE Using MPI
Proceedings of the 6th European PVM/MPI Users' Group Meeting on Recent Advances in Parallel Virtual Machine and Message Passing Interface
Hybrid Monte Carlo Methods for Matrix Computation
NMA '02 Revised Papers from the 5th International Conference on Numerical Methods and Applications
Coarse Grained Parallel Monte Carlo Algorithms for Solving SLAE Using PVM
Proceedings of the 5th European PVM/MPI Users' Group Meeting on Recent Advances in Parallel Virtual Machine and Message Passing Interface
ICCS '02 Proceedings of the International Conference on Computational Science-Part II
Hybrid Monte Carlo Methods for Matrix Computation
NMA '02 Revised Papers from the 5th International Conference on Numerical Methods and Applications
Monte Carlo methods for matrix computations on the grid
Future Generation Computer Systems
A sparse parallel hybrid monte carlo algorithm for matrix computations
ICCS'05 Proceedings of the 5th international conference on Computational Science - Volume Part III
Parallel hybrid monte carlo algorithms for matrix computations
ICCS'05 Proceedings of the 5th international conference on Computational Science - Volume Part III
Finding the smallest eigenvalue by the inverse monte carlo method with refinement
ICCS'05 Proceedings of the 5th international conference on Computational Science - Volume Part III
On scalability behaviour of Monte Carlo sparse approximate inverse for matrix computations
ScalA '13 Proceedings of the Workshop on Latest Advances in Scalable Algorithms for Large-Scale Systems
Hi-index | 0.00 |
In this paper we consider mixed (fast stochastic approximation and deterministic refinement) algorithms for Matrix Inversion (MI) and Solving Systems of Linear Equations (SLAE). Monte Carlo methods are used for the stochastic approximation, since it is known that they are very efficient in finding a quick rough approximation of the element or a row of the inverse matrix or finding a component of the solution vector. In this paper we show how the stochastic approximation of the MI can be combined with a deterministic refinement procedure to obtain MI with the required precision and further solve the SLAE using MI. We employ a splitting A = D - C of a given non-singular matrix A, where D is a diagonal dominant matrix and matrix C is a diagonal matrix. In our algorithm for solving SLAE and MI different choices of D can be considered in order to control the norm of matrix T = D-1C; of the resulting SLAE and to minimize the number of the Markov Chains required to reach given precision. Experimental results with dense and sparse matrices are presented.