rd,α,γ | ψd,α,γ and ωd,α,γ | |
SPT | sγ<∞ | sγ<∞ |
PT | sγ<∞ | sγ<∞ |
QPT | γI<1 | suff.: γI<1 |
UWT | γI<1 | ? |
(σ,τ)-WT, σ∈(0,1] | γI<1 | ? |
WT | γI<1 | ? |
(σ,τ)-WT, σ>1 | no extra condition on γ | no extra condition on γ |
This paper is devoted to the study of tractability of the L2-approximation and integration from weighted Korobov spaces of increasing smoothness in the worst-case setting. The considered algorithms use information from the class Λall, including all continuous linear functionals, and from the class Λstd, including function evaluations. Necessary and sufficient conditions on the weights of the function space for strong polynomial tractability, polynomial tractability, quasi-polynomial tractability, uniform weak tractability, weak tractability and (σ,τ)-weak tractability, are provided. Our results give a comprehensive picture of the weight conditions for all standard notions of algebraic tractability. It may be helpful to study the tractability of nonhomogeneous tensor product spaces.
Citation: Weiran Ding, Jiansong Li, Yumeng Jia, Wanru Yang. Tractability of L2-approximation and integration over weighted Korobov spaces of increasing smoothness in the worst case setting[J]. Electronic Research Archive, 2025, 33(2): 1160-1184. doi: 10.3934/era.2025052
[1] | Akeel A. AL-saedi, Jalil Rashidinia . Application of the B-spline Galerkin approach for approximating the time-fractional Burger's equation. Electronic Research Archive, 2023, 31(7): 4248-4265. doi: 10.3934/era.2023216 |
[2] | María Guadalupe Morales, Zuzana Došlá, Francisco J. Mendoza . Riemann-Liouville derivative over the space of integrable distributions. Electronic Research Archive, 2020, 28(2): 567-587. doi: 10.3934/era.2020030 |
[3] | Sida Lin, Lixia Meng, Jinlong Yuan, Changzhi Wu, An Li, Chongyang Liu, Jun Xie . Sequential adaptive switching time optimization technique for maximum hands-off control problems. Electronic Research Archive, 2024, 32(4): 2229-2250. doi: 10.3934/era.2024101 |
[4] | Munirah Aali Alotaibi, Bessem Samet . A nonlinear delay integral equation related to infectious diseases. Electronic Research Archive, 2023, 31(12): 7337-7348. doi: 10.3934/era.2023371 |
[5] | Dandan Zuo, Wansheng Wang, Lulu Zhang, Jing Han, Ling Chen . Non-fragile sampled-data control for synchronizing Markov jump Lur'e systems with time-variant delay. Electronic Research Archive, 2024, 32(7): 4632-4658. doi: 10.3934/era.2024211 |
[6] | Victor Ginting . An adjoint-based a posteriori analysis of numerical approximation of Richards equation. Electronic Research Archive, 2021, 29(5): 3405-3427. doi: 10.3934/era.2021045 |
[7] | Zhengmao Chen . A priori bounds and existence of smooth solutions to a $ L_p $ Aleksandrov problem for Codazzi tensor with log-convex measure. Electronic Research Archive, 2023, 31(2): 840-859. doi: 10.3934/era.2023042 |
[8] | Ling-Xiong Han, Wen-Hui Li, Feng Qi . Approximation by multivariate Baskakov–Kantorovich operators in Orlicz spaces. Electronic Research Archive, 2020, 28(2): 721-738. doi: 10.3934/era.2020037 |
[9] | Shaoqiang Shang, Yunan Cui . Weak approximative compactness of hyperplane and Asplund property in Musielak-Orlicz-Bochner function spaces. Electronic Research Archive, 2020, 28(1): 327-346. doi: 10.3934/era.2020019 |
[10] | Ming-Kai Wang, Cheng Wang, Jun-Feng Yin . A class of fourth-order Padé schemes for fractional exotic options pricing model. Electronic Research Archive, 2022, 30(3): 874-897. doi: 10.3934/era.2022046 |
This paper is devoted to the study of tractability of the L2-approximation and integration from weighted Korobov spaces of increasing smoothness in the worst-case setting. The considered algorithms use information from the class Λall, including all continuous linear functionals, and from the class Λstd, including function evaluations. Necessary and sufficient conditions on the weights of the function space for strong polynomial tractability, polynomial tractability, quasi-polynomial tractability, uniform weak tractability, weak tractability and (σ,τ)-weak tractability, are provided. Our results give a comprehensive picture of the weight conditions for all standard notions of algebraic tractability. It may be helpful to study the tractability of nonhomogeneous tensor product spaces.
Solving multivariate continuous problems is a classical topic in applications, and there are thousands of papers to study this problem, usually for a fixed space of functions. Such problems can almost never be solved analytically. Since they must be solved numerically, they can only be solved approximately to within a threshold ε. To deal with these problems, we often use algorithms based on finitely many information evaluations, either from the class Λall of general linear information consisting of all continuous linear functionals, or from the class Λstd of standard information consisting of function evaluations only. This motivates us to study the tractability of multivariate problems referring to S={Sd:Fd→Gd}, where Fd is a Banach space of functions and Gd is another Banach space. For ε∈(0,1), the information complexity nX(ε,Sd,Fd;Λ), X∈{ABS,NOR}, can be defined as the least number of linear functionals that are sufficient to obtain an ε-approximation for the information class Λ under the absolute error criterion (ABS) and normalized error criterion (NOR). Tractability describes how the information complexity nX(ε,Sd,Fd;Λ) behaves as a function of d and ε−1, which mainly includes strong polynomial tractability (SPT), polynomial tractability (PT), quasi-polynomial tractability (QPT), uniformly weak tractability (UWT), weak tractability (WT), and (σ,τ)-weak tractability ((σ,τ)-WT).
The two most important and widely studied such problems are multivariate approximation and multivariate integration (see e.g., [1,2,3]). In many applications, including the solution of important computational problems such as differential and integral equations, and problems in financial mathematics, people are particularly interested in Sobolev spaces. A variant of Sobolev spaces is the Korobov space, which is often called the Sobolev spaces of dominating mixed smoothness. Such spaces are probably the most important spaces for the study of computational problems for periodic smooth functions. Moreover, these spaces also have many interesting applications for non-periodic functions due to interesting relations and estimates between the complexity for the periodic and non-periodic computational problems, such as multivariate integration, see e.g., [2]. This paper is devoted to discussing the tractability of multivariate L2-approximation
APP:={APPd:H(KR)→L2([0,1]d)}d∈N |
with APPd(f)=f, and multivariate integration
INT:={INTd:H(KR)→R}d∈N |
with INTd(f)=∫[0,1]df(x)dx, for both Λall and Λstd in the worst-case setting, where H(KR) denotes a weighted Korobov space with increasing smoothness defined over [0,1]d with Fourier weight R∈{rd,α,γ,ψd,α,γ,ωd,α,γ} (see Section 2.2 for details). Here, γ={(1≥)γ1≥γ2≥…} is a positive weight sequence and α={(1<)α1≤α2≤…} is an increasing smoothness sequence (see Section 2 for details). We remark that in our considered case the initial error is 1; therefore, the results under ABS and NOR coincide. The related problem has already been discussed in a large number papers, see e.g., [1,3,4,5,7]. Particularly, for R=rd,α,γ with 1<α1=α2=…, the complete picture of multivariate L2-approximation for Λ∈{Λall,Λstd} and multivariate integration was given in [8]. For R=rd,α,γ with 1<α1≤α2≤…, the results on multivariate L2-approximation for Λall were given in [9]. For a more general case, R∈{ψd,α,γ,ωd,α,γ}, some partial results of multivariate L2-approximation for Λall were found in [10]. The related results will be formulated in Section 2.3. In this paper, our aim is to provide the conditions that are both sufficient and necessary for all tractability notions of multivariate L2-approximation and multivariate integration for both information classes Λall and Λstd.
The remainder of this article is organized as follows. In Section 2 we recall some basic facts on tractability and weighted Korobov spaces, and briefly formulate previous results and the main results of the current paper. In Sections 3 and 4, we discuss the tractability of L2-approximation and integration in the worst-case setting, respectively.
In this section, we first introduce the fundamental concepts related to tractability in multivariate problems and weighted Korobov spaces, and then we recall the previous results on the tractability of L2-approximation and integration. Finally, our results are summarized in a table.
Let {Fd}d∈N, {Gd}d∈N be two sequences of Hilbert spaces, and let {Sd:Fd→Gd}d∈N be a family of continuous linear operators. We will consider two special choices of S, namely:
● L2-approximation of functions f∈Fd. In this case, we have Sd(f)=f and Gd=L2([0,1]d).
● Integration of functions f∈Fd. In this case, we have Sd(f)=∫[0,1]df(x)dx and Gd=R.
We approximate Sd by algorithms An,d of the form
An,d(f)=n∑i=1Li(f)gi, | (2.1) |
where gi∈Gd and Li∈Λ for i=1,⋯,n. We will assume that the considered functionals Li belong to Λall=F∗d or Λstd. The worst-case error for the algorithm An,d of the form (2.1) is defined as
e(An,d,Sd,Fd):=sup‖f‖Fd≤1‖Sd(f)−An,d(f)‖Gd. |
Then, the n-th minimal worst-case error is defined by
e(n,Sd,Fd;Λ):=infAn,dL1,…,Ln∈Λe(An,d,Sd,Fd), |
where the infimum is taken over all linear algorithms of the form (2.1). Particularly, for n=0, the initial error of the problem Sd in the worst-case setting is defined by
e(0,Sd,Fd;Λ)≡e(0,Sd,Fd)=sup‖f‖Fd≤1‖Sd(f)‖Gd=‖Sd‖Fd→Gd. |
It is interesting to see how the worst-case errors of An,d depend on the number n and on the dimension d under ABS or NOR. For this purpose, we introduce the so-called information complexity as
nX(ε,Sd,Fd;Λ):=min{n∈N0:e(n,Sd,Fd;Λ)≤εCRIXd} |
where
CRIXd:={1,if X= ABS,e(0,Sd,Fd),if X= NOR. |
The following notions have been frequently studied. For more about tractability we refer to [1,2,3] by Novak and Woźniakowski.
Consider the multivariate problem S:={Sd}d∈N using the information class Λ∈{Λall,Λstd} under ABS or NOR. We say that S is
1) polynomially tractable (PT) if there are C>0 and τ,σ≥0 satisfying
n(ε,Sd,Fd;Λ)≤Cε−τdσ for all d∈N,ε∈(0,1). |
2) strong polynomially tractable (SPT) if there are C>0 and τ≥0 satisfying
n(ε,Sd,Fd;Λ)≤Cε−τ for all d∈N,ε∈(0,1). | (2.2) |
The infimum of τ≥0 for which the inequality (2.2) is satisfied for a certain C>0 is referred to as the exponent of SPT, denoted by τ∗(Λ).
3) quasi-polynomially tractable (QPT) if there are t≥0 and C>0 satisfying
n(ε,Sd,Fd;Λ)≤Cexp(t(1+lnd)(1+lnε−1)) for all d∈N,ε∈(0,1). | (2.3) |
The infimum of t≥0 for which the inequality (2.3) is satisfied for a certain C>0 is referred to as the exponent of QPT, denoted by τ∗(Λ).
4)weak tractable (WT) if
limd+ε−1→∞lnn(ε,Sd,Fd;Λ)d+ε−1=0. |
5) (σ,τ)-weak tractable ((σ,τ) -WT) for σ,τ>0 if
limd+ε−1→∞lnn(ε,Sd,Fd;Λ)dσ+ε−τ=0. |
6)uniform weak tractable (UWT) if (σ,τ)-WT holds for all σ,τ>0.
It is evident that the following relationship holds:
SPT⇒PT⇒QPT⇒UWT⇒(σ,τ)-WTfor all (σ,τ)∈(0,∞)2. |
Clearly, the notions of WT and (1,1)-WT coincide.
Now we briefly recall some facts on weighted Korobov spaces of increasing smoothness. Weighted Korobov spaces are special types of the reproducing kernel Hilbert spaces H(KR) with a reproducing kernel of the form
KR(x,y):=∑k∈ZdR(k)exp(2πik⋅(x−y)), for all x,y∈[0,1]d, |
for some summable function R:Zd→(0,+∞), i.e., ∑k∈ZdR(k)<+∞, which is often called the Fourier weight, and the corresponding inner product
⟨f,g⟩H(KR):=∑k∈Zd1R(k)ˆf(k)¯ˆg(k)and‖f‖H(KR)=√⟨f,f⟩H(KR) . |
Here, the Fourier coefficients {ˆf(k)}k∈Zd are given by
ˆf(k)=∫[0,1]df(x)exp(−2πik⋅x)dx. |
We remark that the reproducing kernel KR is well defined, since for all x,y∈[0,1]d,
|KR(x,y)|≤∑k∈ZdR(k)<+∞. |
We will consider three possible Fourier weights. Assume that γ={γj}j∈N is a sequence of so-called product weights and α={αj}j∈N is the smoothness parameter sequence satisfying
1<α1≤α2≤⋯. | (2.4) |
1) Weighted Korobov space.
In this case, we consider R(k) as
rd,α,γ(k)=d∏j=1rαj,γj(kj), |
where
rα,γ(k):={1,for k=0,γ1|k|α,for k≠0, |
for γ>0,α>1.
Furthermore, when αj∈N for all j∈N, we define the norm
‖f‖2rd,α,γ:=∑τi∈{0,…,αi}i=1,…,d(d∏j=1τj≠0γ−1j)∫[0,1]d(∂τxf(x))2dx, |
where ∂τx:=∂τ1x1…∂τdxd for τ:=(τ1,…,τd)∈Zd. It can be checked that the norm ‖f‖H(Krd,α,γ) is finite if and only if the norm ‖f‖rd,α,γ is finite with the existence of ∂τxf for all τ≤α:=(α1,…,αd). This also indicates a relation between smoothness and α.
2) A first variant of the weighted Korobov space.
In this case, we consider R(k) as
ψd,α,γ(k)=d∏j=1ψαj,γj(kj), |
where
ψα,γ(k):={1,for k=0,γ1|k|!,for 1≤|k|<α,γ(|k|−⌈α⌉)!|k|!,for |k|≥α, |
for γ>0,α>1.
3) A second variant of the weighted Korobov space.
In this case, we consider R(k) as
ωd,α,γ(k)=d∏j=1ωαj,γj(kj), |
where
ωα,γ(k):=(1+1γ⌈α⌉∑s=1βs(k))−1, βs(k):={|k|!(|k|−s)!,if |k|≥s,0,otherwise, |
for γ>0,α>1.
Without loss of generality, we usually assume that the weights γ={γj}j∈N are ordered by
1≥γ1≥γ2≥⋯>0. | (2.5) |
Remark 1 ([[10]]). Let Rαj,γj∈{rαj,γj,ψαj,γj,ωαj,γj} for all j∈N. Then
Rαj,γj(0)=1 and Rαj,γj(1)≥γj2. | (2.6) |
Furthermore, we have the following proposition from [11] (also [10, Remark 4]).
Proposition 1 ([10,11]). Let α={αj}j∈N be the smoothness parameter sequence satisfying (2.4) and let γ={γj}j∈N be the weight sequence satisfying (2.5). Set Rd,α,γ={Rαj,γj}j≥1∈{rd,α,γ,ψd,α,γ,ωd,α,γ}. Then, for all j∈N and k∈Z, we have
rαj,γj(k)≤ψαj,γj(k)≤⌈αj⌉⌈αj⌉rαj,γj(k), |
and
13rαj,γj(k)≤ωαj,γj(k)≤ψαj,γj(k)≤⌈αj⌉⌈αj⌉rαj,γj(k). | (2.7) |
Hence, we obtain
13rαj,γj(k)≤Rαj,γj(k)≤⌈α1⌉⌈α1⌉rα1,γj(k), | (2.8) |
which implies that for all k∈Zd,
rd,α,γ/3(k)≤Rd,α,γ(k)≤Rd,α1,γ(k)≤rd,α1,⌈α1⌉⌈α1⌉γ(k), | (2.9) |
and
rd,α,γ(k),ωd,α,γ(k)≤ψd,α,γ(k), | (2.10) |
where rd,α,γ/3=∏dj=1rαj,γj/3 and rd,α1,⌈α1⌉⌈α1⌉γ=∏dj=1rα1,⌈α1⌉⌈α1⌉γj.
Remark 2. From the definition of H(KR), we know that if R1,R2:Zd→(0,+∞) are two Fourier weights for weighted Korobov spaces H(KR1) and H(KR2) with
R1(k)≤CR2(k), for all k∈Zd, |
then
C‖f‖H(KR1)≥‖f‖H(KR2), for all f∈H(KR2). |
That is, H(KR1) is continuously embedded into H(KR2), denoted by H(KR1)↪H(KR2). Meanwhile, by (2.7), we have
rd,α,γ(k)≤ψd,α,γ(k)≤(d∏j=1⌈αj⌉⌈αj⌉)rd,α,γ(k), |
and
rd,α,γ/3(k)≤ωd,α,γ(k)≤(d∏j=1⌈αj⌉⌈αj⌉)rd,α,γ(k), |
which implies that
H(Krd,α,γ)↪H(Kψd,α,γ)↪H(Krd,α,γ), |
and
H(Krd,α,γ/3)↪H(Kωd,α,γ)↪H(Krd,α,γ). |
These also indicate a relation between smoothness and α.
Let γ={γj}j∈N be a sequence of product weights and let α={αj}j∈N be the smoothness parameter sequence satisfying 1<α1≤α2≤⋯. For the weight space H(Krd,α,γ) with 1<α1=α2=…, the complete picture of tractability of L2-approximation APP for both Λall and Λstd, and integration INT, has been obtained in [12]. Besides, [5,10] obtained some partial results of weighted Korobov space H(KR) with R∈{rd,α,γ,ψd,α,γ,ωd,α,γ}. The characterization of tractability will be given in terms of decay conditions on the weight sequence γ={γj}j∈N. Before presenting our findings, we will first outline some key notations. (We use the convention that inf∅=∞.)
● The infimum of the sequence γ is denoted by
γI:=infj≥1γj. |
● The infimum of the sequence (lnγ−1j/lnj)j∈N is denoted by
δγ:=lim infj→∞lnγ−1jlnj. |
● The sum exponent sγ is defined as
sγ:=inf{κ>0:∞∑j=1γκj<∞}. |
● The exponent tγ is defined as
tγ:=inf{κ>0:lim supd→∞1lndd∑j=1γκj<∞}. |
We summarize the above results in the following two tables (Tables 1 and 2).
Remark 3. We remark that the condition δγ>0 is equivalent to the condition sγ<∞ which can be seen from the condition 1<α=α1=α2=⋯ in Tables 1 and 2.
rd,α,γ | ψd,α,γ and ωd,α,γ | |
SPT | sγ<∞ | sγ<∞ |
PT | sγ<∞ | sγ<∞ |
QPT | γI<1 | suff.: γI<1 |
UWT | γI<1 | ? |
(σ,τ)-WT, σ∈(0,1] | γI<1 | ? |
WT | γI<1 | ? |
(σ,τ)-WT, σ>1 | no extra condition on γ | no extra condition on γ |
rd,α,γ | ψd,α,γ and ωd,α,γ | |
SPT | δγ>0 | δγ>0 |
PT | δγ>0 | δγ>0 |
QPT | γI<1 | ? |
UWT | γI<1 | ? |
(σ,τ)-WT, σ∈(0,1] | γI<1 | ? |
WT | γI<1 | ? |
(σ,τ)-WT, σ>1 | no extra condition on γ | no extra condition on γ |
In this paper, we outline the conditions necessary for the tractability of L2-approximation in H(KR) for both Λall and Λstd, as demonstrated in Theorems 1 and 2. The results of our findings are presented in Table 3.
Λall | Λstd | |
SPT | sγ<∞ | ∑∞j=1γj<∞ |
PT | sγ<∞ | lim supd→∞1lnd∑dj=1γj<∞ |
QPT | γI<1 | lim supd→∞1lnd∑dj=1γj<∞ |
UWT | γI<1 | limd→∞1dσ∑dj=1γj=0 ∀σ∈(0,1] |
(σ,τ)-WT for σ∈(0,1] | γI<1 | limd→∞1dσ∑dj=1γj=0 |
WT | γI<1 | limd→∞1dd∑j=1γj=0 |
(σ,τ)-WT for σ>1 | no extra condition on γ | no extra condition on γ |
This paper also investigates the tractability of integration INT. We obtain the sufficient and necessary conditions on all tractability notions, which are the same as for the tractability of the L2-approximation APP for the class Λall.
This section focuses on the tractability of L2-approximation on the weighted Korobov space H(KR). We will study this problem for the classes Λall and Λstd, respectively.} More precisely, we want to approximate the embedding operators
APPd:H(KR)→L2([0,1]d), APPd(f)=f, |
in the worst-case setting. It can be seen that APPd is the compact embedding from the weighted space H(KR) to L2([0,1]d). To approximate APPd with respect to the L2-norm ‖⋅‖L2 over [0,1]d, it follows from [1,13] that it suffices to use linear algorithms Aappn,d of the form
Aappn,d(f)=n∑i=1Li(f)gi | (3.1) |
with gi∈L2([0,1]d) and Li∈{Λall,Λstd} for i=1,…,n. Remember that H(KR) is a reproducing kernel Hilbert space and Λstd⊆Λall.
The n-th minimal worst-case error with respect to the information class Λ is given by
e(n,APPd,H(KR);Λ):=infAappn,dL1,…,Ln∈Λsup‖f‖H(KR)≤1‖f−Aappn,d(f)‖L2, |
where the infimum is extended over all algorithms of the form (3.1) with information from Λ. For n=0, the initial approximation error is given by
e(0,APPd,H(KR)):=sup‖f‖H(KR)≤1|APPd(f)|=‖APPd‖. |
In the following, we assume that 0<R(k)≤1 for all k∈Zd, which implies that the norm of APPd is 1, since for all f∈H(KR),
‖APPd(f)‖2L2=‖f‖2L2=∑k∈Zd|ˆf(k)|2≤∑k∈Zd1R(k)|ˆf(k)|2=‖f‖2H(KR)<+∞, |
and the above inequality is sharp if f≡1. In the remaining part of this paper, we always assume that 0<R(k)≤1, which implies that there is no need to distinguish between ABS and NOR. For abbreviation, we write n(ε,APPd,H(KR);Λ)≡nX(ε,APPd,H(KR);Λ).
First we present the results on L2-approximation in the weighted Korobov space H(KR) for the information class Λall.
Theorem 1. Let α={αj}j∈N be the smoothness parameter sequence satisfying (2.4) and let γ={γj}j∈N be the weight sequence satisfying (2.5). Consider the L2-approximation problem APP={APPd}d∈N for the weighted spaces H(KR), R∈{rd,α,γ,ψd,α,γ,ωd,α,γ} for the information class Λall under ABS and NOR. We have the following conditions:
(i) (Cf. [10]) SPT and PT are equivalent and hold iff sγ<∞. In this case, the exponent of SPT is
τ∗(Λall)=2max(1α1,sγ). |
(ii) QPT, UWT, WT, and (σ,τ)-WT with σ∈(0,1], are are equivalent and hold iff γI<1. In this case, the exponent of QPT satisfies
t∗(Λall)≤2max(1α1,1lnγ−1I), |
and the equality holds for R∈{rd,α,γ,ψd,α,γ}. (In particular, if γI=0, then we set (lnγ−1I)−1:=0 and we have that t∗(Λall)=2/α1.)
(iii) (Cf. [10]) For σ>1, (σ,τ)-WT always holds.
To obtain the above theorem, we recall some basic knowledge about L2-approximation for the information class Λall. It is well known, see e.g., [1], that the n-th minimal worst-case errors e(n,APPd,H(KR);Λall) and the information complexity n(ε,APPd,H(KR);Λall) depend on the eigenvalues of the continuous linear operator
Wd=APP∗dAPPd:H(KR)→H(KR). |
Let (λd,j,ηd,j)j∈N be the eigenpairs of Wd, i.e.,
Wdηd,j=λd,jηd,j, for all j∈N, |
where the eigenvalues (λd,j) are ordered by λd,1≥λd,2≥⋯≥0, and the eigenvectors (ηd,j)j∈N are orthonormal,
⟨ηd,i,ηd,j⟩H(KR)=δi,j, for all i,j∈N. |
Then, the n-th minimal error is obtained (see [1, Corollary 4.12]) for the algorithm
A∗n,df=n∑j=1⟨f,ηd,j⟩H(KR)ηd,j, for all f∈H(KR), |
and
e(n,APPd,H(KR);Λall)=e(A∗n,d,H(KR),Λall)=λ1/2d,n+1. | (3.2) |
Hence, the information complexity is equal to
n(ε,APPd,H(KR);Λall)=|{n∈N0:λd,n>ε2}|, |
with ε∈(0,1) and d∈N. Since it follows from [1] that the eigenpairs of the operator Wd are (R(k),ek) with k∈Zd, where
ek(x):=√R(k)exp(2πik⋅x), |
we have
n(ε,APPd,H(KR);Λall)=|{k∈Zd:R(k)>ε2}|. | (3.3) |
Now we begin to prove Theorem 1.
Proof of Theorem 1. (ⅰ) and (ⅲ) have been proved in [10]. We include the proof of item (ⅰ) and (ⅱ) in a different way as a warm-up.
(ⅰ) From [1, Theorem 5.2], we know that APP is PT for Λall if and only if there exist q≥0 and τ>0 such that
Cτ,q:=supd∈N(∞∑j=1λτd,j)1/τd−q<∞, | (3.4) |
and τ∗(Λall)=inf{2τ:τ satisfies (3.4)}.
Since SPT implies PT, it suffices to demonstrate the sufficiency of sγ<∞ for SPT, and the necessary of sγ<∞ for PT.
● Sufficiency of sγ<∞ for SPT.
Now we assume that sγ<∞. Then, for any τ>max(sγ,1α1), we have ∑∞j=1γτj<∞. Using (2.10) and (2.8), we obtain
∞∑j=1λτd,j=d∏j=1(1+2∞∑k=1(Rαj,γj(k))τ)≤d∏j=1(1+2⌈α1⌉⌈α1⌉τ∞∑k=1(rαj,γj(k))τ)≤d∏j=1(1+2⌈α1⌉⌈α1⌉τγτjζ(αjτ))=exp{d∑j=1ln(1+2⌈α1⌉⌈α1⌉τγτjζ(αjτ))}≤exp{2⌈α1⌉⌈α1⌉τd∑j=1ζ(αjτ)γτj}≤exp{2⌈α1⌉⌈α1⌉τζ(α1τ)∞∑j=1γτj}<∞, |
where we also used that ln(1+x)≤x for x≥0 and ζ(α1τ)<∞ for τ>1/α1. This leads to Cτ,q<∞ for all q≥0, thus we have SPT, PT and
τ∗(Λall)≤2max(sγ,1α1). | (3.5) |
● Necessity of sγ<∞ for PT.
Now we assume that PT holds. Then there exists τ>0 such that (3.4) holds with q=0. Clearly, we need τ>1/α1. Again using (2.8) and ζ(αjτ)>1 for all j, we obtain
∞∑j=1λτd,j=∑k∈Zd(Rd,α,γ(k))τ=d∏j=1(1+2∞∑k=1(Rαj,γj(k))τ)≥d∏j=1(1+23τ∞∑k=1(rαj,γj(k))τ)=d∏j=1(1+23τγτjζ(αjτ))≥13τd∑j=1γτj. |
Meanwhile, (3.4) implies that ∑∞j=1γτj<∞ and hence sγ≤τ<∞. Combining both results yields that τ≥max(sγ,1α1) and hence also
τ∗(Λall)≥2max(sγ,1α1). | (3.6) |
Then (3.5) and (3.6) imply τ∗(Λall)=2max(sγ,1α1).
(ⅱ) From [3, Theorem 23.2] (see also [14]), we know that APP is QPT if and only if there exists a τ>0 such that
Cτ:=supd∈NCτ,d:=supd∈N1d2(∞∑j=1λτ(1+lnd)d,j)1/τ<∞. | (3.7) |
First, assume that γI<1. For R∈{rd,α,γ,ψd,α,γ,ωd,α,γ}, we have, by (2.7) and similar to the proof of [11, Theorem 14], that
∞∑j=1λτ(1+lnd)d,j=∑k∈Zd(Rd,α,γ(k))τ(1+lnd)≤∑k∈Zd(ψd,α,γ(k))τ(1+lnd)=d∏j=1(1+2∞∑k=1(ψαj,γj(k))τ(1+lnd))=d∏j=1(1+2γτ(1+lnd)j(⌈αj⌉−1∑k=1(1k!)τ(1+lnd)+∞∑k=⌈αj⌉((k−⌈αj⌉)!k!)τ(1+lnd)))≤d∏j=1(1+2γτ(1+lnd)j(⌈αj⌉−1∑k=1(1k!)τ(1+lnd)+∞∑k=⌈αj⌉(1(k−⌈αj⌉+1)αj)τ(1+lnd)))=d∏j=1(1+2γτ(1+lnd)j(⌈αj⌉−1∑k=1(1k!)τ(1+lnd)+∞∑k=1(1kαj)τ(1+lnd)))=d∏j=1(1+2γτ(1+lnd)j(⌈αj⌉−1∑k=1(1k!)τ(1+lnd)+ζ(αjτ(1+lnd))))≤d∏j=1(1+2γτ(1+lnd)j(∞∑k=1(1k!)τ(1+lnd)+ζ(α1τ(1+lnd)))). |
In order that ζd:=ζ(α1τ(1+lnd))<∞ for all d∈N, we need τ>1/α1. Next, we have
∞∑k=1(1k!)τ(1+lnd)=1+12τ(1+lnd)+16τ(1+lnd)+∞∑k=4(1k!)τ(1+lnd)≤1+12τ(1+lnd)+16τ(1+lnd)+∞∑k=4(12τ(1+lnd))k=1+12τ(1+lnd)+16τ(1+lnd)+124τ(1+lnd)11−1/2τ(1+lnd)=1+12τ(1+lnd)+16τ(1+lnd)+124τ(1+lnd)−23τ(1+lnd)≤1+12τ(1+lnd)+16τ(1+lnd)+123τlnd 124τ−23τ≤1+12τlnd(2+124τ−23τ)=1+cτdτln2, |
where cτ:=2+1/(24τ−23τ). This gives
Cτ,d≤1d2{d∏j=1(1+2γτ(1+lnd)j(cτdτln2+ζd))}1/τ=exp{1τd∑j=1ln(1+2γτ(1+lnd)j(ζd+cτdτln2))−2lnd}≤exp{2τ(ζd+cτdτln2)d∑j=1γτ(1+lnd)j−2lnd}, |
where we used that ln(1+x)≤x for all x≥0. Using the inequality
ζ(x)≤1+1x−1, for all x>1, |
we have
ζd≤1+1(α1τ−1)+α1τlnd. |
Then we obtain
Cτ,d≤exp{2τ(1+1(α1τ−1)+α1τlnd+cτdτln2)d∑j=1γτ(1+lnd)j−2lnd}. |
Now we distinguish two cases:
● Case 1: γI=0. In this case, we have limj→∞γj=0 and hence, there is a J∈N such that γj≤e−1/τ for all j≥J. Then
d∑j=1γτ(1+lnd)j≤J−1∑j=11+d∑j=Je−lnd=J. |
Note that J depends on τ, and it is finite for any fixed τ. Thus, if τ>1/α1 and γI=0 then
Cτ,d≤exp{2τ(1+1(α1τ−1)+α1τlnd+cτdτln2)J−2lnd}→0 if d→∞. |
By the characterization in (3.7), this implies QPT.
● Case 2: γI∈(0,1). In this case, for any γ∗∈(γI,1), there exists a j0∈N such that γj≤γ∗<1 for all j>j0. Then for every d∈N,
d∑j=1γτ(1+lnd)j≤j0+γτ(1+lnd)∗max(d−j0,0)=j0+γτ∗max(d−j0,0)dτlnγ−1∗≤j0+1, |
whenever τ≥(lnγ−1∗)−1. Thus, if τ>1/α1 and τ≥(lnγ−1∗)−1, then
Cτ,d≤exp{2τ(1+1(α1τ−1)+α1τlnd+cτdτln2)(j0+1)−2lnd}→0,if d→∞. |
Again, by the characterization in (3.7), this implies QPT.
Clearly, QPT implies UWT, WT, and (σ,τ)-WT with σ∈(0,1]. So it remains to prove that WT implies γI<1. To this end, we let γI=1, i.e., γj≡1. Then we derive from (2.6) that for k∈{0,1}d,
Rd,α,γ(k)=d∏j=1Rαj,γj(kj)≥d∏j=1(γj2)=:η, |
which, by (3.3), implies that for any ε∈(0,√η),
n(ε,APPd,H(KR);Λall)≥|{k∈{0,1}d:Rd,α,γ(k)>ε2}|=2d. |
This means that APP suffers from the curse of dimensionality, and thus, we cannot have (σ,τ)-WT for any σ∈(0,1]. This deduces that QPT, UWT, WT, and (σ,τ)-WT with σ∈(0,1] are equivalent, and hold iff γI<1.
Now we turn to calculate the exponent of QPT. Again from [3], Theorem 23.2], we know that the exponent of QPT is
t∗(Λall)=inf{2τ:τ satisfies (3.7)}. |
From the above part of the proof, it follows that τ satisfies (3.7) as long as τ>1/α1 and τ>(lnγ−1I)−1, where we put (lnγ−1I)−1:=0 whenever γI=0. Therefore,
t∗(Λall)≤2max(1α1,1lnγ−1I). |
Assume now that we have QPT for R∈{rd,α,γ,ψd,α,γ}. Then, (3.7) holds true for some τ>0. Consideration of the special case d=1 together with (3.7) yields that
Cτ≥(∞∑j=1λτ1,j)1/τ=(1+2∞∑k=1Rτα1,γ1(k))1/τ≥(1+2∞∑k=1rτα1,γ1(k))1/τ=(1+2ζ(α1τ)γτ1)1/τ, |
where we used (2.10), and hence we must have τ>1/α1. This already implies the result t∗(Λall)=2/α1 whenever γI=0.
It remains to study the case γI>0. Now, again according to (3.7) and (2.10), there exists a τ>1/α1 such that for all d∈N we have
Cτ≥1d2{d∏j=1(1+∞∑k=1(Rαj,γj(k))τ(1+lnd))}1/τ≥1d2{d∏j=1(1+∞∑k=1(rαj,γj(k))τ(1+lnd))}1/τ=1d2{d∏j=1(1+γτ(1+lnd)jζ(αjτ(1+lnd)))}1/τ≥1d2{d∏j=1(1+γτ(1+lnd)j)}1/τ=exp{1τd∑j=1ln(1+γτ(1+lnd)j)−2lnd}. |
Taking the logarithm leads to
lnCτ≥1τd∑j=1ln(1+γτ(1+lnd)j)−2lnd≥dτln(1+γτ(1+lnd)I)−2lnd |
for all d∈N. Since γI∈(0,1) and since ln(1+x)≥xln2 for all x∈[0,1], it follows that for all d∈N we have
lnCτ≥dln2τγτ(1+lnd)I−2lnd=γτIln2τd1−τlnγ−1I−2lnd, |
which implies that τ≥(lnγ−1I)−1. Therefore, we will get
t∗(Λall)≥2max(1α1,1lnγ−1I) |
and the claimed result follows.
(ⅲ) The result for (σ,τ)-weak tractability for σ>1 for the class Λall follows from the corresponding result for the class Λstd from Theorem 2.
The proof is completed.
In the next theorem, we present the respective conditions for tractability of L2-approximation in the weighted Korobov space H(KR) for the information class Λstd.
Theorem 2. Let α={αj}j∈N be the smoothness parameter sequence satisfying (2.4) and let γ={γj}j∈N be the weight sequence satisfying (2.5). Consider the L2-approximation problem APP={APPd}d∈N for the weighted spaces H(KR), R∈{rd,α,γ,ψd,α,γ,ωd,α,γ} for the information class Λstd under ABS and NOR. Then we have the following conditions:
(i) SPT holds iff
∞∑j=1γj<∞ |
(which implies sγ≤1). In this case, the exponent of SPT is
τ∗(Λstd)=2max(1α1,sγ). |
(ii) PT, QPT are equivalent and hold iff
lim supd→∞1lndd∑j=1γj<∞. | (3.8) |
(iii) WT holds iff
limd→∞1dd∑j=1γj=0. | (3.9) |
(iv) UWT holds iff
limd→∞1dσd∑j=1γj=0, for all σ∈(0,1]. | (3.10) |
(v) For σ∈(0,1], (σ,τ)-WT holds iff
limd→∞1dσd∑j=1γj=0. | (3.11) |
For σ>1, (σ,τ)-WT always holds.
The proof of Theorem 2 is based on relations between the minimal errors of Λstd and Λall, in particular on [17, Theorem 1] and on [6, Theorem 1] (see also [3, Theorem 26.10]). These results provide that the trace of the operator Wd=APP∗dAPPd is finite. Recall that the trace of Wd is given by the sum of its eigenvalues, i.e.,
trace(Wd)=∞∑j=1λd,j=∑k∈ZdRd,α,γ(k)=d∏j=1(1+2∞∑k=1Rαj,γj(k)). |
Using (2.10), we obtain
trace(Wd)≤d∏j=1(1+2⌈α1⌉⌈α1⌉∞∑k=1rαj,γj(k))=d∏j=1(1+2⌈α1⌉⌈α1⌉ζ(αj)γj)≤d∏j=1(1+2⌈α1⌉⌈α1⌉ζ(α1)γj) | (3.12) |
and
trace(Wd)≥d∏j=1(1+23∞∑k=1rαj,γj(k))=d∏j=1(1+2γj3ζ(αj)). |
Hence trace(Wd) is infinite if and only if α1=α2=⋯=1. Note that in general there is no relation between the power of Λall and Λstd whenever the trace of Wd is infinite. For a discussion of this issue we refer to [3, Section 26.3].
In order to prove the sufficiency of PT for Λstd, we use a criterion from [3] as stated below.
Lemma 2. ([3, Theorem 26.13]). Consider multivariate approximation APP={APPd}d∈N in the worst-case setting for ABS or NOR, where
APPd:Fd→Gd with APPd(f)=f, |
and Fd is a reproducing kernel Hilbert space continuously embedded in Gd. Assume that the trace of Wd is finite for all d, and there are two constants C,q≥0 such that
trace(Wd)(CRIXd)2≤Cdq, for all d∈N. |
Then polynomial tractabilities of APP for Λall and Λstd are equivalent. Particularly, if q=0, then strong polynomial tractabilities of APP for Λall and Λstd are equivalent. More precisely,
nX(ε,APPd,Fd;Λall)≤Callε−palldqall for all ε∈(0,1),d∈N, |
implies that there exists a constant Cstd≥0 such that
nX(ε,APPd,Fd;Λstd)≤Cstdε−pstddqstd for all ε∈(0,1),d∈N, |
where
pstd=pall+2 and qstd=qall+q. |
Now we prove Theorem 2.
Proof of Theorem 2. It follows from (4.2) in Section 4 that necessary conditions of tractability for integration INT are also necessary for L2 approximation APP. Consequently, the necessity of the conditions outlined in items (ⅰ)–(ⅴ) is derived from Theorem 3. It is enough to examine whether the conditions in items (ⅰ)–(ⅴ) are sufficient.
(ⅰ) Clearly, α1>1 implies that trace(Wd) is finite for all d∈N. Then according to [17, Theorem 1], there is a positive integer c satisfying for all n∈N
e(cn,APPd,H(KR);Λstd)2≤1n∞∑k=ne(k,APPd,H(KR);Λall)2. | (3.13) |
Now let ∑∞j=1γj<∞. Then, the sum exponent sγ≤1.
● First we consider sγ<1. Then, using Theorem 1, we have SPT for Λall with exponent
τ∗(Λall)=2max(1α1,sγ)<2. |
Hence, for every τ>τ∗(Λall) there is C>0 such that
n(ε,APPd,H(KR);Λall)≤Cε−τ, |
which yields that
e(k,APPd,H(KR);Λall)≤Ck1/τ. |
Inserting into (3.13) gives
e(cn,APPd,H(KR);Λstd)2≤Cn∞∑k=n1k2/τ≤Cn∫∞n1x2/τdx=Cnτ2−τ1n2/τ−1≤Cτ2−τ1n2/τ. |
Hence, there exists a number aτ>0 such that
e(cn,APPd,H(KR);Λstd)≤aτn1/τ, |
which implies that
n(ε,APPd,H(KR);Λstd)≤⌈caττε−τ⌉. |
Hence, since τ>τ∗(Λstd) is arbitrary, we have SPT with
τ∗(Λstd)=2max(1α1,sγ). |
(Note that trivially τ∗(Λstd)≥τ∗(Λall)=2max(1/α1,sγ).)
● Next we consider sγ=1. From (3.13) and (3.2) we obtain
e(cn,APPd,H(KR);Λstd)2≤1n∞∑k=nλd,k+1≤1n∞∑k=1λd,k=trace(Wd)n. | (3.14) |
Since ∑∞j=1γj<∞ and α1>1, we derive from (3.12) that
trace(Wd)≤exp{d∑j=1ln(1+2⌈α1⌉⌈α1⌉ζ(α1)γj)}≤exp{2⌈α1⌉⌈α1⌉ζ(α1)∞∑j=1γj}=:Γ<∞. |
Hence, insertion into (3.14) gives
e(cn,APPd,H(KR);Λstd)2≤Γn. |
From this, we obtain in the same way as above SPT with
τ∗(Λstd)=2=2max(1α1,sγ). |
(ⅱ) Assume that (3.8) holds. This implies that there exists a constant M>0 such that
1lndd∑j=1γj<M,for all d∈N. |
Then, by (3.12) we have
trace(Wd)≤exp{2⌈α1⌉⌈α1⌉ζ(α1)d∑j=1γj}≤exp(2⌈α1⌉⌈α1⌉ζ(α1)lndM)=d2⌈α1⌉⌈α1⌉ζ(α1)M. |
Meanwhile, since the sequence (γj)j∈N is non-increasing, we have
dγdlnd≤1lndd∑j=1γj<M, for all d∈N, |
which implies that* γj=O(j−1lnj) and sγ=1. By Theorem 1 (this implies that APP is SPT for Λall, i.e., there exist two constants Call,pall≥0 such that
Here, Aj=O(Bj) means that there exists a constant C>0 independent of j such that Aj≤CBj.
n(ε,APPd,H(KR);Λall)≤Callε−pall for all ε∈(0,1) and d∈N. |
Now Lemma 2 implies that there is a constant Cstd>0 such that
n(ε,APPd,H(KR);Λstd)≤Cstdε−pstddqstd for all ε∈(0,1) and d∈N, |
where
pstd=pall+2 and qstd=⌈α1⌉⌈α1⌉ζ(α1)M. |
Hence, we have PT and QPT for Λstd.
(ⅲ)–(ⅴ) If any of the three conditions (3.9), (3.10) or (3.11) holds, then this leads to γI<1, since otherwise, for every σ∈(0,1],
limd→∞1dσd∑j=1γj=limd→∞ddσ=limd→∞d1−σ≥1. |
Therefore, we know from Theorem 1 that WT, UWT, and (σ,τ)-WT with σ∈(0,1] hold for Λall.
We only prove the sufficiency of (3.11) for (σ,τ)-WT with σ∈(0,1] since the other two weak tractabiliies can be deduced similarly. Thus, it suffices to show that (3.11) implies (σ,τ)-WT with σ∈(0,1]. Indeed, we observe from (3.12) that
ln(trace(Wd))dσ≤1dσln(d∏j=1(1+2⌈α1⌉⌈α1⌉ζ(α1)γj))=1dσd∑j=1ln(1+2⌈α1⌉⌈α1⌉ζ(α1)γj)≤2⌈α1⌉⌈α1⌉ζ(α1)dσd∑j=1γj, |
where we used that ln(1+x)≤x for all x≥0. Thus (3.11) implies
limd→∞ln(trace(Wd))dσ=0. |
Following the proof of [3, Theorem 26.11], we can obtain (σ,τ)-WT with σ∈(0,1] for Λstd.
The proof is completed.
In this section, we investigate the tractability of integration operators defined on the weighted Korobov space H(KR).} More precisely, we want to approximate the integral operators
INTd:H(KR)→R, INTd(f)=∫[0,1]df(x)dx=ˆf(0), |
in the worst-case setting. We approximate INTd by means of linear algorithms Aintn,d of the form
Aintn,d(f):=n∑i=1λif(xi), | (4.1) |
with nodes x1,…,xn∈[0,1]d and integration weights λ1,…,λn∈R.
The n-th minimal error for integration in H(KR) is defined as
e(n,INTd,H(KR)):=infAintn,dsup‖f‖H(KR)≤1|INTd(f)−Aintn,d(f)| |
where the infimum is taken over all algorithms of the form (4.1). For n=0, the initial integration error is given by
e(0,INTd,H(KR)):=sup‖f‖H(KR)≤1|INTd(f)|=‖INTd‖. |
In the following, we assume that R(0)=1, which leads to the norm of INTd is 1, since
‖INTd‖=sup0≠f∈H(KR)|ˆf(0)|‖f‖H(KR)=sup0≠f∈H(KR)|ˆf(0)|√∑k∈ZdR−1(k)|ˆf(k)|2≤sup0≠f∈H(KR)|ˆf(0)|√|ˆf(0)|2=1 |
and the upper bound can be attained by choosing f=1. Therefore, there is no need to distinguish between ABS and NOR. For abbreviation, we write n(ε,INTd,H(KR))≡nX(ε,INTd,H(KR)).
Obviously, we have
e(n,INTd,H(KR))≤e(n,APPd,H(KR);Λstd), |
which implies
n(ε,INTd,H(KR))≤n(ε,APPd,H(KR);Λstd). | (4.2) |
Now we formulate our result of the tractability of integration INT in the weighted Korobov space H(KR) for R∈{rd,α,γ,ψd,α,γ,ωd,α,γ}.
Theorem 3. Let α={αj}j∈N be the smoothness parameter sequence satisfying (2.4), and let γ={γj}j∈N be the weight sequence satisfying (2.5). Consider the integration problem INT={INTd}d∈N for the weighted spaces H(KR), R∈{rd,α,γ,ψd,α,γ,ωd,α,γ}. Then we have the following conditions:
(i) SPT holds iff
∞∑j=1γj<∞ |
(which is equivalent to sγ≤1). In this case the exponent of SPT satisfies
τ∗(Λstd)=2max(1α1,sγ). |
(ii) PT, QPT are equivalent and hold iff
lim supd→∞1lndd∑j=1γj<∞. |
(iii) WT holds iff
limd→∞1dd∑j=1γj=0. |
(iv) UWT holds iff
limd→∞1dσd∑j=1γj=0, for all σ∈(0,1]. |
(v) For σ∈(0,1] (σ,τ)-WT holds iff
limd→∞1dσd∑j=1γj=0. |
For σ>1, (σ,τ)-WT always holds.
To verify Theorem 3, we shall give the lower bound of the n-th minimal error for INT in H(KR). we will employ the subsequent proposition and lemma.
Proposition 3 ([15]). For j=1,…,d, let βj∈(0,1] and let Hj be a reproducing kernel space on [0,1] such that (1,βjcos(2πx),βjsin(2πx)) are orthonormal in Hj. Consider the integration problem INT for f∈Hd=H1⊗⋯⊗Hd. We have
e(n,INTd,Hd)2≥1−nd∏j=1(1+β2j)−1. |
Lemma 4. Let α={αj}j∈N be the smoothness parameter sequence satisfying (2.4), and let γ={γj}j∈N be the weight sequence satisfying (2.5). Set R∈{rd,α,γ,ψd,α,γ,ωd,α,γ}. Then, for all n∈N,
e(n,INTd,H(KR))≥e(n,INTd,H(Krd,α,γ/3)). |
Proof. This can be deduced from Remark 2, (2.9), and the definition of the n-th minimal error for integration.
Applying Proposition 3 and Lemma 4 to the weighted space H(Krd,α,γ/3), we obtain the lower bound of e(n,INTd,H(KR)) for R∈{rd,α,γ,ψd,α,γ,ωd,α,γ}.
Lemma 5. Let α={αj}j∈N be the smoothness parameter sequence satisfying (2.4), and let γ={γj}j∈N be the weight sequence. Set R∈{rd,α,γ,ψd,α,γ,ωd,α,γ}. Then we have
e(n,INTd,H(KR))2≥1−nd∏j=1(1+13γj)−1. | (4.3) |
This implies that for all ε∈(0,1),
n(ε,INTd,H(KR))≥(1−ε2)d∏j=1(1+13γj). | (4.4) |
Proof. Clearly, the weighted Korobov space H(Krd,α,γ/3) satisfies the conditions of Proposition 3 with βj=√γj/3, j=1,…,d, which, by Proposition 3 and Lemma 4, implies that for R∈{rd,α,γ,ψd,α,γ,ωd,α,γ},
e(n,INTd,H(KR))2≥e(n,INTd,H(Krd,α,γ/3))2≥1−nd∏j=1(1+13γj)−1, |
which gives (4.3). For all ε∈(0,1), letting e(n,INTd,H(KR))≤ε, we will obtain (4.4). This completes the proof.
Remark 4. We remark that for R∈{rd,α,γ,ψd,α,γ}, the constant 1/3 in (4.3) and (4.4) can be replaced by 1. That is,
e(n,INTd,H(KR))2≥1−nd∏j=1(1+γj)−1, |
and for all ε∈(0,1),
n(ε,INTd,H(KR))≥(1−ε2)d∏j=1(1+γj). |
Now we turn to give the proof of Theorem 3.
Proof of Theorem 3. According to (4.2), we know that the sufficient condition on the tractability of the L2-approximation problem for Λstd also works for the integration problem in H(KR). Therefore, we only need to verify the necessary conditions.
First, we assert that limj→∞γj=0. If it does not hold, there would exist a constant γ∗>0 such that γj≥γ∗>0 for every j∈N. According to Lemma 5, it can be derived that
n(ε,INTd,H(KR))≥(1−ε2)(1+13γ∗)d. |
Thus n(ε,INTd,H(KR)) grows exponentially fast in d, which implies that INT suffers from the curse of dimensionality. This proves our claim.
Now we assume that ∑∞j=1γj=∞. Then we can derive from limj→∞γj=0 and the inequality xln2≤ln(1+x)≤x,x∈[0,1] that†
†Here, Ad=Θ(Bd) means that there exsits a constant C>0 independent of d such that C−1Bd≤Ad≤CBd.
d∏j=1(1+13γj)=exp{d∑j=1ln(1+13γj)}=Θ(exp(13d∑j=1γj)). | (4.5) |
Then it follows from (4.4) and (4.5) that
limd→∞n(ε,INTd,H(KR))=∞, |
which implies that INT cannot be SPT. Thus ∑∞j=1γj<∞ is necessary for SPT.
Next assume that lim supd→∞1lnd∑dj=1γj=∞. Then we can derive from limj→∞γj=0 that
d∏j=1(1+13γj)=Θ(d13lndd∑j=1γj). |
Combining with (4.4), we obtain that n(ε,INTd,H(KR)) goes to infinity faster than any power of d, and INT cannot be PT. Thus, lim supd→∞1lnd∑dj=1γj<∞ is necessary for PT.
Finally, assume that
limd+ε−1→∞lnn(ε,INTd,H(KR))dσ+ε−τ=0, σ∈(0,1]. |
Then we can derive from (4.4) and (4.5) that
limd→∞1dσd∑j=1γj=0. |
This implies the condition is necessary for the three WT notions.
The proof is completed.
This paper gives a complete picture of the tractability of multivariate L2-approximation for both Λall and Λstd, and multivariate integration from weighted Korobov spaces of increasing smoothness in the worst-case setting. According to the results in [16], the corresponding tractability conditions of L2-approximation in the randomized case setting are the same as in the worst-case setting. Moreover, our results may be helpful to study the tractability of nonhomogeneous tensor product spaces.
The authors declare they have not used Artificial Intelligence (AI) tools in the creation of this article.
The authors declare there is no conflict of interest.
[1] | E. Novak, H. Woźniakowski, Tractability of Multivariate Problems, Volume I: Linear Information, EMS, Zurich, 2008. https://doi.org/10.4171/026 |
[2] | E. Novak, H. Woźniakowski, Tractability of Multivariate Problems, Volume II: Standard Information for Functionals, EMS, Zurich, 2010. https://doi.org/10.4171/084 |
[3] | E. Novak, H. Woźniakowski, Tractability of Multivariate Problems, Volume III: Standard Information for Operators, EMS, Zurich, 2012. https://doi.org/10.4171/116 |
[4] |
E. Novak, I. H. Sloan, H. Woźniakowski, Tractability of approximation for weighted Korobov spaces on classical and quantum computers, Found. Comput. Math., 4 (2004), 121–156. https://doi.org/10.1007/s10208-002-0074-6 doi: 10.1007/s10208-002-0074-6
![]() |
[5] |
G. W. Wasilkowski, H. Woźniakowski, Weighted tensor product algorithms for linear multivariate problems, J. Complexity, 15 (1999), 402–447. https://doi.org/10.1006/jcom.1999.0512 doi: 10.1006/jcom.1999.0512
![]() |
[6] |
G. W. Wasilkowski, H. Woźniakowski, On the power of standard information for weighted approximation, Found. Comput. Math., 1 (2001), 417–434. https://doi.org/10.1007/s102080010016 doi: 10.1007/s102080010016
![]() |
[7] | H. Woźniakowski, Tractability of multivariate integration for weighted Korobov spaces: My 15 year partnership with Ian Sloan, in Monte Carlo and Quasi-Monte Carlo Methods 2008, (eds. Henryk Woźniakowski), Springer, Berlin, Heidelberg, (2009), 637–653. |
[8] | J. Dick, P. Kritzer, F. Pillichshammer, Lattice Rules: Numerical Integration, Approximation, and Discrepancy, Springer, 2022. https://doi.org/10.1007/978-3-031-09951-9 |
[9] | R. Guo, H. Wang, Tractability of non-homogeneous tensor product problems in the worst case setting, preprint, arXiv: 1901.05752. https://doi.org/10.48550/arXiv.1901.05752 |
[10] |
H. Yan, J. Chen, Tractability of approximation of functions defined over weighted Hilbert spaces, Axioms, 13 (2024). https://doi.org/10.3390/axioms13020108 doi: 10.3390/axioms13020108
![]() |
[11] |
G. Leobacher, F. Pillichshammer, A. Ebert, Tractability of L2-approximation and integration in weighted Hermite spaces of finite smoothness, J. Complexity, 78 (2023) 101768. https://doi.org/10.1016/j.jco.2023.101768 doi: 10.1016/j.jco.2023.101768
![]() |
[12] |
A. Ebert, F. Pillichshammer, Tractability of approximation in the weighted Korobov space in the worst-case setting–a complete picture. J. Complexity, 67 (2021), 15. https://doi.org/10.1016/j.jco.2021.101571 doi: 10.1016/j.jco.2021.101571
![]() |
[13] | J. F. Traub, G. W. Wasilkowski, H. Woźniakowski, Information-Based Complexity, Academic Press, New York, 1988. |
[14] |
P. Kritzer, H. Woźniakowski, Simple characterizations of exponential tractability for linear multivariate problems, J. Complexity, 51 (2019), 110–128. https://doi.org/10.1016/j.jco.2018.10.004 doi: 10.1016/j.jco.2018.10.004
![]() |
[15] |
A. Hinrichs, D. Krieg, E. Novak, J. Vybíral, Lower bounds for the error of quadrature formulas for Hilbert spaces, J. Complexity, 65 (2021), 101544. https://doi.org/10.1016/j.jco.2020.101544 doi: 10.1016/j.jco.2020.101544
![]() |
[16] |
W. Lu, H. Wang, On the power of standard information for tractability for L2-approximation in the randomized setting, Contemp. Math., 3 (2022), 1–29. https://doi.org/10.37256/cm.3120221229 doi: 10.37256/cm.3120221229
![]() |
[17] |
M. Dolbeault, D. Krieg, M. Ullrich, A sharp upper bound for sampling numbers in L2, Appl. Comput. Harmon. Anal., 63 (2023), 113–134. https://doi.org/10.1016/j.acha.2022.12.001 doi: 10.1016/j.acha.2022.12.001
![]() |
rd,α,γ | ψd,α,γ and ωd,α,γ | |
SPT | sγ<∞ | sγ<∞ |
PT | sγ<∞ | sγ<∞ |
QPT | γI<1 | suff.: γI<1 |
UWT | γI<1 | ? |
(σ,τ)-WT, σ∈(0,1] | γI<1 | ? |
WT | γI<1 | ? |
(σ,τ)-WT, σ>1 | no extra condition on γ | no extra condition on γ |
rd,α,γ | ψd,α,γ and ωd,α,γ | |
SPT | δγ>0 | δγ>0 |
PT | δγ>0 | δγ>0 |
QPT | γI<1 | ? |
UWT | γI<1 | ? |
(σ,τ)-WT, σ∈(0,1] | γI<1 | ? |
WT | γI<1 | ? |
(σ,τ)-WT, σ>1 | no extra condition on γ | no extra condition on γ |
Λall | Λstd | |
SPT | sγ<∞ | ∑∞j=1γj<∞ |
PT | sγ<∞ | lim supd→∞1lnd∑dj=1γj<∞ |
QPT | γI<1 | lim supd→∞1lnd∑dj=1γj<∞ |
UWT | γI<1 | limd→∞1dσ∑dj=1γj=0 ∀σ∈(0,1] |
(σ,τ)-WT for σ∈(0,1] | γI<1 | limd→∞1dσ∑dj=1γj=0 |
WT | γI<1 | limd→∞1dd∑j=1γj=0 |
(σ,τ)-WT for σ>1 | no extra condition on γ | no extra condition on γ |
rd,α,γ | ψd,α,γ and ωd,α,γ | |
SPT | sγ<∞ | sγ<∞ |
PT | sγ<∞ | sγ<∞ |
QPT | γI<1 | suff.: γI<1 |
UWT | γI<1 | ? |
(σ,τ)-WT, σ∈(0,1] | γI<1 | ? |
WT | γI<1 | ? |
(σ,τ)-WT, σ>1 | no extra condition on γ | no extra condition on γ |
rd,α,γ | ψd,α,γ and ωd,α,γ | |
SPT | δγ>0 | δγ>0 |
PT | δγ>0 | δγ>0 |
QPT | γI<1 | ? |
UWT | γI<1 | ? |
(σ,τ)-WT, σ∈(0,1] | γI<1 | ? |
WT | γI<1 | ? |
(σ,τ)-WT, σ>1 | no extra condition on γ | no extra condition on γ |
Λall | Λstd | |
SPT | sγ<∞ | ∑∞j=1γj<∞ |
PT | sγ<∞ | lim supd→∞1lnd∑dj=1γj<∞ |
QPT | γI<1 | lim supd→∞1lnd∑dj=1γj<∞ |
UWT | γI<1 | limd→∞1dσ∑dj=1γj=0 ∀σ∈(0,1] |
(σ,τ)-WT for σ∈(0,1] | γI<1 | limd→∞1dσ∑dj=1γj=0 |
WT | γI<1 | limd→∞1dd∑j=1γj=0 |
(σ,τ)-WT for σ>1 | no extra condition on γ | no extra condition on γ |