Galerkin--Chebyshev approximation of Gaussian random fields on compact Riemannian manifolds

Abstract

A new numerical approximation method for a class of Gaussian random fields on compact connected oriented Riemannian manifolds is introduced. This class of random fields is characterized by the Laplace--Beltrami operator on the manifold. A Galerkin approximation is combined with a polynomial approximation using Chebyshev series. This so-called Galerkin--Chebyshev approximation scheme yields efficient and generic sampling algorithms for Gaussian random fields on manifolds. Strong and weak orders of convergence for the Galerkin approximation and strong convergence orders for the Galerkin--Chebyshev approximation are shown and confirmed through numerical experiments.

arXiv commentVersion submitted to a peer-reviewed journal. Changes: fixed residual typos, new outline. 33 pages, 5 figures

References in corpus (7)

• Strictly and non-strictly positive definite functions on spheres
• Isotropic Gaussian random fields on the sphere: Regularity, fast simulation and stochastic partial differential equations
• Numerical solution of fractional elliptic stochastic PDEs with spatial white noise
• The rational SPDE approach for Gaussian random fields with general smoothness
• Weak convergence of Galerkin approximations for fractional elliptic stochastic PDEs with spatial white noise
• Fast generation of isotropic Gaussian random fields on the sphere
• Surface finite element approximation of spherical Whittle--Matérn Gaussian random fields