Randomized block Krylov space methods for trace and log-determinant estimators

02/29/2020
by   Hanyu Li, et al.
0

We present randomized algorithms based on block Krylov space method for estimating the trace and log-determinant of Hermitian positive semi-definite matrices. Using the properties of Chebyshev polynomial and Gaussian random matrix, we provide the error analysis of the proposed estimators and obtain the expectation and concentration error bounds. These bounds improve the corresponding ones given in the literature. Numerical experiments are presented to illustrate the performance of the algorithms and to test the error bounds.

READ FULL TEXT

page 1

page 2

page 3

page 4

research
12/19/2022

Computation of the von Neumann entropy of large matrices via trace estimators and rational Krylov methods

We consider the problem of approximating the von Neumann entropy of a la...
research
07/05/2023

Suboptimal subspace construction for log-determinant approximation

Variance reduction is a crucial idea for Monte Carlo simulation and the ...
research
05/03/2023

Estimating the error in CG-like algorithms for least-squares and least-norm problems

In [Meurant, Papež, Tichý; Numerical Algorithms 88, 2021], we presented ...
research
05/20/2020

On randomized trace estimates for indefinite matrices with an application to determinants

Randomized trace estimation is a popular and well studied technique that...
research
12/10/2020

A generalised log-determinant regularizer for online semi-definite programming and its applications

We consider a variant of online semi-definite programming problem (OSDP)...
research
02/06/2022

Monte Carlo Methods for Estimating the Diagonal of a Real Symmetric Matrix

For real symmetric matrices that are accessible only through matrix vect...
research
01/03/2018

Randomized Linear Algebra Approaches to Estimate the Von Neumann Entropy of Density Matrices

The von Neumann entropy, named after John von Neumann, is the extension ...

Please sign up or login with your details

Forgot password? Click here to reset