Sparse Compression of Expected Solution Operators
Numerical Analysis
2020-03-17 v2 Numerical Analysis
Abstract
We show that the expected solution operator of prototypical linear elliptic partial differential operators with random coefficients is well approximated by a computable sparse matrix. This result is based on a random localized orthogonal multiresolution decomposition of the solution space that allows both the sparse approximate inversion of the random operator represented in this basis as well as its stochastic averaging. The approximate expected solution operator can be interpreted in terms of classical Haar wavelets. When combined with a suitable sampling approach for the expectation, this construction leads to an efficient method for computing a sparse representation of the expected solution operator.
Cite
@article{arxiv.1807.01741,
title = {Sparse Compression of Expected Solution Operators},
author = {Michael Feischl and Daniel Peterseim},
journal= {arXiv preprint arXiv:1807.01741},
year = {2020}
}