Quasi-Monte Carlo for Partial Differential Equations with Generalized Gaussian Input Uncertainty

SIAM Journal on Numerical Analysis, Volume 63, Issue 4, Page 1666-1690, August 2025.
Abstract. There has been a surge of interest in uncertainty quantification for parametric partial differential equations (PDEs) with Gevrey regular inputs. The Gevrey class contains functions that are infinitely smooth with a growth condition on the higher-order partial derivatives, but which are nonetheless not analytic in general. Recent studies by Chernov and Lê [Comput. Math. Appl., 164 (2024), pp. 116–130; SIAM J. Numer. Anal., 62 (2024), pp. 1874–1900] as well as Harbrecht, Schmidlin, and Schwab [Math. Models Methods Appl. Sci., 34 (2024), pp. 881–917] analyze the setting wherein the input random field is assumed to be uniformly bounded with respect to the uncertain parameters. In this paper, we relax this assumption and allow for parameter-dependent bounds. The parametric inputs are modeled as generalized Gaussian random variables, and we analyze the application of quasi-Monte Carlo (QMC) integration to assess the PDE response statistics using randomly shifted rank-1 lattice rules. In addition to the QMC error analysis, we also consider the dimension truncation and finite element errors in this setting.