Efficient generation of Gaussian random fields on metric graphs via domain decomposition and mass matrix lumping
For researchers and practitioners using Gaussian random fields on metric graphs, this method enables scalable sampling on large graphs without sacrificing accuracy.
This work tackles the computational bottleneck of sampling Gaussian random fields on large metric graphs by combining Neumann-Neumann domain decomposition with mass matrix lumping, achieving multi-order speedups and massive memory reductions while preserving exact theoretical convergence rates.
We consider Gaussian Random Fields on metric graphs defined implicitly as the stationary solution to a fractional SPDE driven by Gaussian white noise. Sampling from the finite element approximation requires the Cholesky factorization of the mass matrix, causing non-linear execution time explosions and massive memory fill-in on large graphs. Hence, we combine Neumann-Neumann graph decomposition with mass matrix lumping and demonstrate empirically, that our approach preserves exact theoretical convergence rates established in [8] while achieving multi-order speedups and massive memory reductions.