Abstract:Gaussian random fields are popular models for spatially varying uncertainties, arising for instance in geotechnical engineering, hydrology or image processing. A Gaussian random field is fully characterised by its mean function and covariance operator. In more complex models these can also be partially unknown. In this case we need to handle a family of Gaussian random fields indexed with hyperparameters. Sampling for a fixed configuration of hyperparameters is already very expensive due to the nonlocal nature of many classical covariance operators. Sampling from multiple configurations increases the total computational cost severely. In this report we employ parameterised Karhunen-Loève expansions for sampling. To reduce the cost we construct a reduced basis surrogate built from snapshots of Karhunen-Loève eigenvectors. In particular, we consider Matérn-type covariance operators with unknown correlation length and standard deviation. We suggest a linearisation of the covariance function and describe the associated online-offline decomposition. In numerical experiments we investigate the approximation error of the reduced eigenpairs. As an application we consider forward uncertainty propagation and Bayesian inversion with an elliptic partial differential equation where the logarithm of the diffusion coefficient is a parameterised Gaussian random field. In the Bayesian inverse problem we employ Markov chain Monte Carlo on the reduced space to generate samples from the posterior measure. All numerical experiments are conducted in 2D physical space, with non-separable covariance operators, and finite element grids with $\sim 10^4$ degrees of freedom.

gapsplit: Efficient random sampling for non-convex constraint-based models

HOPS: high-performance library for (non-)uniform sampling of convex-constrained models

Random Multi-Constraint Projection: Stochastic Gradient Methods for Convex Optimization with Many Constraints

Randomized subspace gradient method for constrained optimization

Gradient and Uncertainty Enhanced Sequential Sampling for Global Fit

A Randomized Nonlinear Rescaling Method in Large-Scale Constrained Convex Optimization

Randomised subspace methods for non-convex optimization, with applications to nonlinear least-squares

A Fast, Robust Elliptical Slice Sampling Implementation for Linearly Truncated Multivariate Normal Distributions

Learning to sample fibers for goodness-of-fit testing

Latin hypercube sampling with inequality constraints

Sample-and-Bound for Non-Convex Optimization

Random minibatch subgradient algorithms for convex problems with functional constraints

Large-scale nonconvex optimization: randomization, gap estimation, and numerical resolution

Random constraint sampling and duality for convex optimization

A New Split Algorithm for 3D Gaussian Splatting

An efficient multi-core implementation of a novel HSS-structured multifrontal solver using randomized sampling

Constrained Sampling with Primal-Dual Langevin Monte Carlo

hopsy — a methods marketplace for convex polytope sampling in Python

Fast sampling of parameterised Gaussian random fields

GAP: A General Framework for Information Pooling in Two-Sample Sparse Inference

P-split formulations: A class of intermediate formulations between big-M and convex hull for disjunctive constraints