Preprint
Article

Tractability of Approximation of Functions Defined over Weighted Hilbert Spaces

Altmetrics

Downloads

80

Views

14

Comments

0

A peer-reviewed article of this preprint also exists.

This version is not peer-reviewed

Submitted:

28 December 2023

Posted:

02 January 2024

You are already at the latest version

Alerts
Abstract
We investigate \texorpdfstring{$L_2$}{}-approximation problems in the worst case setting in the weighted Hilbert spaces \texorpdfstring{$H(K_{R_{d,\a,\g}})$}{} with weights \texorpdfstring{$R_{d,\a,{\bm \ga}}$}{} under parameters \texorpdfstring{$1\ge \ga_1\ge \ga_2\ge \cdots\ge 0$ and $1<\az_1\le \az_2\le \cdots$}{}. Several interesting weighted Hilbert spaces \texorpdfstring{$H(K_{R_{d,\a,\g}})$}{} appear in this paper. We consider the worst case error of algorithms that use finitely many arbitrary continuous linear functionals. Multivariate approximation; information complexity; tractability; weighted Hilbert spacesWe discuss tractability of \texorpdfstring{$L_2$}{}-approximation problems for the involved Hilbert spaces, which describes how the information complexity depends on \texorpdfstring{$d$}{} and \texorpdfstring{$\va^{-1}$}{}. As a consequence we study the strongly polynomial tractability (SPT), polynomial tractability (PT), weak tractability (WT), and \texorpdfstring{$(t_1,t_2)$}{}-weak tractability (\texorpdfstring{$(t_1,t_2)$}{}-WT) for all \texorpdfstring{$t_1>1$}{} and \texorpdfstring{$t_2>0$}{} in terms of the introduced weights under the absolute error criterion or the normalized error criterion.
Keywords: 
Subject: Computer Science and Mathematics  -   Mathematics

1. Introduction

We investigate multivariate approximation problems S d with large or even huge d. Examples include these problems in statistics, computational finance and physics. In order to solve these problems we usually consider algorithms using finitely many evaluations of arbitrary continuous linear functionals. We use either the absolute error criterion (ABS) or the normalized error criterion (NOR). For X { ABS , NOR } we define the information complexity n X ( ε , S d ) to be the minimal number of linear functionals which are needed to find an algorithm whose worst case error is at most ε . The behavior of the information complexity n X ( ε , S d ) is the major concern when the accuracy ε of approximation goes to zero and the number d of variables goes to infinity. For small ε and large d, tractability is aimed at studying how the information complexity n X ( ε , S d ) behaves as a function of d and ε 1 , while the exponential convergence-tractability (EC-tractability) is aimed at studying how the information complexity n X ( ε , S d ) behaves as a function of d and ( 1 + ln ( ε 1 ) ) . Recently the study of tractability and EC-tractability in the worst case setting has attracted much interest in analytic Korobov spaces, weighted Korobov spaces and weighted Gaussian ANOVA spaces; see [1,2,3,4,5,6,7,8,9,10,11,12] and the references therein.
Weighted multivariate approximation of functions on space [ 0 , 1 ] d are studied in many problems. We are interested in weighted Hilbert spaces of functions in this paper. We present three examples of weighted Hilbert spaces, which are similar but also different. We devote to discussing worst case tractability of L 2 -approximation problem
APP = APP d : H ( K R d , α , γ ) L 2 ( [ 0 , 1 ] d ) d N
with APP d ( f ) = f for all f H ( K R d , α , γ ) in weighted Hilbert spaces H ( K R d , α , γ ) with three weights R d , α , γ under positive parameter sequences γ = { γ j } j N and α = { α j } j N . The tractability and EC-tractability of such problem APP in weighted Korobov spaces with parameters 1 γ 1 γ 2 0 and 1 < α 1 = α 2 = were discussed in [2,4,6,11] and in [13], respectively. Additionally, [4] considered the tractability of the L 2 -approximation in several weighted Hilbert spaces for permissible information class consisting of arbitrary continuous linear functionals and consisting of functions evaluations.
In this paper we study SPT, PT, WT and ( t 1 , t 2 ) -WT for all t 1 > 1 and t 2 > 0 of the above problem APP with parameters
1 γ 1 γ 2 0 ,
and
1 < α 1 α 2
for the ABS or the NOR under the information class consisting of arbitrary continuous linear functionals. Especially, we get a compete sufficient and necessary condition for SPT, PT and WT, respectively, and the exponent of SPT.
The paper is organized as follows. In Section 2 we give preliminaries about multivariate approximation problems in Hilbert spaces for information class consisting of arbitrary continuous linear functionals in the worst case setting, and definitions of tractability. In Section 3 we present several examples of weighted Hilbert spaces and study some facts and relations between them. In Section 4 we discuss the tractability properties of L 2 -approximation problems in the above weighted Hilbert spaces, then state out main result Theorem 4.1.

2. Approximation and tractability in Hilbert spaces

2.1. Approximation in Hilbert spaces

Let F d and G d be two sequences of Hilbert spaces. Consider a sequence of compact linear operators
S d : F d G d
for all d N . We approximation S d by algorithm A n , d of the form
A n , d ( f ) = i = 1 n T i ( f ) g i , for f S d ,
where functions g i G d and continuous linear functionals T i F d * for i = 1 , , n . The worst case error for the algorithm A n , d of the form (2.1) is defined as
e ( A n , d ) : = sup f F d , | | f | | F d 1 | | S d ( f ) A n , d ( f ) | | G d .
The nth minimal worst-case error, for n 1 , is defined by
e ( n , S d ) : = inf A n , d e ( A n , d ) ,
where the infimum is taken over all linear algorithms of the form (2.1). For n = 0 , we use A 0 , d = 0 . We call
e ( 0 , S d ) = sup f F d , | | f | | F d 1 | | S d ( f ) | | G d
the initial error of the problem S d .
The information complexity for S d can be studied using either the absolute error criterion (ABS), or the normalized error criterion (NOR). The information complexity n X ( ε , S d ) for X { ABS , NOR } is defined by
n X ( ε , S d ) : = min { n N 0 : e ( n , S d ) ε CRI d } ,
where
CRI d : = 1 , for X = ABS , e ( 0 , S d ) , for X = NOR .
Here, N 0 = { 0 , 1 , } and N = { 1 , 2 , } .
It is well known, see e.g., [7,14], that the nth minimal worst case errors e ( n , S d ) and the information complexity n X ( ε , S d ) depend on the eigenvalues of the continuously linear operator W d = S d * S d : F d F d . Let ( λ d , j , η d , j ) be the eigenpairs of W d , i.e.,
W d η d , j = λ d , j η d , j for all j N ,
where the eigenvalues λ d , j are ordered,
λ d , 1 λ d , 2 0 ,
and the eigenvectors η d , j are orthonormal,
η d , i , η d , j F d = δ i , j for all i , j N .
Then the nth minimal error is obtained for the algorithm
A n , d f = j = 1 n f , η d , j F d η d , j for all f F d ,
and
e ( n , S d ) = e ( A n , d ) = λ d , n + 1 for all n N 0 .
Hence the information complexity is equal to
n X ( ε , S d ) = min { n N 0 : λ d , n + 1 ε CRI d } = min { n N 0 : λ d , n + 1 ε 2 CRI d 2 } = | { n N : λ d , n > ε 2 CRI d 2 } | ,
with ε ( 0 , 1 ) and d N . We focus on the rate of the information complexity when the error threshold ε tends to 0 and the problem dimension d grows to infinity.

2.2. Tractability

In order to characterize the dependency of the information complexity n X ( ε , S d ) for the absolute error criterion and the normalized error criterion on the dimension d and the error threshold ε , we will briefly recall some of the basic tractability and exponential convergence-tractability (EC-tractability) notions.
Let S = { S d } d N . For X { ABS , NOR } , we say S is
• strongly polynomially tractable (SPT) iff there exist non-negative numbers C and p such that for all d N , ε ( 0 , 1 ) ,
n X ( ε , S d ) C ( ε 1 ) p .
The exponent p str of SPT is defined to be the infimum of all p for which the above inequality holds.
• polynomially tractable (PT) iff there exist non-negative numbers C, p and q such that for all d N , ε ( 0 , 1 ) ,
n X ( ε , S d ) C d q ( ε 1 ) p .
• quasi-polynomially tractable (QPT) iff there exist two constants C , t > 0 such that for all d N , ε ( 0 , 1 ) ,
n X ( ε , S d ) C exp ( t ( 1 + ln ε 1 ) ( 1 + ln d ) ) .
The exponent t pol of QPT is defined to be the infimum of all t for which the above inequality holds.
• uniformly weakly tractable (UWT) iff for all t 1 , t 2 > 0 ,
lim ε 1 + d ln n X ( ε , S d ) d t 1 + ( ε 1 ) t 2 = 0 ;
• weakly tractable (WT) iff
lim ε 1 + d ln n X ( ε , S d ) d + ε 1 = 0 .
( t 1 , t 2 ) -weakly tractable ( ( t 1 , t 2 ) -WT) for fixed positive t 1 and t 2 iff
lim ε 1 + d ln n X ( ε , S d ) d t 1 + ( ε 1 ) t 2 = 0 .
We call that S suffers from the curse of dimensionality if there exist positive numbers C 1 , C 2 , ε 0 such that for all 0 < ε ε 0 and infinitely many d N ,
n ( ε , d ) C 1 ( 1 + C 2 ) d .
• Exponential convergence-strongly polynomially tractable (EC-SPT) iff there exist non-negative numbers C and p such that for all d N , ε ( 0 , 1 ) ,
n X ( ε , S d ) C 1 + ln ( ε 1 ) p .
The exponent of SPT is defined to be the infimum of all p for which the above inequality holds.
• Exponential convergence-polynomially tractable (EC-PT) iff there exist non-negative numbers C, p and q such that for all d N , ε ( 0 , 1 ) ,
n X ( ε , S d ) C d q 1 + ln ( ε 1 ) p .
• Exponential convergence-uniformly weakly tractable (EC-UWT) iff for all t 1 , t 2 > 0
lim ε 1 + d ln n X ( ε , S d ) d t 1 + 1 + ln ( ε 1 ) t 2 = 0 .
• Exponential convergence-weakly tractable (EC-WT) iff
lim ε 1 + d ln n X ( ε , APP d ) d + ln ( ε 1 ) = 0 .
• Exponential convergence- ( t 1 , t 2 ) -weakly tractable (EC- ( t 1 , t 2 ) -WT) for fixed positive t 1 and t 2 iff
lim ε 1 + d ln n X ( ε , S d ) d t 1 + 1 + ln ( ε 1 ) t 2 = 0 .
Clearly, (1,1)-WT is the same as WT, and EC-(1,1)-WT is the same as EC-WT. Obviously, in the definitions of SPT, PT, QPT, UWT, WT and ( t 1 , t 2 ) -WT, if we replace ε 1 by ( 1 + ln ( ε 1 ) ) , we get the definitions of EC-SPT, EC-PT, EC-QPT, EC-UWT, EC-WT and EC- ( t 1 , t 2 ) -WT, respectively. We also have
SPT PT QPT UWT WT ,
EC - SPT EC - PT EC - QPT EC - UWT EC - WT ,
EC - SPT SPT , EC - PT PT , EC - QPT QPT ,
and
EC - ( t 1 , t 2 ) - WT ( t 1 , t 2 ) - WT , EC - UWT UWT , EC - WT WT .
We can learn more information about tractability of multivariate problems in the volumes [7,8,9] by Novak and Woźniakowski.
Lemma 2.1. 
([7] Theorem 5.2) Consider the non-zero problem S = { S d } for compact linear problems S d defined over Hilbert spaces. Then S is PT for NOR iff there exist q 0 and τ > 0 such that
C τ , q : = sup d N j = 1 ( λ d , j λ d , 1 ) τ 1 τ d q < .
Expecially, S is SPT for NOR iff (2.3) holds with q=0. The exponent of SPT is
p str = inf { 2 τ | τ satisfies ( 2.3 ) with q = 0 } .

3. Weighted Hilbert spaces

Let the space H ( K R d , α , γ ) with weight R d , α , γ under positive parameter sequences γ = { γ j } j N and α = { α j } j N satisfying
1 γ 1 γ 2 0 ,
and
1 < α 1 α 2 .
be a reproducing kernel Hilbert space. The reproducing kernel function K R d , α , γ : [ 0 , 1 ] d × [ 0 , 1 ] d C of the space H ( K R d , α , γ ) is given by
K R d , α , γ ( x , y ) : = k = 1 d K R α k , γ k ( x k , y k ) ,
x = ( x 1 , x 2 , , x d ) , y = ( y 1 , y 2 , , y d ) [ 0 , 1 ] d , where
K R α , γ ( x , y ) : = k N 0 R α , γ ( k ) exp ( 2 π i k · ( x y ) ) , x , y [ 0 , 1 ]
is a universal weighted function. Here Fourier weight R α , γ : N 0 R + be a summable function, i.e., k N 0 R α , γ ( k ) < . We will consider weight R α , γ later on in some examples.
Then we have
K d , α , γ ( x , y ) = k N 0 d R d , α , γ ( k ) exp ( 2 π i k · ( x y ) ) , x , y [ 0 , 1 ] d ,
and the corresponding inner product
f , g H ( K R d , α , γ ) = k N 0 d 1 R d , α , γ ( k ) f ^ ( k ) g ^ ¯ ( k )
and
| | f | | H ( K R d , α , γ ) = f , f H ( K R d , α , γ ) ,
where
R d , α , γ ( k ) : = j = 1 d R α j , γ j ( k j ) , k = ( k 1 , k 2 , , k d ) N 0 d ,
x · y : = k = 1 d x k · y k , x = ( x 1 , x 2 , , x d ) , y = ( y 1 , y 2 , , y d ) [ 0 , 1 ] d ,
and
f ^ ( k ) = [ 0 , 1 ] d f ( x ) exp ( 2 π i k · x ) d x .
We note that the kernel K d , α , γ ( x , y ) is well defined for 1 < α 1 α 2 and for all x , y [ 0 , 1 ] d , since | K d , α , γ ( x , y ) | k N 0 d R d , α , γ ( k ) = j = 1 d ( k N 0 R α j , γ j ( k ) ) < . If γ 1 = γ 2 = = 1 and α 1 = α 2 = > 1 then the space is called unweighted space.
The weights are introduced to model the importance of the functions from the space. The idea can be seen in the reference [15] by Sloan and Woźniakowski. There are various ways to introduce weighted Hilbert spaces. We consider possible choices for Fourier weights R d , α , γ on three examples.

3.1. A Korobov space

Let α = { α j } j N and γ = { γ j } j N satisfy (3.1) and (3.2), respectively. We are interesting in the weighted Korobov space H ( K R d , α , γ ) defined by Irrgeher and Leobacher (see [16]) with kernel (3.3) and corresponding inner product (3.4), where weight R d , α , γ ( k ) = r d , α , γ ( k ) : = j = 1 d r α j , γ j ( k j ) with
r α , γ ( k ) : = 1 , for k = 0 , γ k α , for k 1 ,
for α > 1 and γ ( 0 , 1 ] . Note that we have r α , γ ( k ) ( 0 , 1 ] for all k N 0 .
The space H ( K R d , α , γ ) : = H ( K r d , α , γ ) is a reproducing kernel Hilbert space with parameter sequences α = { α j } j N and γ = { γ j } j N .

3.2. A first variant of the Korobov space

Let α = { α j } j N and γ = { γ j } j N satisfy (3.1) and (3.2), respectively. We consider the reproducing kernel Hilbert space H ( K R d , α , γ ) with kernel (3.3) and corresponding inner product (3.4) determined by R d , α , γ ( k ) = ψ d , α , γ ( k ) : = j = 1 d ψ α j , γ j ( k j ) with
ψ α , γ ( k ) : = 1 , for k = 0 , γ k ! , for 1 k < α , γ ( k α ) ! k ! , for k α ,
for α > 1 and γ ( 0 , 1 ] .
The following lemma gives the upper bound and the lower bound of the weight ψ α , γ ( k ) , which shows that ψ α , γ ( k ) has the same decay rate as the weight r α , γ ( k ) of the Korobov space H ( K r d , α , γ ) under the same parameter sequences α and γ .
Lemma 3.1. 
For all j , k N we have
r α j , γ j ( k ) ψ α j , γ j ( k ) α j α j r α j , γ j ( k ) .
Proof. 
First for all j , k N we want to prove
ψ α j , γ j ( k ) α j α j r α j , γ j ( k ) .
For 1 k < α j we have
ψ α j , γ j ( k ) = γ j k ! γ j γ j α j k α j .
For k α j we have
ψ α j , γ j ( k ) = γ j ( k α j ) ! k ! = γ j k ( k 1 ) ( k α j + 1 ) γ j ( k α j + 1 ) α j = γ j k α j ( 1 α j 1 k ) α j γ j k α j ( 1 α j 1 α j ) α j = α j α j γ j k α j .
We find for all k N that
ψ α j , γ j ( k ) α j α j γ j k α j = α j α j r α j , γ j ( k ) .
Next, for all j , k N we need to prove
ψ α j , γ j ( k ) r α j , γ j ( k ) .
For 1 k < α j we have
ψ α j , γ j ( k ) = γ j k ! γ j k k γ j k α j .
For k α j we have
ψ α j , γ j ( k ) = γ j ( k α j ) ! k ! = γ j k ( k 1 ) ( k α j + 1 ) γ j k α j .
Hence for all j , k N we obtain
ψ α j , γ j ( k ) γ j k α j = r α j , γ j ( k ) .
This finishes the proof. □

3.3. A second variant of the Korobov space

In [17], the reproducing kernel Hilbert space H ( K R d , α , γ ) was considered with kernel (3.3) and corresponding inner product (3.4). Here R d , α , γ ( k ) = ω d , α , γ ( k ) : = j = 1 d ω α j , γ j ( k j ) was defined as
ω α , γ ( k ) : = 1 + 1 γ l = 1 α θ l ( k ) 1 ,
for α > 1 and γ ( 0 , 1 ] , where
θ l ( k ) : = k ! ( k l ) ! , for k l , 0 , for 0 k < l .
Note that for k N we have
l = 1 α θ l ( k ) 2 k α .
Indeed, for k = 1 we have
l = 1 α θ l ( k ) = 1 2 k α ,
for 2 k α we have
l = 1 α θ l ( k ) = l = 1 k k ! ( k l ) ! = k ! l = 0 k 1 1 l ! k ! l = 0 1 l ! k ! e 2 k k 2 k α ,
and for k > α we have
l = 1 α θ l ( k ) = l = 1 α k ! ( k l ) ! l = 1 α k l = k α + k α k k 1 2 k α .
Lemma 3.2. 
For all j , k N we have
1 3 r α j , γ j ( k ) ω α j , γ j ( k ) α j α j r α j , γ j ( k ) .
Proof. 
First for all j , k N we want to prove
ω α j , γ j ( k ) α j α j r α j , γ j ( k ) .
For 1 k < α j we have
ω α j , γ j ( k ) = 1 + 1 γ j l = 1 α j θ l ( k ) 1 = 1 + 1 γ j l = 1 k θ l ( k ) 1 1 γ j θ k ( k ) 1 = γ j k ! .
For k α j we have
ω α j , γ j ( k ) = 1 + 1 γ j l = 1 α j θ l ( k ) 1 1 γ j θ α j ( k ) 1 = γ j ( k α j ) ! k ! .
Hence for all j , k N we get
ω α j , γ j ( k ) ψ α j , γ j ( k ) ,
and thus by Lemma 3.1
ω α j , γ j ( k ) ψ α j , γ j ( k ) α j α j r α j , γ j ( k )
holds.
Next, for all j , k N we need to prove
ω α j , γ j ( k ) 1 3 r α j , γ j ( k ) .
It follows from (3.5) that for all j , k N we have
ω α j , γ j ( k ) = 1 + 1 γ j l = 1 α j θ l ( k ) 1 1 + 2 k α j γ j 1 3 k α j γ j 1 = 1 3 r α j , γ j ( k ) .
This proof is complete. □
Remark 3.3. 
Set R d , α , γ { r d , α , γ , φ d , α , γ , ω d , α , γ } for all j , k N . From Lemma 3.1 and Lemma 3.2 we have for all j , k N ,
1 3 r α j , γ j ( k ) R α j , γ j ( k ) α j α j r α j , γ j ( k ) .
Note that for all j , k N we have ψ α j , γ j ( k ) ψ α 1 , γ j ( k ) , r α j , γ j ( k ) r α 1 , γ j ( k ) , and ω α j , γ j ( k ) ω α 1 , γ j ( k ) , which means that
R α j , γ j ( k ) R α 1 , γ j ( k ) , for all j , k N .
Combining with (3.6) and (3.7), we conclude
1 3 r α j , γ j ( k ) R α j , γ j ( k ) R α 1 , γ j ( k ) α 1 α 1 r α 1 , γ j ( k ) ,
for all j , k N .
Remark 3.4. 
The weight R d , α , γ are used to describe the importance of the different coordinates for the functions from the space H ( K R d , α , γ ) . According to (3.6) we have the weight ψ d , α , γ and the weight ω d , α , γ have the same decay rate as the weight r d , α , γ of the Korobov space H ( K r d , α , γ ) . Hence the above reproducing kernel Hilbert spaces H ( K r d , α , γ ) , H ( K ψ d , α , γ ) and H ( K ω d , α , γ ) are different but also similar.

4. L 2 -approximation in weighted Hilbert spaces and main results

In this section we consider L 2 -approximation
APP d : H ( K R d , α , γ ) L 2 ( [ 0 , 1 ] d )
with APP d ( f ) = f for all f H ( K R d , α , γ ) in Hilbert spaces H ( K R d , α , γ ) with weights R d , α , γ { r d , α , γ , φ d , α , γ , ω d , α , γ } . It is well known from [6] that this embedding APP d is compact with 1 < α 1 α 2 . The kernel K R d , α , γ ( x , y ) is well defined for α 1 > 1 and for all x , y [ 0 , 1 ] d , since by (3.7)
| K R d , α , γ ( x , y ) | k N d R d , α , γ ( k ) = j = 1 d ( 1 + α 1 α 1 ζ ( α 1 ) γ j ) < ,
where ζ ( · ) is the Riemann zeta function.
In the worst case setting the tractability and EC-tractability of L 2 -approximation problems S d with G d = L 2 ( [ 0 , 1 ] d ) were investigated in analytic Korobov spaces and weighted Korobov spaces; see [1,2,3,6,10,11,12,13]. Additionally, [2,6,11,13] discussed tractability and EC-tractability in weighted Korobov spaces.
From subsection 2.1 the information complexity of APP d depends on the eigenvalues of the operator W d = APP d * APP d : H ( K R d , α , γ ) H ( K R d , α , γ ) . Let ( λ d , j , η d , j ) be the eigenpairs of W d ,
W d η d , j = λ d , j η d , j for all j N ,
where the eigenvalues λ d , j are ordered,
λ d , 1 λ d , 2 0 ,
and the eigenvectors η d , j are orthonormal,
η d , i , η d , j H ( K R d , α , γ ) = δ i , j for all i , j N .
Obviously, we have e ( 0 , APP d ) = 1 (or see [6]). Hence the NOR and the ABS for the problem APP d coincide in the worst case setting. We abbreviate n X ( ε , APP d ) as n ( ε , APP d ) , i.e.,
n ( ε , APP d ) : = n X ( ε , APP d ) .
It is well known that the eigenvalues of the operator W d are R d , α , γ ( k ) with k N d ; see, e.g., [7, p. 215]. Hence by (2.2) we have
n ( ε , APP d ) = | { n N : λ d , n > ε 2 } | = | { k N 0 d : R d , α , γ ( k ) > ε 2 } | = | { k N 0 d : j = 1 d R α j , γ j ( k j ) > ε 2 } | .
Tractability such as SPT, PT, WT, and ( t 1 , t 2 ) -WT for t 1 > 1 , and EC-tractability such as EC-WT and EC- ( t 1 , 1 ) -WT for t 1 < 1 of the above problem APP = { APP d } with the parameter sequences γ = { γ j } j N and α = { α j } j N satisfying
1 γ 1 γ 2 0
and
1 < α = α 1 = α 2 =
have been solved by [2,4,11] and [13], respectively. The following conditions have been obtained therein:
• For R d , α , γ { r d , α , γ , φ d , α , γ , ω d , α , γ } , PT holds iff SPT holds iff
s γ : = inf κ > 0 : j = 1 γ j κ < < ,
and the exponent of SPT is
p str = 2 max s γ , 1 α .
• For R d , α , γ = r d , α , γ , QPT, UWT and WT are equivalent and hold iff
γ I : = inf j N γ j < 1 .
For R d , α , γ { φ d , α , γ , ω d , α , γ } ,
γ I <
implies QPT.
In those cases the exponent of QPT is
t pol : = 2 max 1 α , 1 log γ I 1 , for γ I 0 , 2 α , for γ I = 0 .
• For R d , α , γ { r d , α , γ , φ d , α , γ , ω d , α , γ } and t 1 > 1 , ( t 1 , t 2 ) -WT holds for all 1 γ 1 γ 2 0 .
• For R d , α , γ = r d , α , γ , EC-WT holds iff
lim j γ j = 0 .
• For R d , α , γ = r d , α , γ and t 1 < 1 , EC- ( t 1 , 1 ) -WT holds iff
lim j ln j ln ( γ j 1 ) = 0 .
We will research the worst case tractability of the problem APP with sequences satisfying (3.1) and (3.2).
Theorem 4.1. 
Let the sequences γ = { γ j } j N and α = { α j } j N satisfy (3.1) and (3.2). Consider the L 2 -approximation APP for the weighted Hilbert spaces H R d , α , γ , R d , α , γ { r d , α , γ , φ d , α , γ , ω d , α , γ } . Then we have the following tractability results:
(1) SPT and PT are equivalent and hold iff
δ : = lim inf j ln γ j 1 ln j > 0 .
The exponent of SPT is
p str = 2 max 1 δ , 1 α 1 .
(2) For R d , α , γ = r d , α , γ , WT holds iff
lim j γ j < 1 .
(3) For t 1 > 1 , ( t 1 , t 2 ) -WT holds.
Proof. (1) For the problem APP we have λ d , 1 = 1 . Assume that APP is PT. From Lemma 2.1 there exist q 0 and τ > 0 such that
C τ , q : = sup d N j = 1 λ d , j τ 1 τ d q < .
It follows from
j = 1 λ d , j τ = j = 1 d k = 0 ( R α j , γ j ( k ) ) τ = j = 1 d 1 + k = 1 ( R α j , γ j ( k ) ) τ ,
and (3.8) that
> C τ , q j = 1 λ d , j τ 1 τ d q j = 1 d 1 + k = 1 ( 1 3 r α j , γ j ( k ) ) τ 1 τ d q = j = 1 d 1 + 1 3 τ γ j τ ζ ( α j τ ) 1 τ d q j = 1 d 1 + γ j τ 3 τ 1 τ d q ( 1 + γ d τ 3 τ ) d τ d q .
We conclude that
ln C τ , q + q ln d d τ ln ( 1 + γ d τ 3 τ ) d 2 τ · γ d τ 3 τ ,
where we used ln ( 1 + x ) x 2 for all x [ 0 , 1 ] . We further get
ln ( ln C τ , q + q ln d ) ln d τ ln γ d 1 ln ( 2 τ · 3 τ ) ,
i.e.,
ln γ d 1 ln d ln d ln ( ln C τ , q + q ln d ) ln ( 2 τ · 3 τ ) τ · ln d .
Hence we obtain
δ = lim inf d ln γ d 1 ln d 1 τ > 0 .
Note that if APP is SPT, then it is PT. It implies that if APP is SPT, then (4.3) holds and the exponent
p str 2 max { 1 δ , 1 α 1 } .
On the other hand, assume that (4.1) holds. For an arbitrary ε ( 0 , δ 2 ) , there exists an integer N > 0 such that for all j N we have
ln γ j 1 ln j δ ε .
It means that for all j N
γ j j ( δ ε ) .
Choosing τ = 1 δ 2 ε and noting that δ ε δ 2 ε > 1 , we have
j = N γ j τ j = N j ( δ ε ) τ = j = N j δ ε δ 2 ε < ,
which yields that
j = 1 γ j τ < .
From (3.8) we get
j = 1 λ d , j τ 1 τ d q = j = 1 d 1 + k = 1 ( R α j , γ j ( k ) ) τ 1 τ d q j = 1 d 1 + k = 1 ( α 1 α 1 r α 1 , γ j ( k ) ) τ 1 τ d q = d q · exp ln j = 1 d 1 + α 1 α 1 τ γ j τ ζ ( α 1 τ ) 1 τ = d q · exp 1 τ j = 1 d ln 1 + α 1 α 1 τ γ j τ ζ ( α 1 τ ) d q · exp 1 τ j = 1 d α 1 α 1 τ γ j τ ζ ( α 1 τ ) = d q · exp α 1 α 1 τ ζ ( α 1 τ ) τ · j = 1 d γ j τ d q · exp α 1 α 1 τ ζ ( α 1 τ ) τ · j = 1 γ j τ <
for any q 0 and τ > 1 α 1 . We further get
C τ , q = sup d N j = 1 λ d , j τ 1 τ d q <
for any q 0 and τ > 1 α 1 . It follows from Lemma 2.1 that APP is SPT or PT and the exponent p str 2 τ . Setting ε 0 , we obtain
p str 2 τ 2 max 1 δ , 1 α 1 .
Hence the exponent of SPT is p str = 2 max 1 δ , 1 α 1 .
(2) Let τ > 0 . Due to
n λ d , n τ j = 1 n λ d , j τ j = 1 λ d , j τ ,
we have
λ d , n ( j = 1 λ d , j τ ) 1 τ n 1 τ .
Noting that λ d , n ε 2 holds for
n = j = 1 λ d , j τ ε 2 τ ,
we get
n ( ε , APP d ) = min { n | λ d , n + 1 ε 2 } j = 1 λ d , j τ ε 2 τ 1 + ε 2 τ j = 1 λ d , j τ = 1 + ε 2 τ j = 1 d 1 + k = 1 ( R α j , γ j ( k ) ) τ 2 ε 2 τ j = 1 d 1 + k = 1 ( R α j , γ j ( k ) ) τ ,
where we used (4.2).
Set R d , α , γ = r d , α , γ . Assume that lim j γ j < 1 . Then we have from (4.5) that
ln n ( ε , APP d ) d + ε 1 ln 2 ε 2 τ j = 1 d 1 + k = 1 ( R α j , γ j ( k ) ) τ d + ε 1 = ln 2 ε 2 τ j = 1 d 1 + k = 1 ( r α j , γ j ( k ) ) τ d + ε 1 = ln 2 ε 2 τ j = 1 d 1 + γ j τ ζ ( α j τ ) d + ε 1 ln 2 ε 2 τ j = 1 d 1 + γ j τ ζ ( α 1 τ ) d + ε 1 ln 2 + 2 τ ln ( ε 1 ) ε 1 + j = 1 d ln 1 + γ j τ ζ ( α 1 τ ) ) d ln 2 + 2 τ ln ( ε 1 ) ε 1 + j = 1 d γ j τ ζ ( α 1 τ ) d .
We will consider two cases:
  • Case lim j γ j = 0 : It means that for any δ > 0 there exists a positive integer J = J ( δ ) such that
    γ j < δ for all j J .
Then we conclude from (4.6) that
ln n ( ε , APP d ) d + ε 1 ln 2 + 2 τ ln ( ε 1 ) ε 1 + j = 1 d γ j τ ζ ( α 1 τ ) d ln 2 + 2 τ ln ( ε 1 ) ε 1 + ( J 1 ) ζ ( α 1 τ ) + j = J max ( d , J ) δ τ ζ ( α 1 τ ) d ,
which deduces that
lim d + ε 1 ln n ( ε , APP d ) d + ε 1 δ τ ζ ( α 1 τ ) .
Setting δ 0 , we have lim d + ε 1 ln n ( ε , APP d ) d + ε 1 = 0 . This yields WT.
  • Case lim j γ j ( 0 , 1 ) : Then, for every lim j γ j < γ * < 1 there exists a positive integer J 0 = J ( γ * ) such that
    γ j < γ * for all j J 0 .
We have from (4.6) that
ln n ( ε , APP d ) d + ε 1 ln 2 + 2 τ ln ( ε 1 ) ε 1 + j = 1 d γ j τ ζ ( α 1 τ ) d ln 2 + 2 τ ln ( ε 1 ) ε 1 + j = 1 J 0 1 γ j τ ζ ( α 1 τ ) d + j = J 0 max ( J 0 , d ) γ j τ ζ ( α 1 τ ) d ln 2 + 2 τ ln ( ε 1 ) ε 1 + ( J 0 1 ) ζ ( α 1 τ ) d + j = J 0 max ( J 0 , d ) γ * τ ζ ( α 1 τ ) d ,
which means
lim d + ε 1 ln n ( ε , APP d ) d + ε 1 γ * τ ζ ( α 1 τ ) .
Noting that
ζ ( α ) = 1 + k = 2 1 k α 1 + 1 1 x α d x = 1 + 1 α 1 for all α > 1 ,
and setting τ , we obtain
lim d + ε 1 ln n ( ε , APP d ) d + ε 1 lim τ γ * τ ζ ( α 1 τ ) = 0 .
This implies WT.
On the other hand, it suffices to show that WT yields lim j γ j < 1 . Assume on the contrary that lim j γ j = 1 . It yields that γ j 1 for all j N . It follows that
1 = r d , α , γ ( k ) > ε 2
for all k { 0 , 1 } d . Then we have
n ( ε , APP d ) = | { k N 0 d : r d , α , γ ( k ) > ε 2 } | 2 d .
Hence APP suffers from the curse of dimensionality. We cannot have WT.
(3) Let τ > 0 . Due to (4.5) and (3.8) we have
n ( ε , APP d ) 2 ε 2 τ j = 1 d 1 + k = 1 ( R α j , γ j ( k ) ) τ 2 ε 2 τ j = 1 d 1 + α 1 α 1 τ γ j τ ζ ( α 1 τ ) .
It follows that
ln n ( ε , APP d ) d t 1 + ε t 2 ln 2 ε 2 τ j = 1 d 1 + α 1 α 1 τ γ j τ ζ ( α 1 τ ) d t 1 + ε t 2 = ln 2 + 2 τ ln ( ε 1 ) + j = 1 d ln 1 + α 1 α 1 τ γ j τ ζ ( α 1 τ ) d t 1 + ε t 2 ln 2 + 2 τ ln ( ε 1 ) + j = 1 d α 1 α 1 τ γ j τ ζ ( α 1 τ ) d t 1 + ε t 2 ln 2 + 2 τ ln ( ε 1 ) + α 1 α 1 τ ζ ( α 1 τ ) j = 1 d γ j τ d t 1 + ε t 2 ln 2 + 2 τ ln ( ε 1 ) + α 1 α 1 τ d ζ ( α 1 τ ) d t 1 + ε t 2 ln 2 + 2 τ ln ( ε 1 ) ε t 2 + α 1 α 1 τ d ζ ( α 1 ) d t 1 .
We obtain for all t 1 > 1 and t 2 > 0 ,
lim ε 1 + d ln n ( ε , APP d ) d t 1 + ε t 2 = 0 ,
which means APP is ( t 1 , t 2 ) -WT for all t 1 > 1 and t 2 > 0 .
In this paper we consider the SPT, PT, WT and ( t 1 , t 2 ) -WT for all t 1 < 1 and t 2 > 0 for worst case L 2 -approximation in weighted Hilbert spaces H R d , α , γ with parameters 1 γ 1 γ 2 0 and 1 < α 1 α 2 . We get the matching necessary and sufficient conditions on SPT or PT for R d , α , γ { φ d , α , γ , r d , α , γ , ω d , α , γ } and WT for R d , α , γ = r d , α , γ . In particular, it is ( t 1 , t 2 ) -WT for all t 1 > 1 and t 2 > 0 .

Author Contributions

Conceptualization, Yan, H. and Chen, J.; methodology, Yan, H. and Chen, J.; validation, Yan, H.; formal analysis, Yan, H.; investigation, Chen, J.; resources, Chen, J.; data curation, Chen, J.; writing—original draft preparation, Yan, H.; writing—review and editing, Yan, H.; visualization, Yan, H.; supervision, Yan, H. and Chen, J.; project administration, Yan, H. and Chen, J. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Scientific and Technological Innovation Project of Colleges and Universities in Shanxi Province of found grant number 2022L438, National Natural Science Foundation of China of found grant number 12001342, Basic Youth Research Found Project of Shanxi Datong University of found grant number 2022Q10, and Doctoral Foundation Project of Shanxi Datong University of found grant number 2021-B-17 and found grant number 2019-B-10.

Data Availability Statement

No new data were created or analyzed in this study. Data sharing is not applicable to this article.

Acknowledgments

The authors would like to thank all those for important and very useful comments on this paper. The authors were supported by the Scientific and Technological Innovation Project of Colleges and Universities in Shanxi Province (Project no. 2022L438), National Natural Science Foundation of China (Project no. 12001342), Basic Youth Research Found Project of Shanxi Datong University (Project no. 2022Q10), and Doctoral Foundation Project of Shanxi Datong University (Project no. 2021-B-17 and 2019-B-10).

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

References

  1. Dick, J.; Kritzer, P.; Pillichshammer, F.; Woźniakowski, H. Approximation of analytic functions in Korobov spaces. J. Complex. 2014, 30, 2–28. [Google Scholar] [CrossRef]
  2. Eberta, A.; Pillichshammer, F. Tractability of approximation in the weighted Korobov space in the worst-case setting—a complete picture. J. Complex. 2021, 67, 101571. [Google Scholar] [CrossRef]
  3. Irrgeher, C.; Kritzer, P.; Pillichshammer, F.; Woźniakowski, H. Tractability of multivariate approximation defined over Hilbert spaces with exponential weights. J. Approx. Theory 2016, 207, 301–338. [Google Scholar] [CrossRef]
  4. Leobacher, G.; Pillichshammer, F.; Ebert, A. Tractability of L2-approximation and integration in weighted Hermite spaces of finite smoothness, J. Complex. 2023, 78, 101768. [Google Scholar] [CrossRef]
  5. Liu, Y.; Xu, G. Average case tractability of a multivariate approximation problem. J. Complex 2017, 43, 76–102. [Google Scholar] [CrossRef]
  6. Novak, E.; Sloan, I.H.; Woźniakowski, H. Tractability of approximation for weighted Korobov spaces on classical and quantum computers, Found. Comput. Math. 2004, 4(2), 121–156. [Google Scholar]
  7. Novak, E.; Woźniakowski, H. Tractability of Multivariate Problems, Volume I: Liner Information; Publisher: EMS, Zürich, 2008. [Google Scholar]
  8. Novak, E.; Woźniakowski, H. Tractability of Multivariate Problems, Volume II: Standard Information for Functionals; Publisher: EMS, Zürich, 2010. [Google Scholar]
  9. Novak, E.; Woźniakowski, H. Tractability of Multivariate Problems, Volume III: Standard Information for Operators; Publisher: EMS, Zürich, 2012. [Google Scholar]
  10. Wang, H. A note about EC-(s,t)-weak tractability of multivariate approximation with analytic Korobov kernels, J. Complex. 2019, 55, 101412. [Google Scholar] [CrossRef]
  11. Wasilkowski, G.W.; Woźniakowski, H. Weighted tensor product algorithms for linear multivariate problems. J. Complex. 1999, 15(3), 402–447. [Google Scholar] [CrossRef]
  12. Xu, G. EC-tractability of Lp-approximation in Korobov spaces with exponential weights. J. Approx. Theory 2020, 249, 1–20. [Google Scholar] [CrossRef]
  13. Chen, J. EC-(t1,t2)-tractability of approximation in weighted Korobov spaces in the worst case setting. J. Complex. 2022, 73, 101680. [Google Scholar] [CrossRef]
  14. Traub, J.F.; Wasilkowski, G.W.; Woźniakowski, H. Information-Based Complexity; Publisher: Academic Press, New York, 1988. [Google Scholar]
  15. Sloan, I.H.; Woźniakowski, H. When are quasi-Monte Carlo algorithms efficient for high-dimensional integrals? J. Complex. 1998, 14(1), 1–33. [Google Scholar] [CrossRef]
  16. Irrgeher, C.; Leobacher, G. High-dimensional integration on the R d , weighted Hermite spaces, and orthogonal transforms. J. Complex. 2015, 31, 174–205. [Google Scholar] [CrossRef]
  17. Dick, J.; Irrgeher, C.; Leobacher, G; Pillichshammer, F. On the optimal order of integration in Hermite spaces with finite smoothness, SIAM J. Numer. Anal. 2018, 56 (2) (2018), 684-707.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.
Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.
Prerpints.org logo

Preprints.org is a free preprint server supported by MDPI in Basel, Switzerland.

Subscribe

© 2024 MDPI (Basel, Switzerland) unless otherwise stated