On Equivalence of Anchored and ANOVA Spaces; Lower Bounds

We provide lower bounds for the norms of embeddings between $\boldsymbol{\gamma}$-weighted Anchored and ANOVA spaces of $s$-variate functions with mixed partial derivatives of order one bounded in $L_p$ norm ($p\in[1,\infty]$). In particular we show that the norms behave polynomially in $s$ for Finite Order Weights and Finite Diameter Weights if $p>1$, and increase faster than any polynomial in $s$ for Product Order-Dependent Weights and any $p$

On efficient weighted integration via a change of variables

In this paper, we study the approximation of $d$-dimensional $\rho$-weighted integrals over unbounded domains $\mathbb{R}_+^d$ or $\mathbb{R}^d$ using a special change of variables, so that quasi-Monte Carlo (QMC) or sparse grid rules can be applied to the transformed integrands over the unit cube. We consider a class of integrands with bounded $L_p$ norm of mixed partial derivatives of first order, where $p\in[1,+\infty].$ The main results give sufficient conditions on the change of variables $\nu$ which guarantee that the transformed integrand belongs to the standard Sobolev space of functions over the unit cube with mixed smoothness of order one. These conditions depend on $\rho$ and $p$. The proposed change of variables is in general different than the standard change based on the inverse of the cumulative distribution function. We stress that the standard change of variables leads to integrands over a cube; however, those integrands have singularities which make the application of QMC and sparse grids ineffective. Our conclusions are supported by numerical experiments

Very Low Truncation Dimension for High Dimensional Integration Under Modest Error Demand

We consider the problem of numerical integration for weighted anchored and ANOVA Sobolev spaces of $s$-variate functions. Here $s$ is large including $s=\infty$. Under the assumption of sufficiently fast decaying weights, we prove in a constructive way that such integrals can be approximated by quadratures for functions $f_k$ with only $k$ variables, where $k=k(\varepsilon)$ depends solely on the error demand $\varepsilon$ and is surprisingly small when $s$ is sufficiently large relative to $\varepsilon$. This holds, in particular, for $s=\infty$ and arbitrary $\varepsilon$ since then $k(\varepsilon)<\infty$ for all $\varepsilon$. Moreover $k(\varepsilon)$ does not depend on the function being integrated, i.e., is the same for all functions from the unit ball of the space

Effective dimension of some weighted pre-Sobolev spaces with dominating mixed partial derivatives

This paper considers two notions of effective dimension for quadrature in weighted pre-Sobolev spaces with dominating mixed partial derivatives. We begin by finding a ball in those spaces just barely large enough to contain a function with unit variance. If no function in that ball has more than $\varepsilon$ of its variance from ANOVA components involving interactions of order $s$ or more, then the space has effective dimension at most $s$ in the superposition sense. A similar truncation sense notion replaces the cardinality of the ANOVA component by the largest index it contains. Some Poincar\'e type inequalities are used to bound variance components by multiples of these space's squared norm and those in turn provide bounds on effective dimension. Very low effective dimension in the superposition sense holds for some spaces defined by product weights in which quadrature is strongly tractable. The superposition dimension is $O( \log(1/\varepsilon)/\log(\log(1/\varepsilon)))$ just like the superposition dimension used in the multidimensional decomposition method. Surprisingly, even spaces where all subset weights are equal, regardless of their cardinality or included indices, have low superposition dimension in this sense. This paper does not require periodicity of the integrands

Efficient implementations of the Multivariate Decomposition Method for approximating infinite-variate integrals

In this paper we focus on efficient implementations of the Multivariate Decomposition Method (MDM) for approximating integrals of $\infty$-variate functions. Such $\infty$-variate integrals occur for example as expectations in uncertainty quantification. Starting with the anchored decomposition $f = \sum_{\mathfrak{u}\subset\mathbb{N}} f_\mathfrak{u}$, where the sum is over all finite subsets of $\mathbb{N}$ and each $f_\mathfrak{u}$ depends only on the variables $x_j$ with $j\in\mathfrak{u}$, our MDM algorithm approximates the integral of $f$ by first truncating the sum to some `active set' and then approximating the integral of the remaining functions $f_\mathfrak{u}$ term-by-term using Smolyak or (randomized) quasi-Monte Carlo (QMC) quadratures. The anchored decomposition allows us to compute $f_\mathfrak{u}$ explicitly by function evaluations of $f$. Given the specification of the active set and theoretically derived parameters of the quadrature rules, we exploit structures in both the formula for computing $f_\mathfrak{u}$ and the quadrature rules to develop computationally efficient strategies to implement the MDM in various scenarios. In particular, we avoid repeated function evaluations at the same point. We provide numerical results for a test function to demonstrate the effectiveness of the algorithm