Analysis of Circulant Embedding Methods for Sampling Stationary Random Fields

DOI10.1137/17M1149730zbMATH Open1434.60114arXiv1710.00751OpenAlexW2964295897WikidataQ129646751 ScholiaQ129646751MaRDI QIDQ4572025FDOQ4572025

I. G. Graham, F. Y. Kuo, Robert Scheichl, Dirk Nuyens, Ian H. Sloan

Publication date: 4 July 2018

Published in: SIAM Journal on Numerical Analysis (Search for Journal in Brave)

Abstract: In this paper we prove, under mild conditions, that the positive definiteness of the circulant matrix appearing in the circulant embedding method is always guaranteed, provided the enclosing cube is sufficiently large. We examine in detail the case of the Mat'ern covariance, and prove (for fixed correlation length) that, as $h_{0}ightarrow0$, positive definiteness is guaranteed when the random field is sampled on a cube of size order $(1+u^{1/2}log{h}_0^{-1})$ times larger than the size of the physical domain. (Here $h_0$ is the mesh spacing of the regular grid and $u$ the Mat'ern smoothness parameter.) We show that the sampling cube can become smaller as the correlation length decreases when $h_0$ and $u$ are fixed. Our results are confirmed by numerical experiments. We prove several results about the decay of the eigenvalues of the circulant matrix. These lead to the conjecture, verified by numerical experiment, that they decay with the same rate as the Karhunen--Lo`{e}ve eigenvalues of the covariance operator.

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

zbMATH Keywords

Gaussian random fieldsfast Fourier transformcirculant embeddingstatistical homogeneityMatérn covariance

Mathematics Subject Classification ID

Computational methods for problems pertaining to statistics (62-08) Monte Carlo methods (65C05) Random fields (60G60) Stationary stochastic processes (60G10)

Cites Work

• An Algorithm for Simulating Stationary Gaussian Random Fields
• Title not available (Why is that?)
• Title not available (Why is that?)
• Scattered Data Approximation
• Application of hierarchical matrices for computing the Karhunen-Loève expansion
• Karhunen-Loève approximation of random fields by generalized fast multipole methods
• An introduction to computational stochastic PDEs
• Title not available (Why is that?)
• Fast and Exact Simulation of Stationary Gaussian Processes through Circulant Embedding of the Covariance Matrix
• Title not available (Why is that?)
• Quasi-Monte Carlo methods for elliptic PDEs with random coefficients and applications
• Some Extensions of W. Gautschi's Inequalities for the Gamma Function
• Computational aspects of the stochastic finite element method
• Efficient approximation of random fields for numerical applications
• Quasi-Monte Carlo finite element methods for elliptic PDEs with lognormal random coefficients
• Representations of Gaussian random fields and approximation of elliptic PDEs with lognormal coefficients
• Fast random field generation with \(H\)-matrices

Cited In (28)

• Fast sampling of parameterised Gaussian random fields
• Goal-oriented adaptive finite element multilevel Monte Carlo with convergence rates
• Circulant embedding with QMC: analysis for elliptic PDE with lognormal coefficients
• Multilevel quasi-Monte Carlo for optimization under uncertainty
• A Statistical Framework for Generating Microstructures of Two-Phase Random Materials: Application to Fatigue Analysis
• Multilevel Hierarchical Decomposition of Finite Element White Noise with Application to Multilevel Markov Chain Monte Carlo
• Multilevel Hierarchical Decomposition of Finite Element White Noise with Application to Multilevel Markov Chain Monte Carlo
• Multilevel Quasi Monte Carlo Methods for Elliptic PDEs with Random Field Coefficients via Fast White Noise Sampling
• Unified Analysis of Periodization-Based Sampling Methods for Matérn Covariances
• Parallel Multilevel Monte Carlo Algorithms for Elliptic PDEs with Random Coefficients
• Hurst exponent estimation of fractional surfaces for mammogram images analysis
• Analysis of Boundary Effects on PDE-Based Sampling of Whittle--Matérn Random Fields
• Advances in Gaussian random field generation: a review
• Periodic vector processes with an internal reciprocal dynamics
• Finite elements for Matérn-type random fields: uncertainty in computational mechanics and design optimization
• Recycling Samples in the Multigrid Multilevel (Quasi-)Monte Carlo Method
• Multilevel Monte Carlo estimators for elliptic PDEs with Lévy-type diffusion coefficient
• On properties and applications of Gaussian subordinated Lévy fields
• Stochastic modeling of geometrical uncertainties on complex domains, with application to additive manufacturing and brain interface geometries
• Multilevel quasi-Monte Carlo integration with product weights for elliptic PDEs with lognormal coefficients
• Sparse Cholesky Factorization by Kullback--Leibler Minimization
• A Multilevel Monte Carlo Algorithm for Parabolic Advection-Diffusion Problems with Discontinuous Coefficients
• Initialization of the circulant embedding method to speed up the generation of Gaussian random fields
• Subordinated Gaussian random fields in elliptic partial differential equations
• Optimization-Based Markov Chain Monte Carlo Methods for Nonlinear Hierarchical Statistical Inverse Problems
• Regularity and convergence analysis in Sobolev and Hölder spaces for generalized Whittle-Matérn fields
• Kryging: geostatistical analysis of large-scale datasets using Krylov subspace methods
• A fast direct solver for nonlocal operators in wavelet coordinates

Uses Software

• AS 312

   Recommendations

• Initialization of the circulant embedding method to speed up the generation of Gaussian random fields 👍 👎
• Spectral and circulant approximations to the likelihood for stationary Gaussian random fields 👍 👎
• Spectral density estimation for random fields via periodic embeddings 👍 👎
• Title not available (Why is that?) 👍 👎
• Fast and Exact Simulation of Stationary Gaussian Processes through Circulant Embedding of the Covariance Matrix 👍 👎
• Sampling From Gaussian Markov Random Fields Using Stationary and Non-Stationary Subgraph Perturbations 👍 👎
• Title not available (Why is that?) 👍 👎
• Conditional sampling for spectrally discrete max-stable random fields 👍 👎

This page was built for publication: Analysis of Circulant Embedding Methods for Sampling Stationary Random Fields

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