Some kinds of eigensolver for large sparse matrices require specification of parameters that are based on rough estimates of the desired eigenvalues. In this paper, we propose a stochastic estimation method of eigenvalue distribution using the combination of a stochastic estimator of the matrix trace and contour integrations. The proposed method can be easily parallelized and applied to matrices for which factorization is infeasible. Numerical experiments are executed to show that the method can perform rough estimates at a low computational cost.
2010 The Japan Society for Industrial and Applied Mathematics