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.

Stable evaluation of Gaussian radial basis functions using Hermite polynomials

Stabilized interpolation using radial basis functions augmented with selected radial polynomials

On a generalized Gaussian radial basis function: Analysis and applications

Local Error Estimates for Radial Basis Function Interpolation of Scattered Data

Extending error bounds for radial basis function interpolation to measuring the error in higher order Sobolev norms

An algorithmic approach to Hermite-Birkhoff-interpolation

A Discrete Adapted Hierarchical Basis Solver For Radial Basis Function Interpolation

On Fast Implementation of Higher Order Hermite-Fejer Interpolation

Radial polynomials as alternatives to smooth radial basis functions and their applications

Hermite-Birkhoff interpolation on arbitrarily distributed data on the sphere and other manifolds

Stability estimates for radial basis function methods applied to linear scalar conservation laws

Hermite–Birkhoff interpolation on scattered data on the sphere and other manifolds

Stable and Accurate Least Squares Radial Basis Function Approximations on Bounded Domains

Hermite-Birkhoff interpolation on scattered data on the sphere and other manifolds

Univariate interpolation by exponential functions and gaussian RBFs for generic sets of nodes

Improved Laguerre Spectral Methods with Less Round-off Errors and Better Stability

The Role of Hilbert–Schmidt SVD basis in Hermite–Birkhoff interpolation in fractional sense

Fast sampling of parameterised Gaussian random fields

Hermite-Birkhoff interpolation of scattered data by radial basis functions