Sufficient conditions for QMC analysis of finite elements for parametric differential equations

Andreas Rupp; Jay Gopalakrishnan; Vesa Kaarnioja

arxiv: 2511.01703 · v2 · submitted 2025-11-03 · 🧮 math.NA · cs.NA

Sufficient conditions for QMC analysis of finite elements for parametric differential equations

Vesa Kaarnioja , Andreas Rupp , Jay Gopalakrishnan This is my paper

Pith reviewed 2026-05-18 01:11 UTC · model grok-4.3

classification 🧮 math.NA cs.NA

keywords Quasi-Monte CarloFinite element methodsParametric PDEsParametric regularityFlux approximationDiffusion equationGevrey regularityUncertainty quantification

0 comments

The pith

Quasi-Monte Carlo methods achieve dimension-independent convergence faster than Monte Carlo for quantities of interest that depend continuously on the primal variable, flux, or gradient in finite element discretizations of parametric PDEs.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper establishes sufficient conditions under which quasi-Monte Carlo methods applied to finite element discretizations of parametric differential equations achieve convergence rates that do not deteriorate as the number of parameters increases and are better than those of plain Monte Carlo sampling. It models the random parameter as a Gevrey-regular function of infinitely many independent random variables and derives error bounds for quantities of interest that depend continuously on the primal solution, its flux, or its gradient. These conditions are verified for standard conforming elements, mixed methods, and hybridizable discontinuous Galerkin schemes applied to the diffusion equation. If the claims hold, practitioners could use QMC to compute statistics of fluxes and gradients in high-dimensional parametric problems with fewer samples than Monte Carlo requires. The analysis highlights the importance of accurate flux approximation in the discretizations.

Core claim

The authors prove that under assumptions ensuring parametric regularity of the flux vector field in discretizations of balance laws with Gevrey-regular random parameters expressible via sequences of independent random variables, the QMC quadrature error for a quantity of interest is bounded by a rate independent of dimension and superior to Monte Carlo, whenever that quantity depends continuously on the primal solution, the flux, or the gradient. The assumptions are verified for the diffusion equation discretized by conforming finite elements, mixed finite elements, and hybridizable discontinuous Galerkin methods.

What carries the argument

Parametric regularity estimates for the discrete flux vector field, obtained from the Gevrey regularity of the random parameter modeled as a function of countably infinite sequences of independent uniform or normal random variables.

If this is right

If the quantity of interest depends continuously on the flux, the QMC convergence rate is dimension-independent and faster than Monte Carlo.
The same dimension-independent improvement holds when the quantity depends continuously on the primal variable or its gradient.
Conforming finite elements, mixed methods, and hybridizable discontinuous Galerkin schemes for the stationary diffusion equation satisfy the required assumptions.
Accurate approximation of the flux in the discretization is necessary to realize the improved QMC rates.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

The same regularity framework could support QMC analysis for other linear parametric PDEs such as linear elasticity or Stokes flow.
Flux-recovery post-processing techniques may become more valuable in finite element codes aimed at uncertainty quantification tasks.
The approach might extend to mildly nonlinear problems if the parametric regularity estimates can be carried through the nonlinearity.

Load-bearing premise

The random parameter must be modeled as a Gevrey-regular random field expressible as a function of countably infinite sequences of independent random variables.

What would settle it

A numerical test on a discretization satisfying the paper's assumptions where a quantity of interest depends continuously on the flux but the observed QMC convergence rate fails to exceed the Monte Carlo rate or deteriorates as the parameter dimension grows.

Figures

Figures reproduced from arXiv: 2511.01703 by Andreas Rupp, Jay Gopalakrishnan, Vesa Kaarnioja.

read the original abstract

Parametric regularity of discretizations of flux vector fields satisfying a balance law is studied under some assumptions on a random parameter that links the flux with an unknown primal variable (often through a constitutive law). In the primary example of the stationary diffusion equation, the parameter corresponds to the inverse of the diffusivity. The random parameter is modeled here as a Gevrey-regular random field. Specific focus is on random fields expressible as functions of countably infinite sequences of independent random variables, which may be uniformly or normally distributed. Quasi-Monte Carlo (QMC) error bounds for some quantity of interest that depends on the flux are then derived using the parametric regularity. It is shown that the QMC method achieves a dimension-independent, faster-than-Monte Carlo convergence rate if the quantity of interest depends continuously on the primal variable, its flux, or its gradient. A series of assumptions are introduced with the goal of encompassing a broad class of discretizations by various finite element methods. The assumptions are verified for the diffusion equation discretized using conforming finite elements, mixed methods, and hybridizable discontinuous Galerkin schemes. Numerical experiments confirm the analytical findings, highlighting the role of accurate flux approximation in QMC methods.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

This paper gives a reusable set of assumptions that let QMC bounds apply to conforming, mixed, and HDG discretizations of parametric diffusion problems, with explicit checks and some numerics.

read the letter

This paper gives sufficient conditions that let QMC work on a range of finite element methods for parametric diffusion equations. The authors set up assumptions on how the discrete solutions depend on the random parameter, then show these hold for conforming, mixed, and HDG discretizations of the stationary diffusion problem. They derive QMC convergence rates that are dimension-independent and better than Monte Carlo when the quantity of interest is continuous in the primal variable, flux, or gradient. Numerical experiments are included to check the findings. What stands out is the explicit verification for mixed and hybridizable discontinuous Galerkin methods. Prior work covered conforming elements more often, so this broadens the practical reach for uncertainty quantification in engineering problems where those methods are preferred for flux accuracy. The potential weak point is whether the Gevrey regularity carries through with constants that do not grow with the mesh size. The stress-test note points out that inverse estimates or local Lipschitz conditions on the constitutive law could introduce h-dependent factors in the mixed partial derivatives. If the paper's proofs keep those factors bounded independently of h, the claims hold; otherwise the rates might degrade for refined meshes. The abstract claims verification, but the strength depends on how they control the parameter dependence in the discrete operators. This work is aimed at researchers combining quasi-Monte Carlo with finite element methods for parametric PDEs. Someone already familiar with the regularity theory for random fields will find the discretization assumptions the useful new piece. It deserves a serious referee because the framework is concrete and the extension to additional methods is substantive. I would recommend sending it for peer review.

Referee Report

3 major / 2 minor

Summary. The manuscript derives sufficient conditions for parametric regularity of finite-element discretizations of flux vector fields satisfying a balance law, with the random parameter modeled as a Gevrey-regular function of countably many i.i.d. uniform or normal random variables. It obtains QMC error bounds for quantities of interest that depend continuously on the primal variable, its flux, or its gradient, establishing dimension-independent convergence rates faster than Monte Carlo. The assumptions are verified for the stationary diffusion equation discretized by conforming finite elements, mixed methods, and hybridizable discontinuous Galerkin schemes, with numerical experiments supporting the analysis.

Significance. If the central claims hold, the work supplies a general framework for QMC analysis of parametric PDEs that covers multiple finite-element families and extends to flux- and gradient-dependent quantities of interest. The explicit verification of the assumptions across discretization types and the accompanying numerical confirmation of the predicted rates constitute a concrete contribution to the analysis of high-dimensional parametric problems.

major comments (3)

[§3] §3 (parametric regularity assumptions): The propagation of Gevrey regularity from the random field to the discrete flux and gradient requires that all stability and approximation constants in the discrete balance law remain independent of the mesh size h. The manuscript does not explicitly bound the h-dependence arising from inverse estimates in the flux space; if such factors appear, the mixed partial derivatives grow with negative powers of h, which would make the QMC weights h-dependent and destroy both the dimension-independent rate and the faster-than-MC claim for any fixed discretization.
[§4.3] §4.3 (verification for HDG schemes): The claim that the discrete flux inherits Gevrey-type bounds uniform in h relies on the hybridizable formulation preserving the required Lipschitz and differentiability properties of the constitutive map without h-dependent constants. No explicit estimate is given showing that the stabilization parameter and the trace inequalities used in the HDG flux reconstruction remain analytic with h-independent radii of analyticity.
[§5] §5 (QMC error bounds): The derivation of the QMC convergence rate assumes that the parametric regularity constants are independent of the discretization parameter h once the mesh is chosen to balance discretization and sampling errors. If the constants actually deteriorate as O(h^{-k}) for some k>0, the overall error analysis must be revisited to determine whether the dimension-independent rate survives after h is fixed by the discretization tolerance.

minor comments (2)

[§2] The notation for the Gevrey class and the sequence of multi-indices should be introduced once and used consistently; occasional switches between bold-face and subscript notation for the parameter vector obscure the counting of derivatives.
[Numerical experiments] Figure 2 caption should state the precise value of the mesh size h and the truncation dimension used in the QMC experiments so that the observed rates can be directly compared with the theoretical constants.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for the careful and constructive report. The comments correctly emphasize the necessity of h-independent constants to preserve the dimension-independent QMC rates. We address each major comment below and indicate the revisions planned.

read point-by-point responses

Referee: [§3] §3 (parametric regularity assumptions): The propagation of Gevrey regularity from the random field to the discrete flux and gradient requires that all stability and approximation constants in the discrete balance law remain independent of the mesh size h. The manuscript does not explicitly bound the h-dependence arising from inverse estimates in the flux space; if such factors appear, the mixed partial derivatives grow with negative powers of h, which would make the QMC weights h-dependent and destroy both the dimension-independent rate and the faster-than-MC claim for any fixed discretization.

Authors: We agree that explicit control over h-dependence is required. The assumptions in §3 are stated so that stability and approximation constants (including those from inverse estimates) are independent of h under standard shape-regularity conditions on the mesh family. To address the referee's concern directly, we will insert a short lemma after Assumption 3.2 that explicitly bounds the constants arising from inverse inequalities in the flux space, confirming they remain uniform in h for the finite-element families considered. This addition will make the h-independence of the Gevrey constants fully transparent without altering the existing analysis. revision: yes
Referee: [§4.3] §4.3 (verification for HDG schemes): The claim that the discrete flux inherits Gevrey-type bounds uniform in h relies on the hybridizable formulation preserving the required Lipschitz and differentiability properties of the constitutive map without h-dependent constants. No explicit estimate is given showing that the stabilization parameter and the trace inequalities used in the HDG flux reconstruction remain analytic with h-independent radii of analyticity.

Authors: The referee rightly notes that the HDG verification would benefit from more explicit tracking of constants. In the current §4.3 the stabilization parameter is chosen as a fixed multiple of the local diffusivity bound (independent of h) and the trace inequalities are invoked under the standard assumption of shape-regular meshes. In the revision we will add a dedicated paragraph that recalls the h-independent analyticity radius for the constitutive map and cites the uniform trace-inverse estimates that hold for the HDG flux reconstruction on shape-regular triangulations. This will confirm that the Gevrey constants remain uniform in h. revision: yes
Referee: [§5] §5 (QMC error bounds): The derivation of the QMC convergence rate assumes that the parametric regularity constants are independent of the discretization parameter h once the mesh is chosen to balance discretization and sampling errors. If the constants actually deteriorate as O(h^{-k}) for some k>0, the overall error analysis must be revisited to determine whether the dimension-independent rate survives after h is fixed by the discretization tolerance.

Authors: Under the assumptions and verifications already given in §§3–4, the parametric regularity constants are independent of h. With h fixed by a prescribed discretization tolerance, the QMC sampling error therefore retains its dimension-independent, faster-than-Monte-Carlo rate. In the revised manuscript we will add a clarifying remark at the end of §5 that explicitly states this independence and sketches the combined discretization-plus-sampling error bound, thereby confirming that the claimed rates survive the balancing of h. revision: partial

Circularity Check

0 steps flagged

No significant circularity; bounds derived from explicit regularity assumptions and verified discretization properties.

full rationale

The derivation begins from the modeling of the random parameter as a Gevrey-regular function of countably many i.i.d. uniforms or normals, propagates this regularity through the balance law and finite-element discretization operators under stated stability and approximation assumptions that are independent of the target QMC rate, and obtains the dimension-independent convergence from the resulting mixed-derivative bounds together with continuity of the QoI. These steps are self-contained: the assumptions are verified directly for conforming, mixed, and HDG schemes without reduction to fitted parameters, self-referential definitions, or load-bearing self-citations. The central claim therefore rests on externally stated hypotheses rather than on any quantity defined in terms of the conclusion.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The central claim rests on standard functional-analysis and probability assumptions for random fields in PDEs together with discretization-specific stability and approximation properties; no new free parameters or invented entities are introduced.

axioms (2)

domain assumption The random parameter is a Gevrey-regular random field expressible as a function of countably infinite sequences of independent random variables that may be uniformly or normally distributed.
This modeling choice is invoked to obtain the parametric regularity that underpins the QMC error bounds.
domain assumption The discretizations satisfy a series of assumptions that encompass a broad class of finite element methods for flux vector fields satisfying a balance law.
These assumptions are introduced to allow the QMC analysis to apply to conforming, mixed, and HDG schemes.

pith-pipeline@v0.9.0 · 5744 in / 1503 out tokens · 36499 ms · 2026-05-18T01:11:21.460679+00:00 · methodology

discussion (0)

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

IndisputableMonolith/Cost/FunctionalEquation.lean washburn_uniqueness_aczel unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

It is shown that the QMC method achieves a dimension-independent, faster-than-Monte Carlo convergence rate if the quantity of interest depends continuously on the primal variable, its flux, or its gradient.
IndisputableMonolith/Foundation/RealityFromDistinction.lean reality_from_one_distinction unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

The random parameter is modeled here as a Gevrey-regular random field... (A5) Coefficient of class Gevrey-σ

What do these tags mean?

matches: The paper's claim is directly supported by a theorem in the formal canon.
supports: The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends: The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses: The paper appears to rely on the theorem as machinery.
contradicts: The paper's claim conflicts with a theorem or certificate in the canon.
unclear: Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

36 extracted references · 36 canonical work pages

[1]

S. C. Brenner,Poincar´ e–Friedrichs inequalities for piecewise H 1 functions, SIAM J. Numer. Anal., 41 (2003), pp. 306–324, https://doi.org/10.1137/S0036142902401311. [2]S. B. Chae,Holomorphy and Calculus in Normed Spaces, CRC Press, 1985

work page doi:10.1137/s0036142902401311 2003
[2]

Chernov and T

A. Chernov and T. L ˆe,Analytic and Gevrey class regularity for parametric elliptic eigenvalue prob- lems and applications, SIAM J. Numer. Anal., 62 (2024), pp. 1874–1900, https://doi.org/10.1137/ 23M1596296

work page 2024
[3]

Chernov and T

A. Chernov and T. L ˆe,Analytic and Gevrey class regularity for parametric semilinear reaction-diffusion problems and applications in uncertainty quantification, Comput. Math. Appl., 164 (2024), pp. 116–130, https://doi.org/10.1016/j.camwa.2024.04.007

work page doi:10.1016/j.camwa.2024.04.007 2024
[4]

Cockburn, J

B. Cockburn, J. Gopalakrishnan, and R. Lazarov,Unified hybridization of discontinuous Galerkin, mixed, and continuous Galerkin methods for second order elliptic problems, SIAM J. Numer. Anal., 47 (2009), pp. 1319–1365, https://doi.org/10.1137/070706616

work page doi:10.1137/070706616 2009
[5]

Cohen and R

A. Cohen and R. DeVore,Approximation of high-dimensional parametric PDEs, Acta Numer., 24 (2015), pp. 1–159, https://doi.org/10.1017/S0962492915000033

work page doi:10.1017/s0962492915000033 2015
[6]

Cohen, R

A. Cohen, R. DeVore, and C. Schwab,Convergence rates of best N-term Galerkin approximations for a class of elliptic sPDEs, Found. Comput. Math., 10 (2010), pp. 615–646, https://doi.org/10.1007/ s10208-010-9072-2

work page 2010
[7]

Cools, F

R. Cools, F. Y. Kuo, and D. Nuyens,Constructing embedded lattice rules for multivariate integration, SIAM J. Sci. Comput., 28 (2006), pp. 2162–2188, https://doi.org/10.1137/06065074X

work page doi:10.1137/06065074x 2006
[8]

D. A. Di Pietro, Z. Dong, G. Kanschat, P. Matalon, and A. Rupp,Homogeneous multigrid for hybrid discretizations: application to HHO methods, Numer. Methods Partial Differential Equations, 41 (2025), p. e70023, https://doi.org/10.1002/num.70023

work page doi:10.1002/num.70023 2025
[9]

J. Dick, P. Kritzer, and F. Pillichshammer,Lattice Rules: Numerical Integration, Approximation, and Discrepancy, Springer Series in Computational Mathematics, Springer International Publishing, 2022, https://doi.org/10.1007/978-3-031-09951-9. 28 V. KAARNIOJA, A. RUPP, J. GOPALAKRISHNAN

work page doi:10.1007/978-3-031-09951-9 2022
[10]

J. Dick, F. Y. Kuo, Q. T. Le Gia, D. Nuyens, and C. Schwab,Higher order QMC Petrov–Galerkin discretization for affine parametric operator equations with random field inputs, SIAM J. Numer. Anal., 52 (2014), pp. 2676–2702, https://doi.org/10.1137/130943984

work page doi:10.1137/130943984 2014
[11]

J. Dick, F. Y. Kuo, and I. H. Sloan,High-dimensional integration: the quasi-Monte Carlo way, Acta Numer., 22 (2013), pp. 133–288, https://doi.org/10.1017/S0962492913000044

work page doi:10.1017/s0962492913000044 2013
[12]

J. Dick, Q. T. Le Gia, and C. Schwab,Higher order quasi-Monte Carlo integration for holomorphic, parametric operator equations, SIAM/ASA J. Uncertain. Quantif., 4 (2016), pp. 48–79, https://doi. org/10.1137/140985913

work page doi:10.1137/140985913 2016
[13]

O. G. Ernst and B. Sprungk,Stochastic collocation for elliptic PDEs with random data: the lognormal case, in Sparse Grids and Applications - Munich 2012, J. Garcke and D. Pfl¨ uger, eds., Cham, 2014, Springer International Publishing, pp. 29–53, https://doi.org/10.1007/978-3-319-04537-5 2

work page doi:10.1007/978-3-319-04537-5 2012
[14]

R. N. Gantner, L. Herrmann, and C. Schwab,Multilevel QMC with product weights for affine- parametric, elliptic PDEs, in Contemporary Computational Mathematics - A Celebration of the 80th Birthday of Ian Sloan, J. Dick, F. Y. Kuo, and H. Wo´ zniakowski, eds., Springer International Publishing, 2018, pp. 373–405, https://doi.org/10.1007/978-3-319-72456-0 18...

work page doi:10.1007/978-3-319-72456-0 2018
[15]

I. G. Graham, F. Y. Kuo, J. A. Nichols, R. Scheichl, C. Schwab, and I. H. Sloan,Quasi-Monte Carlo finite element methods for elliptic PDEs with lognormal random coefficients, Numer. Math., 131 (2015), pp. 329–368, https://doi.org/10.1007/s00211-014-0689-y

work page doi:10.1007/s00211-014-0689-y 2015
[16]

I. G. Graham, F. Y. Kuo, D. Nuyens, R. Scheichl, and I. H. Sloan,Quasi-Monte Carlo methods for elliptic PDEs with random coefficients and applications, J. Comput. Phys., 230 (2011), pp. 3668–3694, https://doi.org/10.1016/j.jcp.2011.01.023

work page doi:10.1016/j.jcp.2011.01.023 2011
[17]

I. G. Graham, F. Y. Kuo, D. Nuyens, R. Scheichl, and I. H. Sloan,Circulant embedding with QMC: analysis for elliptic PDE with lognormal coefficients, Numer. Math., 140 (2018), pp. 479–511, https://doi.org/10.1007/s00211-018-0968-0

work page doi:10.1007/s00211-018-0968-0 2018
[18]

I. G. Graham, R. Scheichl, and E. Ullmann,Mixed finite element analysis of lognormal diffusion and multilevel Monte Carlo methods, Stoch. Partial Differ. Equ. Anal. Comput., 4 (2015), pp. 41–75, https://doi.org/10.1007/s40072-015-0051-0

work page doi:10.1007/s40072-015-0051-0 2015
[19]

P. A. Guth and V. Kaarnioja,Quasi-Monte Carlo for partial differential equations with generalized Gaussian input uncertainty, SIAM J. Numer. Anal., 63 (2025), pp. 1666–1690, https://doi.org/10. 1137/24M1708164

work page 2025
[20]

Harbrecht, M

H. Harbrecht, M. Schmidlin, and C. Schwab,The Gevrey class implicit mapping theorem with application to UQ of semilinear elliptic PDEs, Math. Models Methods Appl. Sci., 34 (2024), pp. 881– 917, https://doi.org/10.1142/S0218202524500179. [23]G. Hardy, J. Littlewood, and G. P ´olya,Inequalities, Cambridge University Press, 1934

work page doi:10.1142/s0218202524500179 2024
[21]

Herrmann and C

L. Herrmann and C. Schwab,QMC integration for lognormal-parametric elliptic PDEs: local supports and product weights, Numer. Math., 141 (2019), pp. 63–102, https://doi.org/10.1007/ s00211-018-0991-1

work page 2019
[22]

Kaarnioja, Y

V. Kaarnioja, Y. Kazashi, F. Y. Kuo, F. Nobile, and I. H. Sloan,Fast approximation by periodic kernel-based lattice-point interpolation with application in uncertainty quantification, Numer. Math., 150 (2022), pp. 33–77, https://doi.org/10.1007/s00211-021-01242-3

work page doi:10.1007/s00211-021-01242-3 2022
[23]

Kaarnioja and A

V. Kaarnioja and A. Rupp,Quasi-Monte Carlo and discontinuous Galerkin, Electron. Trans. Numer. Anal., 60 (2024), pp. 589–617, https://doi.org/10.1553/etna vol60s589

work page doi:10.1553/etna 2024
[24]

Sufficient conditions for QMC analysis of finite elements for parametric differential equations

V. Kaarnioja, A. Rupp, and J. Gopalakrishnan,Code to generate the numerical results in “Sufficient conditions for QMC analysis of finite elements for parametric differential equations”, 2025, https: //doi.org/10.5281/zenodo.17273609. Zenodo:17273609. [28]F. Y. Kuo,Lattice rule generating vectors. https://web.maths.unsw.edu.au/ ∼fkuo/lattice/index.html

work page doi:10.5281/zenodo.17273609 2025
[25]

F. Y. Kuo and D. Nuyens,Application of quasi-Monte Carlo methods to elliptic PDEs with random diffusion coefficients: a survey of analysis and implementation, Found. Comput. Math., 16 (2016), pp. 1631–1696, https://doi.org/10.1007/s10208-016-9329-5

work page doi:10.1007/s10208-016-9329-5 2016
[26]

F. Y. Kuo, R. Scheichl, C. Schwab, I. H. Sloan, and E. Ullmann,Multilevel Quasi-Monte Carlo methods for lognormal diffusion problems, Math. Comp., 86 (2017), pp. 2827–2860

work page 2017
[27]

F. Y. Kuo, C. Schwab, and I. H. Sloan,Quasi-Monte Carlo finite element methods for a class of elliptic FEM FOR PARAMETRIC PDES 29 partial differential equations with random coefficients, SIAM J. Numer. Anal., 50 (2012), pp. 3351–3374, https://doi.org/10.1137/110845537

work page doi:10.1137/110845537 2012
[28]

F. Y. Kuo, C. Schwab, and I. H. Sloan,Multi-level quasi-Monte Carlo finite element methods for a class of elliptic PDEs with random coefficients, Found. Comput. Math., 15 (2015), pp. 411–449, https://doi.org/10.1007/s10208-014-9237-5. [33]P. Lu, R. Maier, and A. Rupp,A localized orthogonal decomposition strategy for hybrid discontinuous Galerkin methods, ...

work page doi:10.1007/s10208-014-9237-5 2015
[29]

P. Lu, A. Rupp, and G. Kanschat,Analysis of injection operators in geometric multigrid solvers for HDG methods, SIAM Journal on Numerical Analysis, 60 (2022), pp. 2293–2317, https://doi.org/10. 1137/21M1400110

work page 2022
[30]

J. A. Nichols and F. Y. Kuo,Fast CBC construction of randomly shifted lattice rule achieving O(n−1+δ) convergence for unbounded integrals over Rs in weighted spaces with POD weights, J. Complexity, 30 (2014), pp. 444–468, https://doi.org/10.1016/j.jco.2014.02.004

work page doi:10.1016/j.jco.2014.02.004 2014
[31]

Nuyens and R

D. Nuyens and R. Cools,Fast algorithms for component-by-component construction of rank-1 lattice rules in shift-invariant reproducing kernel Hilbert spaces, Math. Comp., 75 (2006), pp. 903–920, https://doi.org/10.1090/S0025-5718-06-01785-6

work page doi:10.1090/s0025-5718-06-01785-6 2006
[32]

M. Y. Qi and L. Rodino,General Theory of Partial Differential Equations and Microlocal Analysis, Chapman & Hall/CRC Research Notes in Mathematics Series, Taylor & Francis, 1996

work page 1996
[33]

Sch¨oberl,C++ 11 implementation of finite elements in NGSolve, Institute for Analysis and Scientific Computing, Vienna University of Technology, 30 (2014)

J. Sch¨oberl,C++ 11 implementation of finite elements in NGSolve, Institute for Analysis and Scientific Computing, Vienna University of Technology, 30 (2014)

work page 2014
[34]

Schwab,QMC Galerkin discretization of parametric operator equations, in Monte Carlo and Quasi- Monte Carlo Methods 2012, J

C. Schwab,QMC Galerkin discretization of parametric operator equations, in Monte Carlo and Quasi- Monte Carlo Methods 2012, J. Dick, F. Y. Kuo, G. W. Peters, and I. H. Sloan, eds., Springer Berlin Heidelberg, 2013, pp. 613–629, https://doi.org/10.1007/978-3-642-41095-6 32

work page doi:10.1007/978-3-642-41095-6 2012
[35]

Schwab and C

C. Schwab and C. J. Gittelson,Sparse tensor discretizations of high-dimensional parametric and stochastic PDEs, Acta Numer., 20 (2011), pp. 291–467, https://doi.org/10.1017/S0962492911000055

work page doi:10.1017/s0962492911000055 2011
[36]

Schwab and R

C. Schwab and R. A. Todor,Karhunen–Lo` eve approximation of random fields by generalized fast multipole methods, J. Comput. Phys., 217 (2006), pp. 100–122, https://doi.org/10.1016/j.jcp.2006.01. 048. Appendix A. More details. A.1. Proof of (2.2).Using a translational change of variable in the definition (2.1) of I φ s , we have, for any fixedk, I φ s = Z ...

work page doi:10.1016/j.jcp.2006.01 2006

[1] [1]

S. C. Brenner,Poincar´ e–Friedrichs inequalities for piecewise H 1 functions, SIAM J. Numer. Anal., 41 (2003), pp. 306–324, https://doi.org/10.1137/S0036142902401311. [2]S. B. Chae,Holomorphy and Calculus in Normed Spaces, CRC Press, 1985

work page doi:10.1137/s0036142902401311 2003

[2] [2]

Chernov and T

A. Chernov and T. L ˆe,Analytic and Gevrey class regularity for parametric elliptic eigenvalue prob- lems and applications, SIAM J. Numer. Anal., 62 (2024), pp. 1874–1900, https://doi.org/10.1137/ 23M1596296

work page 2024

[3] [3]

Chernov and T

A. Chernov and T. L ˆe,Analytic and Gevrey class regularity for parametric semilinear reaction-diffusion problems and applications in uncertainty quantification, Comput. Math. Appl., 164 (2024), pp. 116–130, https://doi.org/10.1016/j.camwa.2024.04.007

work page doi:10.1016/j.camwa.2024.04.007 2024

[4] [4]

Cockburn, J

B. Cockburn, J. Gopalakrishnan, and R. Lazarov,Unified hybridization of discontinuous Galerkin, mixed, and continuous Galerkin methods for second order elliptic problems, SIAM J. Numer. Anal., 47 (2009), pp. 1319–1365, https://doi.org/10.1137/070706616

work page doi:10.1137/070706616 2009

[5] [5]

Cohen and R

A. Cohen and R. DeVore,Approximation of high-dimensional parametric PDEs, Acta Numer., 24 (2015), pp. 1–159, https://doi.org/10.1017/S0962492915000033

work page doi:10.1017/s0962492915000033 2015

[6] [6]

Cohen, R

A. Cohen, R. DeVore, and C. Schwab,Convergence rates of best N-term Galerkin approximations for a class of elliptic sPDEs, Found. Comput. Math., 10 (2010), pp. 615–646, https://doi.org/10.1007/ s10208-010-9072-2

work page 2010

[7] [7]

Cools, F

R. Cools, F. Y. Kuo, and D. Nuyens,Constructing embedded lattice rules for multivariate integration, SIAM J. Sci. Comput., 28 (2006), pp. 2162–2188, https://doi.org/10.1137/06065074X

work page doi:10.1137/06065074x 2006

[8] [8]

D. A. Di Pietro, Z. Dong, G. Kanschat, P. Matalon, and A. Rupp,Homogeneous multigrid for hybrid discretizations: application to HHO methods, Numer. Methods Partial Differential Equations, 41 (2025), p. e70023, https://doi.org/10.1002/num.70023

work page doi:10.1002/num.70023 2025

[9] [9]

J. Dick, P. Kritzer, and F. Pillichshammer,Lattice Rules: Numerical Integration, Approximation, and Discrepancy, Springer Series in Computational Mathematics, Springer International Publishing, 2022, https://doi.org/10.1007/978-3-031-09951-9. 28 V. KAARNIOJA, A. RUPP, J. GOPALAKRISHNAN

work page doi:10.1007/978-3-031-09951-9 2022

[10] [10]

J. Dick, F. Y. Kuo, Q. T. Le Gia, D. Nuyens, and C. Schwab,Higher order QMC Petrov–Galerkin discretization for affine parametric operator equations with random field inputs, SIAM J. Numer. Anal., 52 (2014), pp. 2676–2702, https://doi.org/10.1137/130943984

work page doi:10.1137/130943984 2014

[11] [11]

J. Dick, F. Y. Kuo, and I. H. Sloan,High-dimensional integration: the quasi-Monte Carlo way, Acta Numer., 22 (2013), pp. 133–288, https://doi.org/10.1017/S0962492913000044

work page doi:10.1017/s0962492913000044 2013

[12] [12]

J. Dick, Q. T. Le Gia, and C. Schwab,Higher order quasi-Monte Carlo integration for holomorphic, parametric operator equations, SIAM/ASA J. Uncertain. Quantif., 4 (2016), pp. 48–79, https://doi. org/10.1137/140985913

work page doi:10.1137/140985913 2016

[13] [13]

O. G. Ernst and B. Sprungk,Stochastic collocation for elliptic PDEs with random data: the lognormal case, in Sparse Grids and Applications - Munich 2012, J. Garcke and D. Pfl¨ uger, eds., Cham, 2014, Springer International Publishing, pp. 29–53, https://doi.org/10.1007/978-3-319-04537-5 2

work page doi:10.1007/978-3-319-04537-5 2012

[14] [14]

R. N. Gantner, L. Herrmann, and C. Schwab,Multilevel QMC with product weights for affine- parametric, elliptic PDEs, in Contemporary Computational Mathematics - A Celebration of the 80th Birthday of Ian Sloan, J. Dick, F. Y. Kuo, and H. Wo´ zniakowski, eds., Springer International Publishing, 2018, pp. 373–405, https://doi.org/10.1007/978-3-319-72456-0 18...

work page doi:10.1007/978-3-319-72456-0 2018

[15] [15]

I. G. Graham, F. Y. Kuo, J. A. Nichols, R. Scheichl, C. Schwab, and I. H. Sloan,Quasi-Monte Carlo finite element methods for elliptic PDEs with lognormal random coefficients, Numer. Math., 131 (2015), pp. 329–368, https://doi.org/10.1007/s00211-014-0689-y

work page doi:10.1007/s00211-014-0689-y 2015

[16] [16]

I. G. Graham, F. Y. Kuo, D. Nuyens, R. Scheichl, and I. H. Sloan,Quasi-Monte Carlo methods for elliptic PDEs with random coefficients and applications, J. Comput. Phys., 230 (2011), pp. 3668–3694, https://doi.org/10.1016/j.jcp.2011.01.023

work page doi:10.1016/j.jcp.2011.01.023 2011

[17] [17]

I. G. Graham, F. Y. Kuo, D. Nuyens, R. Scheichl, and I. H. Sloan,Circulant embedding with QMC: analysis for elliptic PDE with lognormal coefficients, Numer. Math., 140 (2018), pp. 479–511, https://doi.org/10.1007/s00211-018-0968-0

work page doi:10.1007/s00211-018-0968-0 2018

[18] [18]

I. G. Graham, R. Scheichl, and E. Ullmann,Mixed finite element analysis of lognormal diffusion and multilevel Monte Carlo methods, Stoch. Partial Differ. Equ. Anal. Comput., 4 (2015), pp. 41–75, https://doi.org/10.1007/s40072-015-0051-0

work page doi:10.1007/s40072-015-0051-0 2015

[19] [19]

P. A. Guth and V. Kaarnioja,Quasi-Monte Carlo for partial differential equations with generalized Gaussian input uncertainty, SIAM J. Numer. Anal., 63 (2025), pp. 1666–1690, https://doi.org/10. 1137/24M1708164

work page 2025

[20] [20]

Harbrecht, M

H. Harbrecht, M. Schmidlin, and C. Schwab,The Gevrey class implicit mapping theorem with application to UQ of semilinear elliptic PDEs, Math. Models Methods Appl. Sci., 34 (2024), pp. 881– 917, https://doi.org/10.1142/S0218202524500179. [23]G. Hardy, J. Littlewood, and G. P ´olya,Inequalities, Cambridge University Press, 1934

work page doi:10.1142/s0218202524500179 2024

[21] [21]

Herrmann and C

L. Herrmann and C. Schwab,QMC integration for lognormal-parametric elliptic PDEs: local supports and product weights, Numer. Math., 141 (2019), pp. 63–102, https://doi.org/10.1007/ s00211-018-0991-1

work page 2019

[22] [22]

Kaarnioja, Y

V. Kaarnioja, Y. Kazashi, F. Y. Kuo, F. Nobile, and I. H. Sloan,Fast approximation by periodic kernel-based lattice-point interpolation with application in uncertainty quantification, Numer. Math., 150 (2022), pp. 33–77, https://doi.org/10.1007/s00211-021-01242-3

work page doi:10.1007/s00211-021-01242-3 2022

[23] [23]

Kaarnioja and A

V. Kaarnioja and A. Rupp,Quasi-Monte Carlo and discontinuous Galerkin, Electron. Trans. Numer. Anal., 60 (2024), pp. 589–617, https://doi.org/10.1553/etna vol60s589

work page doi:10.1553/etna 2024

[24] [24]

Sufficient conditions for QMC analysis of finite elements for parametric differential equations

V. Kaarnioja, A. Rupp, and J. Gopalakrishnan,Code to generate the numerical results in “Sufficient conditions for QMC analysis of finite elements for parametric differential equations”, 2025, https: //doi.org/10.5281/zenodo.17273609. Zenodo:17273609. [28]F. Y. Kuo,Lattice rule generating vectors. https://web.maths.unsw.edu.au/ ∼fkuo/lattice/index.html

work page doi:10.5281/zenodo.17273609 2025

[25] [25]

F. Y. Kuo and D. Nuyens,Application of quasi-Monte Carlo methods to elliptic PDEs with random diffusion coefficients: a survey of analysis and implementation, Found. Comput. Math., 16 (2016), pp. 1631–1696, https://doi.org/10.1007/s10208-016-9329-5

work page doi:10.1007/s10208-016-9329-5 2016

[26] [26]

F. Y. Kuo, R. Scheichl, C. Schwab, I. H. Sloan, and E. Ullmann,Multilevel Quasi-Monte Carlo methods for lognormal diffusion problems, Math. Comp., 86 (2017), pp. 2827–2860

work page 2017

[27] [27]

F. Y. Kuo, C. Schwab, and I. H. Sloan,Quasi-Monte Carlo finite element methods for a class of elliptic FEM FOR PARAMETRIC PDES 29 partial differential equations with random coefficients, SIAM J. Numer. Anal., 50 (2012), pp. 3351–3374, https://doi.org/10.1137/110845537

work page doi:10.1137/110845537 2012

[28] [28]

F. Y. Kuo, C. Schwab, and I. H. Sloan,Multi-level quasi-Monte Carlo finite element methods for a class of elliptic PDEs with random coefficients, Found. Comput. Math., 15 (2015), pp. 411–449, https://doi.org/10.1007/s10208-014-9237-5. [33]P. Lu, R. Maier, and A. Rupp,A localized orthogonal decomposition strategy for hybrid discontinuous Galerkin methods, ...

work page doi:10.1007/s10208-014-9237-5 2015

[29] [29]

P. Lu, A. Rupp, and G. Kanschat,Analysis of injection operators in geometric multigrid solvers for HDG methods, SIAM Journal on Numerical Analysis, 60 (2022), pp. 2293–2317, https://doi.org/10. 1137/21M1400110

work page 2022

[30] [30]

J. A. Nichols and F. Y. Kuo,Fast CBC construction of randomly shifted lattice rule achieving O(n−1+δ) convergence for unbounded integrals over Rs in weighted spaces with POD weights, J. Complexity, 30 (2014), pp. 444–468, https://doi.org/10.1016/j.jco.2014.02.004

work page doi:10.1016/j.jco.2014.02.004 2014

[31] [31]

Nuyens and R

D. Nuyens and R. Cools,Fast algorithms for component-by-component construction of rank-1 lattice rules in shift-invariant reproducing kernel Hilbert spaces, Math. Comp., 75 (2006), pp. 903–920, https://doi.org/10.1090/S0025-5718-06-01785-6

work page doi:10.1090/s0025-5718-06-01785-6 2006

[32] [32]

M. Y. Qi and L. Rodino,General Theory of Partial Differential Equations and Microlocal Analysis, Chapman & Hall/CRC Research Notes in Mathematics Series, Taylor & Francis, 1996

work page 1996

[33] [33]

Sch¨oberl,C++ 11 implementation of finite elements in NGSolve, Institute for Analysis and Scientific Computing, Vienna University of Technology, 30 (2014)

J. Sch¨oberl,C++ 11 implementation of finite elements in NGSolve, Institute for Analysis and Scientific Computing, Vienna University of Technology, 30 (2014)

work page 2014

[34] [34]

Schwab,QMC Galerkin discretization of parametric operator equations, in Monte Carlo and Quasi- Monte Carlo Methods 2012, J

C. Schwab,QMC Galerkin discretization of parametric operator equations, in Monte Carlo and Quasi- Monte Carlo Methods 2012, J. Dick, F. Y. Kuo, G. W. Peters, and I. H. Sloan, eds., Springer Berlin Heidelberg, 2013, pp. 613–629, https://doi.org/10.1007/978-3-642-41095-6 32

work page doi:10.1007/978-3-642-41095-6 2012

[35] [35]

Schwab and C

C. Schwab and C. J. Gittelson,Sparse tensor discretizations of high-dimensional parametric and stochastic PDEs, Acta Numer., 20 (2011), pp. 291–467, https://doi.org/10.1017/S0962492911000055

work page doi:10.1017/s0962492911000055 2011

[36] [36]

Schwab and R

C. Schwab and R. A. Todor,Karhunen–Lo` eve approximation of random fields by generalized fast multipole methods, J. Comput. Phys., 217 (2006), pp. 100–122, https://doi.org/10.1016/j.jcp.2006.01. 048. Appendix A. More details. A.1. Proof of (2.2).Using a translational change of variable in the definition (2.1) of I φ s , we have, for any fixedk, I φ s = Z ...

work page doi:10.1016/j.jcp.2006.01 2006