Quasi-Monte Carlo and multilevel Monte Carlo methods for computing posterior expectations in elliptic inverse problems

R Scheichl, AM Stuart, AL Teckentrup - SIAM/ASA Journal on Uncertainty …, 2017 - SIAM
We are interested in computing the expectation of a functional of a PDE solution under a
Bayesian posterior distribution. Using Bayes's rule, we reduce the problem to estimating the …

Analysis of circulant embedding methods for sampling stationary random fields

IG Graham, FY Kuo, D Nuyens, R Scheichl… - SIAM Journal on …, 2018 - SIAM
A standard problem in uncertainty quantification and in computational statistics is the
sampling of stationary Gaussian random fields with given covariance in a d-dimensional …

Approximation of high-dimensional periodic functions with Fourier-based methods

D Potts, M Schmischke - SIAM Journal on Numerical Analysis, 2021 - SIAM
We propose an approximation method for high-dimensional 1-periodic functions based on
the multivariate ANOVA decomposition. We provide analysis of classical ANOVA …

Interpretable approximation of high-dimensional data

D Potts, M Schmischke - SIAM Journal on Mathematics of Data Science, 2021 - SIAM
In this paper we apply the previously introduced approximation method based on the
analysis of variance (ANOVA) decomposition and Grouped Transformations to synthetic and …

Multilevel approximation of Gaussian random fields: fast simulation

L Herrmann, K Kirchner, C Schwab - Mathematical Models and …, 2020 - World Scientific
We propose and analyze several multilevel algorithms for the fast simulation of possibly
nonstationary Gaussian random fields (GRFs) indexed, for example, by the closure of a …

Learning multivariate functions with low-dimensional structures using polynomial bases

D Potts, M Schmischke - Journal of Computational and Applied …, 2022 - Elsevier
In this paper we propose a method for the approximation of high-dimensional functions over
finite intervals with respect to complete orthonormal systems of polynomials. An important …

Multilevel quasi-Monte Carlo integration with product weights for elliptic PDEs with lognormal coefficients

L Herrmann, C Schwab - ESAIM: Mathematical Modelling and …, 2019 - esaim-m2an.org
We analyze the convergence rate of a multilevel quasi-Monte Carlo (MLQMC) Finite
Element Method (FEM) for a scalar diffusion equation with log-Gaussian, isotropic …

Quasi-Monte Carlo for partial differential equations with generalized Gaussian input uncertainty

PA Guth, V Kaarnioja - arXiv preprint arXiv:2411.03793, 2024 - arxiv.org
There has been a surge of interest in uncertainty quantification for parametric partial
differential equations (PDEs) with Gevrey regular inputs. The Gevrey class contains …

Generalized dimension truncation error analysis for high-dimensional numerical integration: lognormal setting and beyond

PA Guth, V Kaarnioja - SIAM Journal on Numerical Analysis, 2024 - SIAM
Partial differential equations (PDEs) with uncertain or random inputs have been considered
in many studies of uncertainty quantification. In forward uncertainty quantification, one is …

Quasi-Monte Carlo finite element analysis for wave propagation in heterogeneous random media

M Ganesh, FY Kuo, IH Sloan - SIAM/ASA Journal on Uncertainty Quantification, 2021 - SIAM
We propose and analyze a quasi-Monte Carlo (QMC) algorithm for efficient simulation of
wave propagation modeled by the Helmholtz equation in a bounded region in which the …