BZPEER

preprint2026arXiv

Sparse-grid polynomial interpolation approximation and integration for parametric and stochastic elliptic PDEs with lognormal inputs

By combining a certain approximation property in the spatial domain, and weighted $\ell_2$ -summability of the Hermite polynomial expansion coefficients in the parametric domain obtained in [M. Bachmayr, A. Cohen, R. DeVore and G. Migliorati, ESAIM Math. Model. Numer. Anal. $\bf 51$ (2017), 341-363] and [M. Bachmayr, A. Cohen, D. Dũng and C. Schwab, SIAM J. Numer. Anal. $\bf 55$ (2017), 2151-2186], we investigate linear non-adaptive methods of fully discrete polynomial interpolation approximation as well as fully discrete weighted quadrature methods of integration for parametric and stochastic elliptic PDEs with lognormal inputs. We explicitly construct such methods and prove corresponding convergence rates in $n$ of the approximations by them, where $n$ is a number characterizing computation complexity. The linear non-adaptive methods of fully discrete polynomial interpolation approximation are sparse-grid collocation methods. Moreover, they generate in a natural way discrete weighted quadrature formulas for integration of the solution to parametric and stochastic elliptic PDEs and its linear functionals, and the error of the corresponding integration can be estimated via the error in the Bochner space $L_1({\mathbb R}^\infty,V,γ)$ norm of the generating methods where $γ$ is the Gaussian probability measure on ${\mathbb R}^\infty$ and $V$ is the energy space. We also briefly consider similar problems for parametric and stochastic elliptic PDEs with affine inputs, and by-product problems of non-fully discrete polynomial interpolation approximation and integration. In particular, the convergence rate of non-fully discrete obtained in this paper improves the known one.

preprint2026arXiv

Weighted approximate sampling recovery and integration based on B-spline interpolation and quasi-interpolation

We propose novel methods for approximate sampling recovery and integration of functions in the Freud-weighted Sobolev space $W^r_{p,w}(\mathbb{R})$ . The approximation error of sampling recovery is measured in the norm of the Freud-weighted Lebesgue space $L_{q,w}(\mathbb{R})$ . Namely, we construct equidistant compact-supported B-spline quasi-interpolation and interpolation sampling algorithms $Q_{ρ,m}$ and $P_{ρ,m}$ which are asymptotically optimal in terms of the sampling $n$ -widths $\varrho_n(\boldsymbol{W}^r_{p,w}(\mathbb{R}), L_{q,w}(\mathbb{R}))$ for every pair $p,q \in [1,\infty]$ , and prove the right convergence rate of these sampling $n$ -widths, where $\boldsymbol{W}^r_{p,w}(\mathbb{R})$ denotes the unit ball in $W^r_{p,w}(\mathbb{R})$ . The algorithms $Q_{ρ,m}$ and $P_{ρ,m}$ are based on truncated scaled B-spline quasi-interpolation and interpolation, respectively. We also prove the asymptotical optimality and right convergence rate of the equidistant quadratures generated from $Q_{ρ,m}$ and $P_{ρ,m}$ , for Freud-weighted numerical integration of functions in $W^r_{p,w}(\mathbb{R})$ .

preprint2021arXiv

Computation complexity of deep ReLU neural networks in high-dimensional approximation

The purpose of the present paper is to study the computation complexity of deep ReLU neural networks to approximate functions in Hölder-Nikol'skii spaces of mixed smoothness $H_\infty^α(\mathbb{I}^d)$ on the unit cube $\mathbb{I}^d:=[0,1]^d$ . In this context, for any function $f\in H_\infty^α(\mathbb{I}^d)$ , we explicitly construct nonadaptive and adaptive deep ReLU neural networks having an output that approximates $f$ with a prescribed accuracy $\varepsilon$ , and prove dimension-dependent bounds for the computation complexity of this approximation, characterized by the size and the depth of this deep ReLU neural network, explicitly in $d$ and $\varepsilon$ . Our results show the advantage of the adaptive method of approximation by deep ReLU neural networks over nonadaptive one.

preprint2021arXiv

High-dimensional nonlinear approximation by parametric manifolds in Hölder-Nikol'skii spaces of mixed smoothness

We study high-dimensional nonlinear approximation of functions in Hölder-Nikol'skii spaces $H^α_\infty(\mathbb{I}^d)$ on the unit cube $\mathbb{I}^d:=[0,1]^d$ having mixed smoothness, by parametric manifolds. The approximation error is measured in the $L_\infty$ -norm. In this context, we explicitly constructed methods of nonlinear approximation, and give dimension-dependent estimates of the approximation error explicitly in dimension $d$ and number $N$ measuring computation complexity of the parametric manifold of approximants. For $d=2$ , we derived a novel right asymptotic order of noncontinuous manifold $N$ -widths of the unit ball of $H^α_\infty(\mathbb{I}^2)$ in the space $L_\infty(\mathbb{I}^2)$ . In constructing approximation methods, the function decomposition by the tensor product Faber series and special representations of its truncations on sparse grids play a central role.

preprint2016arXiv

B-spline quasi-interpolation sampling representation and sampling recovery in Sobolev spaces of mixed smoothness

We proved direct and inverse theorems on B-spline quasi-interpolation sampling representation with a Littlewood-Paley-type norm equivalence in Sobolev spaces $W^r_p$ of mixed smoothness $r$ , established estimates of the approximation error of recovery in $L_q$ -norm of functions from the unit ball $U^r_p$ in the spaces $W^r_p$ by linear sampling algorithms based on this representation, the asymptotic optimality of these sampling algorithms in terms of Smolyak sampling width $r^s_n(U^r_p, L_q)$ and sampling width $r_n(U^r_p, L_q)$ .

preprint2015arXiv

High-dimensional periodic sampling on Smolyak grids based on B-spline quasi-interpolation

We constructed linear algorithms of sampling recovery and cubature formulas on Smolyak grids parametrized by $m \in \mathbb{N}$ of periodic $d$ -variate functions having Lipschitz-Hölder mixed smoothness $α> 0$ based on B-spline quasi-interpolation, and studied their optimality. We established lower estimates (for $α\le 2$ ) and upper bounds of the error of the optimal sampling recovery and the optimal integration on Smolyak grids, explicit in $d$ , $m$ and the number $ν$ of active variables of functions when $d$ and $m$ may be large.

preprint2015arXiv

Hyperbolic cross approximation in infinite dimensions

We give tight upper and lower bounds of the cardinality of the index sets of certain hyperbolic crosses which reflect mixed Sobolev-Korobov-type smoothness and mixed Sobolev-analytic-type smoothness in the infinite-dimensional case where specific summability properties of the smoothness indices are fulfilled. These estimates are then applied to the linear approximation of functions from the associated spaces in terms of the $\varepsilon$ -dimension of their unit balls. Here, the approximation is based on linear information. Such function spaces appear for example for the solution of parametric and stochastic PDEs. The obtained upper and lower bounds of the approximation error as well as of the associated $\varepsilon$ -complexities are completely independent of any dimension. Moreover, the rates are independent of the parameters which define the smoothness properties of the infinite-variate parametric or stochastic part of the solution. These parameters are only contained in the order constants. This way, linear approximation theory becomes possible in the infinite-dimensional case and corresponding infinite-dimensional problems get tractable.

preprint2015arXiv

Sampling and cubature on sparse grids based on a B-spline quasi-interpolation

Let $X_n = \{x^j\}_{j=1}^n$ be a set of $n$ points in the $d$ -cube $[0,1]^d$ , and $Φ_n = \{φ_j\}_{j =1}^n$ a family of $n$ functions on $[0,1]^d$ . We consider the approximate recovery functions $f$ on $[0,1]^d$ from the sampled values $f(x^1), ..., f(x^n)$ , by the linear sampling algorithm \begin{equation} \nonumber L_n(X_n,Φ_n,f) \ := \ \sum_{j=1}^n f(x^j)φ_j. \end{equation} The error of sampling recovery is measured in the norm of the space $L_q([0,1]^d)$ -norm or the energy norm of the isotropic Sobolev sapce $W^γ_q([0,1]^d)$ for $0 < q \le \infty$ and $γ> 0$ . Functions $f$ to be recovered are from the unit ball in Besov type spaces of an anisotropic smoothness, in particular, spaces $B^a_{p,θ}$ of a nonuniform mixed smoothness $a \in {\mathbb R}^d_+$ , and spaces $B^{α,β}_{p,θ}$ of a "hybrid" of mixed smoothness $α> 0$ and isotropic smoothness $β\in \mathbb R$ . We constructed optimal linear sampling algorithms $L_n(X_n^*,Φ_n^*,\cdot)$ on special sparse grids $X_n^*$ and a family $Φ_n^*$ of linear combinations of integer or half integer translated dilations of tensor products of B-splines. We computed the asymptotic of the error of the optimal recovery. This construction is based on a B-spline quasi-interpolation representations of functions in $B^a_{p,θ}$ and $B^{α,β}_{p,θ}$ . As consequences we obtained the asymptotic of optimal cubature formulas for numerical integration of functions from the unit ball of these Besov type spaces.

preprint2014arXiv

Kolmogorov n-Widths of Function Classes Induced by a Non-Degenerate Differential Operator: A Convex Duality Approach

Let $P(D)$ be the differential operator induced by a polynomial $P$ , and let ${U^{[P]}_2}$ be the class of multivariate periodic functions $f$ such that $\|P(D)(f)\|_2\leq 1$ . The problem of computing the asymptotic order of the Kolmogorov $n$ -width $d_n({U^{[P]}_2},L_2)$ in the general case when ${U^{[P]}_2}$ is compactly embedded into $L_2$ has been open for a long time. In the present paper, we use convex analytical tools to solve it in the case when $P(D)$ is non-degenerate.

preprint2014arXiv

Lower bounds for the integration error for multivariate functions with mixed smoothness and optimal Fibonacci cubature for functions on the square

We prove lower bounds for the error of optimal cubature formulae for $d$ -variate functions from Besov spaces of mixed smoothness $B^α_{p,θ}({\mathbb G}^d)$ in the case $0 < p, θ\le \infty$ and $α> 1/p$ , where ${\mathbb G}^d$ is either the $d$ -dimensional torus ${\mathbb T}^d$ or the $d$ -dimensional unit cube ${\mathbb I}^d$ . We prove upper bounds for QMC methods of integration on the Fibonacci lattice for bivariate periodic functions from $B^α_{p,θ}({\mathbb T}^2)$ in the case $1\leq p \leq \infty$ , $0 < θ\leq \infty$ , $α>1/p$ . A non-periodic modification of this classical formula yields upper bounds for $B^α_{p,θ}({\mathbb I}^2)$ if $1/p<α<1+1/p$ . In combination these results yield the correct asymptotic error of optimal cubature formulae for functions from $B^α_{p,θ}({\mathbb G}^2)$ and indicate that a corresponding result is most likely also true in case $d>2$ . This is compared to the correct asymptotic of optimal cubature formulae on Smolyak grids which results in the observation that any cubature formula on Smolyak grids is never optimal for the general setting.

preprint2014arXiv

Sampling on energy-norm based sparse grids for the optimal recovery of Sobolev type functions in $H^γ$

We investigate the rate of convergence of linear sampling numbers of the embedding $H^{α,β} (\mathbb{T}^d) \hookrightarrow H^γ(\mathbb{T}^d)$ . Here $α$ governs the mixed smoothness and $β$ the isotropic smoothness in the space $H^{α,β}(\mathbb{T}^d)$ of hybrid smoothness, whereas $H^γ(\mathbb{T}^d)$ denotes the isotropic Sobolev space. If $γ>β$ we obtain sharp polynomial decay rates for the first embedding realized by sampling operators based on "energy-norm based sparse grids" for the classical trigonometric interpolation. This complements earlier work by Griebel, Knapek and Dũng, Ullrich, where general linear approximations have been considered. In addition, we study the embedding $H^α_{mix} (\mathbb{T}^d) \hookrightarrow H^γ_{mix}(\mathbb{T}^d)$ and achieve optimality for Smolyak's algorithm applied to the classical trigonometric interpolation. This can be applied to investigate the sampling numbers for the embedding $H^α_{mix} (\mathbb{T}^d) \hookrightarrow L_q(\mathbb{T}^d)$ for $2<q\leq \infty$ where again Smolyak's algorithm yields the optimal order. The precise decay rates for the sampling numbers in the mentioned situations always coincide with those for the approximation numbers, except probably in the limiting situation $β= γ$ (including the embedding into $L_2(\mathbb{T}^d)$ ). The best what we could prove there is a (probably) non-sharp results with a logarithmic gap between lower and upper bound.

preprint2013arXiv

Whitney's theorem for local anisotropic polynomial L_p-approximation, 0<p<1

Dinh Dũng and T. Ullrich have proven a multivariate Whitney's theorem for the local anisotropic polynomial approximation in $L_p(Q)$ for $1 \le p \le \infty$ , where $Q$ is a $d$ -parallelepiped in $\RR^d$ with sides parallel to the coordinate axes. They considered the error of best approximation of a function $f$ by algebraic polynomials of fixed degree at most $r_i - 1$ in variable $x_i,\ i=1,...,d$ . The convergence rate of the approximation error when the size of $Q$ going to 0 is characterized by a so-called total mixed modulus of smoothness. The method of proof used by these authors is not suitable to the case $0 <p<1$ . In the present paper, by a different method we proved this theorem for $0< p \le \infty$ .

preprint2011arXiv

Adaptive algorithms in sampling recovery

We study optimal algorithms in adaptive sampling recovery of smooth functions defined on the unit $d$ -cube ${\II}^d:= [0,1]^d$. The recovery error is measured in the quasi-norm $\|\cdot\|_q$ of $L_q := L_q(\II^d)$. For $B$ a subset in $L_q,$ we define a sampling algorithm of recovery with the free choice of sample points and recovering functions from $B$ as follows. For each $f$ from the quasi-normed Besov space $B^α_{p,θ}$ , we choose $n$ sample points. This choice defines $n$ sampled values. Based on these sample points and sampled values, we choose a function from $B$ for recovering $f$ . The choice of $n$ sample points and a recovering function from $B$ for each $f \in B^α_{p,θ}$ defines a $n$ -sampling algorithm $S_n^B$ by functions in $B$ . If $Φ= \{ϕ_k\}_{k \in K}$ is a family of elements in $L_q$ , let $Σ_n(Φ)$ be the non-linear set of linear combinations of $n$ free terms from $Φ,$ that is $Σ_n(Φ):= \{\, ϕ= \sum_{j=1}^n a_j ϕ_{k_j}: \ k_j \in K \, \}$ . Denote by ${\mathcal G}$ the set of all families $Φ$ in $L_q$ such that the intersection of $Φ$ with any finite dimensional subspace in $L_q$ is a finite set, and by $\Cc(B^α_{p,θ}, L_q)$ the set of all continuous mappings from $B^α_{p,θ}$ into $L_q$ . We define the quantity $$ν_n(B^α_{p,θ},L_q) := \inf_{Φ\in {\mathcal G}} \inf_{S_n^B \in \Cc(X, L_q): B= Σ_n(Φ)} \sup_{\|f\|_{B^α_{p,θ}} \le 1} \ \|f - S_n^B(f)\|_q.$$ Let $0 < p,q, θ\le \infty$ and $α> d/p$ . Then we prove the asymptotic order $ν_n(B^α_{p,θ},L_q) \asymp n^{- α/ d}.$ We also obtained the asymptotic order of quantities of optimal recovery by $S_n^B$ in terms of best $n$ -term approximation as well of other non-linear $n$ -widths.

preprint2010arXiv

B-spline quasi-interpolant representations and sampling recovery of functions with mixed smoothness

Let $ξ= \{x^j\}_{j=1}^n$ be a grid of $n$ points in the $d$ -cube ${\II}^d:=[0,1]^d$, and $Φ= \{ϕ_j\}_{j =1}^n$ a family of $n$ functions on ${\II}^d$. We define the linear sampling algorithm $L_n(Φ,ξ,\cdot)$ for an approximate recovery of a continuous function $f$ on ${\II}^d$ from the sampled values $f(x^1), ..., f(x^n)$ , by $L_n(Φ,ξ,f)\ := \ \sum_{j=1}^n f(x^j)ϕ_j$ . For the Besov class $B^α_{p,θ}$ of mixed smoothness $α$ (defined as the unit ball of the Besov space $\MB$), to study optimality of $L_n(Φ,ξ,\cdot)$ in $L_q({\II}^d)$ we use the quantity $r_n(B^α_{p,θ})_q \ := \ \inf_{H,ξ} \ \sup_{f \in B^α_{p,θ}} \, \|f - L_n(Φ,xi,f)\|_q$ , where the infimum is taken over all grids $ξ= \{x^j\}_{j=1}^n$ and all families $Φ= \{ϕ_j\}_{j=1}^n$ in $L_q({\II}^d)$. We explicitly constructed linear sampling algorithms $L_n(Φ,ξ,\cdot)$ on the grid $ξ= \ G^d(m):= \{(2^{-k_1}s_1,...,2^{-k_d}s_d) \in \II^d : \ k_1 + ... + k_d \le m\}$, with $Φ$ a family of linear combinations of mixed B-splines which are mixed tensor products of either integer or half integer translated dilations of the centered B-spline of order $r$ . The grid $G^d(m)$ is of the size $2^m m^{d-1}$ and sparse in comparing with the generating dyadic coordinate cube grid of the size $2^{dm}$ . For various $0<p,q,θ\le \infty$ and $1/p < α< r$ , we proved upper bounds for the worst case error $\sup_{f \in B^α_{p,θ}} \, \|f - L_n(Φ,ξ,f)\|_q$ which coincide with the asymptotic order of $r_n(B^α_{p,θ})_q$ in some cases. A key role in constructing these linear sampling algorithms, plays a quasi-interpolant representation of functions $f \in B^α_{p,θ}$ by mixed B-spline series.

Dinh Dũng

Research overview

Is this your work?

14 publications

Sparse-grid polynomial interpolation approximation and integration for parametric and stochastic elliptic PDEs with lognormal inputs

Weighted approximate sampling recovery and integration based on B-spline interpolation and quasi-interpolation

Computation complexity of deep ReLU neural networks in high-dimensional approximation

High-dimensional nonlinear approximation by parametric manifolds in Hölder-Nikol'skii spaces of mixed smoothness

B-spline quasi-interpolation sampling representation and sampling recovery in Sobolev spaces of mixed smoothness

High-dimensional periodic sampling on Smolyak grids based on B-spline quasi-interpolation

Hyperbolic cross approximation in infinite dimensions

Sampling and cubature on sparse grids based on a B-spline quasi-interpolation

Kolmogorov n-Widths of Function Classes Induced by a Non-Degenerate Differential Operator: A Convex Duality Approach

Lower bounds for the integration error for multivariate functions with mixed smoothness and optimal Fibonacci cubature for functions on the square

Sampling on energy-norm based sparse grids for the optimal recovery of Sobolev type functions in $H^γ$

Whitney's theorem for local anisotropic polynomial L_p-approximation, 0<p<1

Adaptive algorithms in sampling recovery

B-spline quasi-interpolant representations and sampling recovery of functions with mixed smoothness

See the researcher in context

Building this map preview

Dinh Dũng

Research overview

Is this your work?

14 publications

Sparse-grid polynomial interpolation approximation and integration for parametric and stochastic elliptic PDEs with lognormal inputs

Weighted approximate sampling recovery and integration based on B-spline interpolation and quasi-interpolation

Computation complexity of deep ReLU neural networks in high-dimensional approximation

High-dimensional nonlinear approximation by parametric manifolds in Hölder-Nikol'skii spaces of mixed smoothness

B-spline quasi-interpolation sampling representation and sampling recovery in Sobolev spaces of mixed smoothness

High-dimensional periodic sampling on Smolyak grids based on B-spline quasi-interpolation

Hyperbolic cross approximation in infinite dimensions

Sampling and cubature on sparse grids based on a B-spline quasi-interpolation

Kolmogorov n-Widths of Function Classes Induced by a Non-Degenerate Differential Operator: A Convex Duality Approach

Lower bounds for the integration error for multivariate functions with mixed smoothness and optimal Fibonacci cubature for functions on the square

Sampling on energy-norm based sparse grids for the optimal recovery of Sobolev type functions in HγH^γHγ

Whitney's theorem for local anisotropic polynomial L_p-approximation, 0<p<1

Adaptive algorithms in sampling recovery

B-spline quasi-interpolant representations and sampling recovery of functions with mixed smoothness

See the researcher in context

Building this map preview

Sampling on energy-norm based sparse grids for the optimal recovery of Sobolev type functions in $H^γ$