On the Local Structure and Approximation Stability of Block Isotropic Gaussian Fields
Munki Jeong, Alexander Strang
TL;DR
The authors develop a rigorous framework for local quadratic approximations of block isotropic Gaussian fields, focusing on skew-symmetric constructions where traditional stationarity is incompatible with skew-symmetry. They show that the local geometry at matched inputs decomposes into a Gaussian gradient and Gaussian Hessian blocks (GOE and GSSE components), and that the second-order Taylor approximation remains Gaussian with an explicit covariance structure. By analyzing the error between the original field and its quadratic surrogate, they derive pointwise, regional (ellipsoidal), and distributional uniform bounds, including detailed asymptotics for small and large neighborhoods and dimension-dependent behavior. The results provide practical guidance for applying local quadratic approximations in high-dimensional settings, including how kernel smoothness and input distribution influence error magnitudes and how region sizes must be adapted with dimension to maintain accuracy. Overall, the work offers a tractable, distributional lens for understanding local approximations of skew-symmetric Gaussian fields and their errors in both theory and potential applications.
Abstract
Skew-symmetric functions are a class of functions defined on a product space $M \times M$ that are antisymmetric with respect to the order of their inputs. In [13], the authors proved that non-deterministic skew-symmetric Gaussian fields cannot be stationary or isotropic and proposed an alternative notion: stationarity (isotropy) in each component space. Our work focuses on local quadratic approximations of the associated Gaussian fields. Local quadratic approximations to random fields are random polynomials parametrized by a jointly sampled gradient vector and Hessian matrix. We characterize the distribution of the corresponding random vectors and random matrices. Then, we study the error in the quadratic approximation, which is also a Gaussian field. We investigate the error induced by the quadratic approximation in three senses: the pointwise error, the maximal error over an ellipsoidal region, and the worst-case error for multivariate Gaussian inputs at a given confidence level. Next, we explore the limiting behavior of the worst-case error as the distance between an expansion point and evaluation points approaches zero and infinity. Finally, we study how, as the input dimension increases, the variance of multivariate Gaussian distributions must be restricted to keep the worst-case error bound constant.
