A highly parallel algorithm for computing the action of a matrix exponential on a vector based on a multilevel Monte Carlo method

04/29/2019
by   Juan A. Acebron, et al.
0

A novel algorithm for computing the action of a matrix exponential over a vector is proposed. The algorithm is based on a multilevel Monte Carlo method, and the vector solution is computed probabilistically generating suitable random paths which evolve through the indices of the matrix according to a suitable probability law. The computational complexity is proved in this paper to be significantly better than the classical Monte Carlo method, which allows the computation of much more accurate solutions. Furthermore, the positive features of the algorithm in terms of parallelism were exploited in practice to develop a highly scalable implementation capable of solving some test problems very efficiently using high performance supercomputers equipped with a large number of cores. For the specific case of shared memory architectures the performance of the algorithm was compared with the results obtained using an available Krylov-based algorithm, outperforming the latter in all benchmarks analyzed so far.

READ FULL TEXT
research
04/29/2019

A Monte Carlo method for computing the action of a matrix exponential on a vector

A Monte Carlo method for computing the action of a matrix exponential fo...
research
08/02/2023

A Fast Monte Carlo algorithm for evaluating matrix functions with application in complex networks

We propose a novel stochastic algorithm that randomly samples entire row...
research
06/17/2021

Scale-invariant multilevel Monte Carlo method

In this paper, the scale-invariant version of the mean and variance mult...
research
03/29/2021

Monte Carlo algorithm for the extrema of tempered stable processes

We develop a novel Monte Carlo algorithm for the vector consisting of th...
research
06/23/2023

Multilevel Monte Carlo methods for the Grad-Shafranov free boundary problem

The equilibrium configuration of a plasma in an axially symmetric reacto...
research
05/25/2023

Computing the Quadratic Numerical Range

A novel algorithm for the computation of the quadratic numerical range i...
research
07/12/2016

Scheduling massively parallel multigrid for multilevel Monte Carlo methods

The computational complexity of naive, sampling-based uncertainty quanti...

Please sign up or login with your details

Forgot password? Click here to reset