# Mitigating the Influence of the Boundary on PDE-based Covariance Operators

arXiv: https://arxiv.org/abs/1610.05280 · HTML: https://mlsanity.com/paper/1610.05280/

## Key facts

- Cited by 4 papers in the map (data through 14 Oct 2026).
- First posted on arXiv 17 Oct 2016.

## Built on by (most cited first, 4 of 4)

- [Bayesian model calibration for block copolymer self-assembly: Likelihood-free inference and expected information gain computation via measure transport](https://mlsanity.com/paper/2206.11343/) (2022)
- [Stochastic PDE representation of random fields for large-scale Gaussian process regression and statistical finite element analysis](https://mlsanity.com/paper/2305.13879/) (2023)
- [Sequential Bayesian Optimal Experimental Design in Infinite Dimensions via Policy Gradient Reinforcement Learning](https://mlsanity.com/paper/2601.05868/) (2026)
- [Enhanced uncertainty quantification variational autoencoders for the solution of Bayesian inverse problems](https://mlsanity.com/paper/2502.13105/) (2025)

## Abstract (arXiv)

Gaussian random fields over infinite-dimensional Hilbert spaces require the definition of appropriate covariance operators. The use of elliptic PDE operators to construct covariance operators allows to build on fast PDE solvers for manipulations with the resulting covariance and precision operators. However, PDE operators require a choice of boundary conditions, and this choice can have a strong and usually undesired influence on the Gaussian random field. We propose two techniques that allow to ameliorate these boundary effects for large-scale problems. The first approach combines the elliptic PDE operator with a Robin boundary condition, where a varying Robin coefficient is computed from an optimization problem. The second approach normalizes the pointwise variance by rescaling the covariance operator. These approaches can be used individually or can be combined. We study properties of these approaches, and discuss their computational complexity. The performance of our approaches is studied for random fields defined over simple and complex two- and three-dimensional domains.