Low-rank solution methods for stochastic eigenvalue problems

From MaRDI portal
Publication:5230664

DOI10.1137/18M122100XzbMATH Open1420.65046arXiv1803.03717OpenAlexW2969984443MaRDI QIDQ5230664FDOQ5230664


Authors: Tengfei Su, Howard C. Elman Edit this on Wikidata


Publication date: 28 August 2019

Published in: SIAM Journal on Scientific Computing (Search for Journal in Brave)

Abstract: We study efficient solution methods for stochastic eigenvalue problems arising from discretization of self-adjoint partial differential equations with random data. With the stochastic Galerkin approach, the solutions are represented as generalized polynomial chaos expansions. A low-rank variant of the inverse subspace iteration algorithm is presented for computing one or several minimal eigenvalues and corresponding eigenvectors of parameter-dependent matrices. In the algorithm, the iterates are approximated by low-rank matrices, which leads to significant cost savings. The algorithm is tested on two benchmark problems, a stochastic diffusion problem with some poorly separated eigenvalues, and an operator derived from a discrete stochastic Stokes problem whose minimal eigenvalue is related to the inf-sup stability constant. Numerical experiments show that the low-rank algorithm produces accurate solutions compared to the Monte Carlo method, and it uses much less computational time than the original algorithm without low-rank approximation.


Full work available at URL: https://arxiv.org/abs/1803.03717




Recommendations




Cites Work


Cited In (14)

Uses Software





This page was built for publication: Low-rank solution methods for stochastic eigenvalue problems

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q5230664)