We begin with isotropic Gaussian random fields, and show how the Bochner–Godement theorem gives a natural way to describe their covariance structure. We continue with a study of Matérn processes on Euclidean space, spheres, manifolds and graphs, using Bessel potentials and stochastic partial differential equations (SPDEs). We then turn from this continuous setting to approximating discrete settings, Gaussian Markov random fields (GMRFs), and the computational advantages they bring in handling large data sets, by exploiting the sparseness properties of the relevant precision (concentration) matrices.