arxiv: v2 [math.na] 15 Aug 2018

Transcription

1 Noname manuscript No. (will be inserted by the editor) Lattice rules in non-periodic subspaces of Sobolev spaces Takashi Goda Kosuke Suzuki Takehito Yoshiki arxiv: v2 [math.na] 5 Aug 28 Received: date / Accepted: date Abstract We investigate quasi-monte Carlo (QMC) integration over the s- dimensional unit cube based on rank- lattice point sets in weighted nonperiodic Sobolev spaces H(K) sob and their subspaces of high order smoothness α >, where γ denotes a set of the weights. A recent paper by Dick, Nuyens and Pillichshammer has studied QMC integration in half-period cosine spaces with smoothness parameter α > /2 consisting of non-periodic smooth functions, denoted by H(K cos ), and also in the sum of half-period cosine spaces and Korobov spaces with common parameter α, denoted by H(K kor+cos ). Motivated by the results shown there, we first study embeddings and norm equivalences on those function spaces. In particular, for an integer α, we provide their corresponding norm-equivalent subspaces of H(K). sob This implies that H(K kor+cos ) is strictly smaller than H(K sob ) as sets for α 2, which solves an open problem by Dick, Nuyens and Pillichshammer. Then we study the worst-case error of tent-transformed lattice rules in H(K2,γ,s sob ) and also the worst-case error of symmetrized lattice rules in an intermediate space between H(K kor+cos ) and H(K). sob We show that the almost optimal rate of convergence can be achieved for both cases, while a weak dependence of the worst-case error bound on the dimension can be obtained for the former case. Mathematics Subject Classification (2) 65C5 65D3 65D32 T. Goda School of Engineering, University of Tokyo, 7-3- Hongo, Bunkyo-ku, Tokyo , Japan goda@frcer.t.u-tokyo.ac.jp K. Suzuki Graduate School of Science, Hiroshima University, -3- Kagamiyama, Higashi-Hiroshima, , Japan. JSPS Research Fellow kosuke-suzuki@hiroshima-u.ac.jp T. Yoshiki Department of Applied Mathematics and Physics, Graduate School of Informatics, Kyoto University, Kyoto , Japan yoshiki.takehito.47x@st.kyoto-u.ac.jp

2 2 Takashi Goda et al. Introduction In this paper we study multivariate integration of smooth functions defined overthe s-dimensional unit cube [,] s. For an integrable function f: [,] s R, we denote the integral of f by I(f) = f(x)dx. [,] s A quasi-monte Carlo (QMC) rule denotes an approximation of I(f) by the average of function evaluations on a finite point set P [,] s : Q P (f) = f(x), P x P where we interpret P as a set in which the multiplicity of elements matters. For a function space V with norm V, the worst-case error of a QMC rule using a point set P [,] s is defined by e wor (P;V) := sup Q P (f) I(f). f V f V The aim of this paper is to construct a good deterministic point set P which makes the worst-case error small for a specific function space. There are two main families for QMC point sets: digital nets and integer lattices. We refer the reader to [7,9,6,9] and the references cited therein for general information on this subject. In this paper we focus on integer lattices, in particular, rank- lattice point sets which are defined as follows: Definition Let N,s N and z {,...,N } s. A rank- lattice point set with N points and generating vector z is defined by {{ nz } } P N,z := n < N, N where {x} = x x denotes the fractional part of x R and is applied component-wise for vectors. A QMC rule using a (rank-) lattice point set is called a (rank-) lattice rule. For weighted Korobov spaces consisting of periodic functions whose Fourier coefficients decay algebraically fast, it is well known that there are good generating vectors such that the corresponding lattice rules achieve the almost optimal rate of convergence of the worst-case error and also hold a good dependence of the worst-case error bound on the dimension s [3,,4,2,22]. Here the weights of function spaces play a role in moderating the relative importance of different variables or groups of variables [2]. It is much less known, however, whether or not there are good lattice rules for function spaces consisting of non-periodic smooth functions. In [22], it was proven that the shift-averaged worst-case error of randomly shifted lattice rules in weighted non-periodic Sobolev spaces of first order smoothness, i.e.,

3 Lattice rules in non-periodic subspaces of Sobolev spaces 3 function spaces which consist of non-periodic functions such that the mixed first partial derivatives are square-integrable, coincides with the worst-case error of (deterministic) lattice rules in Korobov spaces with modified weights. Thereafter it was shown in [3] that there exist good generating vectors such that the randomly shifted and then tent-transformed lattice rules achieve the almostoptimalorderofconvergencen 2+ε,ε >,ofthe shift-averagedworstcase error in weighted non-periodic Sobolev spaces of second order smoothness. The point is, however, that the results shown in these papers rely on a random shift of lattice point sets, so that the algorithm is not completely deterministic. Recently, in order to address this issue, Dick, Nuyens and Pillichshammer [8] introduced so-called weighted half-period cosine spaces which consist of non-periodic smooth functions, and together with a successive paper [2], it has been proved that the worst-case error of (deterministic) tent-transformed lattice rules in those spaces is bounded above by the worst-case error of lattice rules in Korobov spaces with modified weights. This means, there are good deterministic tent-transformed lattice rules which achieve the almost optimal rate of convergence of the worst-case error. Moreover, in [8], the sum of halfperiod cosine space and Korobov space was considered and symmetrized lattice rules were shown to achieve the almost optimal rate of convergence in this function space (but with a stronger dependence of the worst-case error bound on the dimension s). We would emphasize, however, that the smoothness of functions in the half-period cosine space is measured not by the differentiability but by the decay rate of the cosine coefficients of functions. The only known exception is that the half-period cosine space with the smoothness parameter α = and the Sobolev space of first order smoothness coincides. Therefore, the measure of smoothness can be equivalently transformed from the decay rate of the cosine coefficients to the first order differentiability. Also for the sum of halfperiod cosine space and Korobov space, it is unknown whether the smoothness of functions can be interpreted in terms of the differentiability. In fact, the authors of [8] commented We do not know whether H(K sob ) differs from H(K kor+cos ) for α >. Here H(K) sob denotes the weighted non-periodic Sobolev space of α-th order smoothness, and H(K kor+cos ) does the sum of the weighted half-period cosine space and the weighted Korobov space with the common parameter α.. Summary of main findings In the light of the above-mentioned researches, the main contribution of this paper is threefold:. In Section 3 we provide two strict subspaces of H(K) sob which are norm equivalent to H(K cos ) and H(Kkor+cos ), respectively. Here H(K cos ) denotes the weighted half-period cosine space with smoothness parameter α.

4 4 Takashi Goda et al. This implies that H(K kor+cos ) is strictly smaller than H(K sob ) as sets, which solves the above problem. 2. In Section 4. we prove that the worst-case error of tent-transformed lattice rules in H(K2,γ,s sob ) is bounded above by the squared worst-case error of lattice rules in Korobov space with the smoothness parameter and modified weights. 3. In Section 4.2 we consider an intermediate Sobolev space, denoted by H(K sob(odd+α) ), between H(K kor+cos ) and H(K sob ), and show that the worst-case error of symmetrized lattice rules is bounded above by the squared worst-case error of lattice rules in Korobov space with the smoothness parameter α/2 and modified weights. The latter two results imply that there are good generating vectors such that the tent-transformed and symmetrized lattice rules achieve the almost optimal rates of convergence N 2+ε and N α+ε with ε >, respectively, in the corresponding function spaces. In fact, the fast component-by-component algorithm due to [8] is directly available to find such good generating vectors. We note that the number of function evaluations for symmetrized lattice rules grows exponentially in the dimension s, so that the worst-case error bound depends exponentially on s regardless of the weights γ, which does not happen for tent-transformed lattice rules. Whether or not deterministic tent-transformed lattice rules can achieve ) has remained unknown for a while after the work of Hickernell [3]. It can be seen from our first main result that the results in [2,8] cannot reach this question. Our second main result gives an affirmative answer to this question. O(N 2+ε ) convergence in H(K sob 2,γ,s.2 Basic notation Throughout this paper, we denote by Z the set of integers and by N the set of positive integers. We write N = N {}, and : s = {,...,s} for s N. For a vector k Z s and a subset u : s, we write k u = (k j ) and (k u,) = l Z s where l j = k j if j u and l j = otherwise. 2 Reproducing kernel Hilbert spaces In order to study embeddings and norm equivalences of normed function spaces, we introduce several reproducing kernel Hilbert spaces (RKHSs) in this section. Our standard reference on the theory of RKHS is []. 2. Korobov spaces For f: [,] s R and h Z s, the h-th Fourier coefficient of f is defined by ˆf(h) := f(x)e 2πih x dx, [,] s

5 Lattice rules in non-periodic subspaces of Sobolev spaces 5 where the dot product denotes the usual inner product in R s. For a set of weightsγ = (γ u ) u :s,γ u,the reproducingkernelofthe weightedkorobov space H(K) kor with smoothness parameter α R, α > /2, is given by K kor (x,y) := h Z s r (h)e 2πih (x y), where the function r : Z s R is defined by r () = and r (h u,) = γ u h j 2α () for vectors h u (Z\{}) u with a non-empty subset u : s. (Here we note that the smoothness parameter α differs by a factor of 2 from what has been used in some literature, see [,22], so that one needs to transfer the results given in those papers carefully.) In particular, when α is a positive integer, it follows from the Fourier series of the Bernoulli polynomial of degree 2α, denoted by B 2α, that the reproducing kernel reduces to K kor (x,y) = + = + γ u h j Z\{} γ u The inner product in H(K kor ) is given by f,g K kor = h Z s ˆf(h)ĝ(h) r (h) = ˆf()ĝ()+ γ u e 2πihj(xj yj) h j 2α (2π) 2α ( ) α+ (2α)! B 2α( x j y j ). ˆf(hu,)ĝ(h u,) h j 2α, h u (Z\{}) u where, for u : s such that γ u =, we assume ˆf(hu,)ĝ(h u,) h j 2α =, h u (Z\{}) u and formally set / :=. In what follows, we shall make similar assumptions on the projections for other weighted Hilbert spaces without any further notice. It is clear from the form of the inner product that the smoothness parameter α moderates the decay rate of the Fourier coefficients. However, the parameter α is also related to the differentiability but only of periodic functions. For simplicity, let s = and α be a positive integer. Then the squared norm of the space H(Kα,γ,) kor can be simplified to ( )2 f 2 K = f(x)dx + α,γ, kor (2π) 2α γ (f (α) (x)) 2 dx,

6 6 Takashi Goda et al. where f (α) denotes the α-th derivative of f, see for instance [7, Appendix A], and H(Kα,γ, kor ) is equivalent to the set { f: [,] R f (τ) absolutely continuous with f (τ) () = f (τ) () } for τ =,...,α, f (α) L 2 [,]. 2.2 Half-period cosine spaces For f: [,] s R and k N s, the k-th cosine coefficient of f is defined by s f(k) := (2 δ,kj ) /2 cos(πk j x j )dx, [,] s f(x) j= where δ denotes the Kronecker delta function. The reproducing kernel of the weighted half-period cosine space H(K cos ) with the smoothness parameter α R, α > /2, introduced in [8] is given by K(x,y) cos := + r (k u,) 2cos(πk j x j )cos(πk j y j ), k u N u where the function r is defined as in (). The inner product in H(K cos ) is given by f,g K cos = k N s f(k) g(k) r (k). Importantly the space H(K cos ) contains non-periodic functions. For instance, as pointed out in [8], the non-periodic function x 2 = 4 π 2 k: odd cos(πkx) k 2 is included in H(Kα,γ, cos ) for /2 < α < 3/2 but not in H(Kkor α,γ, ) for any α > /2. However let us remark again that the smoothness parameter α measures the decay rate of the cosine coefficients and not the differentiability of functions at this point. 2.3 Sum of Korobov and half-period cosine spaces The sum of the weighted Korobov space and the weighted half-period cosine space, denoted by H(K kor+cos ), also introduced in [8] is defined as follows. The reproducing kernel is defined by K kor+cos (x,y) := + [ γ u +K kor+cos α,, (x j,y j ) ],

7 Lattice rules in non-periodic subspaces of Sobolev spaces 7 where we write Kα,, kor+cos (x,y) = Kkor α,,(x,y)+kα,,(x, cos y) 2 cos(2πk(x y)) + cos(πkx) cos(πky) = + k 2α for x,y [,]. According to the result shown in [, Part I, Section 6], the norm in H(K kor+cos ) is given by f 2 K kor+cos = min f=f kor +f cos 2 s( f kor 2 K kor + f cos 2 K cos where the minimum is taken over all functions f kor H(K) kor and f cos H(K cos ) such that f = f kor +f cos. ), 2.4 Non-periodic Sobolev spaces Let α be a positive integer. The reproducing kernel of the weighted nonperiodic Sobolev space H(K) sob of α-th order smoothness is given by K sob [ (x,y) := + γ u +K sob α,, (x j,y j ) ], where we define K sob α,,(x,y) := + α τ= B τ (x)b τ (y) (τ!) 2 +( ) α+b 2α( x y ). (2α)! Here B τ denotes the Bernoulli polynomial of degree τ. The inner product of H(K sob ) is given by f,g = + K sob γ u [,] v ( v u τ u\v {,...,α } u\v f (τ u\v,α v,) (x)dx :s\v [,] ( s v ) [,] s v g (τ u\v,α v,) (x)dx :s\v where (τ u\v,α v,) denotes the s-dimensional vector whose j-th component is τ j if j u\v, α if j v and otherwise, and then f (τ u\v,α v,) denotes the (τ u\v,α v,)-th partial mixed derivative of f. In the univariate case, the corresponding squared norm is given by ( ) 2 f 2 K = f(x)dx + α,γ, sob γ [ α ( 2 f (x)dx) (τ) + τ= ) dx v (f (α) (x)) 2 dx ],

8 8 Takashi Goda et al. and H(Kα,γ, sob ) is equivalent to { f: [,] R f (τ) absolutely continuous } for τ =,...,α, f (α) L 2 [,]. 2.5 Subspaces of non-periodic Sobolev spaces We also consider the following three subspaces of H(K sob ) in this paper. For s =, a weight γ > and α N, define H(K sob(odd) α,γ, ) := {f H(Kα,γ,) sob f (τ) () = f (τ) () for all odd τ α }, H(K sob(odd+α) α,γ, ) := {f H(Kα,γ,) sob f (τ) () = f (τ) () for all odd τ α 2}, H(K sob(odd bdry) α,γ, ) := {f H(K sob α,γ,) f (τ) () = f (τ) () = for all odd τ α }. These spaces are RKHSs equipped with the inner products inherited from H(Kα,γ,) sob since they are its closed subspaces. For instance, the reproducing kernels of the first and second spaces are given by α K sob(odd) α,γ, (x,y) := +γ and K sob(odd+α) α,γ, (x,y) := +γ τ= τ: odd α τ= τ: odd or τ=α B τ (x)b τ (y) (τ!) 2 B τ (x)b τ (y) (τ!) 2 +( ) α+ γ B 2α( x y ) (2α)! +( ) α+ γ B 2α( x y ), (2α)! respectively. We do not give an explicit formula for the kernel of the third one. For an arbitrary dimension s, a set of the weights γ = (γ u ) u :s, γ u, and α N, we define the respective reproducing kernel by K sob( ) (x,y) := + [ ] γ u +K sob( ) α,, (x j,y j ), for {odd,odd+α,odd bdry}. The following obviously holds: H(K sob(odd bdry),γ,s H(K sob(odd bdry) for odd α >, and H(K sob(odd bdry) ) = H(K sob(odd) ) = H(K sob(odd+α),γ,s,γ,s ) = H(K sob,γ,s ), ) H(K sob(odd) ) = H(K sob(odd+α) ) H(K sob ) ) H(K sob(odd) ) H(K sob(odd+α) ) H(K sob ) forevenα,wherewehaveanexceptionwhenα = 2forwhichH(K sob(odd+α) 2,γ,s ) = H(K2,γ,s sob ) holds.

9 Lattice rules in non-periodic subspaces of Sobolev spaces 9 3 Embeddings and norm equivalences Let H(K ) and H(K 2 ) be RKHSs. We say that H(K ) is continuously embedded in H(K 2 ), which is denoted by H(K ) H(K 2 ), if H(K ) H(K 2 ) holds as sets and there exists a constant C > such that f K2 C f K for all f H(K ). If H(K ) H(K 2 ) and H(K 2 ) H(K ) hold, then we say that H(K ) and H(K 2 ) are norm equivalent and write H(K ) H(K 2 ). If there exists a function f such that f H(K ) and f / H(K 2 ), we write H(K ) H(K 2 ). 3. H(K cos ) and H(K sob(odd bdry) ) This subsection is devoted to prove: Lemma For α N, we have H(K cos ) H(Ksob(odd bdry) ). Proof Itsufficestoprovetheresultforthecases = andaweightγ >.Inthe caseα =,theresultimmediatelyfollowsfromtheequalityh(k sob(odd bdry),γ, ) = H(K,γ, sob )andthenormequivalenceh(kcos,γ, ) H(Ksob,γ, )shownin[8].thus let us consider the case α > in the following. First we prove H(Kα,γ,) cos H(K sob(odd bdry) ). Let α,γ, f(x) = a + k Na k 2cos(πkx) for a real sequence(a k ) k N such that f H(K cos α,γ,), that is, (a k ) k N satisfies f 2 K cos α,γ, = f,f K cos α,γ, = a 2 + γ a k 2 k 2α <. As shown in [8, Section 3.], for an odd integer τ with τ α, we have k N f (τ) (x) = k N( ) (τ+)/2 a k 2(πk) τ sin(πkx), which is pointwise absolutely convergent since ( ) (τ+)/2 a k 2(πk) τ sin(πkx) k N

10 Takashi Goda et al. π τ 2 ( ) /2( ) /2 a k k τ π τ 2 a k 2 k 2α k 2α 2τ k N k N k N π τ 2γζ(2α 2τ) f K cos α,γ, <, where we used Cauchy Schwarz inequality in the second inequality and ζ( ) denotes the zeta function. Hence it follows that f (τ) () = f (τ) () =. Further it is known from [8, Theorem ] that H(Kα,γ, cos ) H(Ksob α,γ, ). Thus we conclude H(Kα,γ, cos ) H(Ksob(odd bdry) α,γ, ). NextweproveH(K sob(odd bdry) α,γ, ) H(Kα,γ, cos ).Letf H(Ksob(odd bdry) We note that such a smooth function can be always represented by a cosine series. For any integer k, applying integration by parts α times while using the boundary condition f (τ) () = f (τ) () = for odd integers τ with τ α gives the k-th cosine coefficient f(k) = 2 f(x)cos(πkx)dx = 2 ( πk) α f (α) (x)cs α (πkx)dx, α,γ, ). where we define cs α (x) := { ( ) (α )/2 sin(x) for odd α, ( ) α/2 cos(x) for even α. Hence f 2 K cos α,γ, = ( f()) 2 + γ k 2α ( f(k)) 2 ( ) 2 = f(x)dx + ( γπ 2α ) 2 f(x)dx + γπ 2α ( 2 ) 2 f (α) (x)cs α (πkx)dx (f (α) (x)) 2 dx (+π 2α ) f 2 K, α,γ, sob where the first inequality follows from Bessel s inequality. Thus we have shown that H(K sob(odd bdry) α,γ, ) H(Kα,γ, cos ). 3.2 H(K kor+cos ) and H(K sob(odd) ) This subsection is devoted to prove: Lemma 2 For α N, we have H(K kor+cos ) H(K sob(odd) ).

11 Lattice rules in non-periodic subspaces of Sobolev spaces Proof It suffices to prove the result for the case s = and a weight γ >. Since the embedding H(Kα,γ, kor+cos ) H(K sob(odd) α,γ, ) can be shown in a way similar to the first part of the proof of Lemma, we omit the proof. We prove the converse embedding H(K sob(odd) α,γ, ) H(Kα,γ, kor+cos ). Recall K sob(odd) α,γ, (x,y) = +γ α τ= τ: odd B τ (x)b τ (y) (τ!) 2 =: +γk (x,y)+γk 2 (x,y). +( ) α+ γ B 2α( x y ) (2α)! This means that H(K sob(odd) α,γ, ) is given by the sum of H(γK ) and H( + γk 2 ). The latter space is nothing but the Korobovspace H(Kα,γ(2π) kor 2α,), see Section 2., and it follows from [8, Theorem ] that H(+γK 2 ) is continuously embedded into H(Kα,γ, kor+cos ). Thus, in order to prove the desired embedding, it suffices to show H(K ) H(Kα,γ, kor+cos ). Following [, Part I, Section 3], H(K ) is a RKHS spanned by {B τ (x) τ is an odd integer with τ α} and in particular is finite dimensional. Thus we only need to show H(K ) H(K kor+cos α,γ, ) as sets. Actually we can prove a bit stronger claim that any finite summation of Bernoulli polynomials of odd degrees is in H(K kor+cos α,γ, ). Let f(x) = q k= a k B 2k+ (x) (2k+)!, where q α/2 is a positive integer and a k R. It suffices to show that f is decomposableasf(x) = g (x)+g 2 (x)withg H(Kα,γ,)andg kor 2 H(Kα,γ,). cos LetV := (( π 2 (2k+) 2 ) l ) q k,l= be anon-singularvandermondematrix.define g 2 H(K cos α,γ,) by q c k g 2 (x) := 2 k= cos(π(2k +)x), where c k R are determined as the unique solution of the equation V(c,...,c q ) = (a,...,a q ). Now let g := f g 2. Then it is straightforward to see that g (τ) () = g(τ) () holds for all integers τ 2q +, and that g is infinitely differentiable, which means g H(Kα,γ, kor ). Thus a desired decomposition f = g +g 2 holds.

12 2 Takashi Goda et al. 3.3 Summary of embedding results Together with [8, Theorem ], here we give a summary on embeddings and norm equivalences of the function spaces introduced in the previous section. Theorem Let α N. For a dimension s and a set of weights γ = (γ u ) u :s, γ u, the following holds true:. For α =, H(K cos 2. For odd α >, 3. For even α,,γ,s) H(K kor+cos ) H(K sob = H(K sob(odd),γ,s H(K cos,γ,s,γ,s) ) = H(K sob(odd+α),γ,s ) = H(K sob(odd bdry),γ,s ). ) H(K sob(odd bdry) ) H(K kor+cos H(K sob(odd) H(K cos ) = H(K sob(odd+α) ) ) H(K sob ). ) H(Ksob(odd bdry) ) H(K kor+cos ) H(K sob(odd) ) H(K sob(odd+α) ) H(K sob wherein H(K sob(odd+α) ) = H(K) sob holds when α = 2. ), For α >, since H(K sob ) H(Ksob(odd) ) is obvious and this theorem shows H(K kor+cos ) H(K sob(odd) ), we see that H(K sob ) is strictly larger than H(K kor+cos ), which solves an open problem in [8]. 4 The worst-case error bounds Herewestudytheworst-caseerroroftent-transformedlatticerulesinH(K2,γ,s) sob and the worst-case error of symmetrized lattice rules in H(K sob(odd+α) ) for even α. Since H(K2,γ,s sob ) and H(Ksob(odd+α) ) are strictly larger than H(K2,γ,s cos ) and H(K kor+cos ), respectively, our results are stronger than those in [2, 8]. As already introduced, B τ denotes the Bernoulli polynomial of degree τ N, and in what follows, we write b τ ( ) = B τ ( )/τ!. Further, we define B τ : R R by extending B τ periodically to R and write b τ ( ) = B τ ( )/τ!. The Fourier series of b τ is given by b τ (x) = (2πi) τ k Z\{} e 2πikx k τ, (2) see [9, Appendix C]. Here the above series converges only conditionally for τ =, whereas it converges pointwise absolutely for τ 2.

13 Lattice rules in non-periodic subspaces of Sobolev spaces 3 Moreover we recall that the squared worst-case error of a QMC rule using a point set P [,] s in a RKHS H(K) with a kernel K: [,] s [,] s R is given by (e wor (P;H(K))) 2 = K(x,y)dxdy 2 K(x,y)dy + [,] P 2s [,] P 2 K(x,y), s x P x,y P see [2, Section 6]. This equality shall be used in the subsequent analysis. (3) 4. Tent-transformed lattice rules in H(K sob 2,γ,s ) The tent transformation φ: [,] [,] is defined by φ(x) := 2x, For a rank- lattice point set P N,z with N N and z {,...,N } s, the tent-transformed rank- lattice point set is simply given by P φ N,z := {φ(x): x P N,z}. where φ is applied componentwise, i.e., we write φ(x) = (φ(x ),...,φ(x s )). A QMC rule using P φ N,z is called a tent-transformed (rank-) lattice rule. Since we have K2,γ,s sob (x,y)dxdy = and K2,γ,s sob (x,y)dy =, [,] s [,] s [,] s for any x [,] s, it follows from (3) that (e wor (P φ N,z ;H(Ksob = + N 2 = N 2 2,γ,s )))2 x,y P N,z K sob x,y P N,z 2,γ,s (φ(x),φ(y)) γ u In Appendix A, we prove the following: Lemma 3 For any x,y [,], we have +K sob 2,, (φ(x),φ(y)) = π 4 [ +K sob 2,,(φ(x j ),φ(y j )) ]. (4) c φ (k,l)cos(2πkx)cos(2πly) k,l=

14 4 Takashi Goda et al. where ( 58 k ) π 2 k 2 k 2 l 2 c φ (k,l) = ( π 2 k 2 6 π 2 l 2 ) k,l: odd,k = l, k,l: odd,k l, 6 k 4 k,l: even,k = l, 4 k 2 l 2 k,l: even,k l, otherwise. In order to prove our main result in this subsection, we need to introduce the definition of the dual lattice of a rank- lattice point set. Definition 2 For N,s N and z {,...,N } s, the dual lattice of the rank- lattice point set P N,z is defined by P N,z := {k Zs k z (mod N)}. For a non-empty subset u : s, we write P,u N,z {k = u (Z\{}) u } (k u,) z (mod N). It is straightforward to see that we have the decomposition PN,z \{} = P,u N,z. Lemma 4 Let N,s N and z {,...,N } s. For k Z s, we have { e 2πik x if k PN,z =, N otherwise. x P N,z Proof See [9, Lemma 2.7]. We can show a bound on the worst-case error of tent-transformed lattice rules in H(K2,γ,s sob ) as follows. Theorem 2 Let N,s N and z {,...,N } s. The squared worst-case error of a tent-transformed rank- lattice rule in H(K2,γ,s sob ) is bounded by 2 (e wor (P φ N,z ;H(Ksob 2,γ,s )))2 γ u k u P,u N,z k 2 j (e wor (P N,z ;H(K kor,γ /2,s )))4, where γ u = γ u(c /(4π 4 )) u with c = 58/3 for u : s, and we write γ /2 = (γ /2 u ) u :s.

15 Lattice rules in non-periodic subspaces of Sobolev spaces 5 Proof Using (4) and Lemma 3, the squared worst-case error is given by (e wor (P φ N,z ;H(Ksob = N 2 = 2,γ,s))) 2 x,y P N,z γ u π 4 u γ u k u,l u N u N 2 [ +K sob 2,, (φ(x j),φ(y j )) ] x,y P N,z c φ (k j,l j )cos(2πk j x j )cos(2πl j y j ). Using the identity cos(2πkx) = (e 2πikx +e 2πikx )/2 and Lemma 4, the innermost sum reduces to N 2 = x,y P N,z 4 u N 2 = 4 u c φ (k j,l j )cos(2πk j x j )cos(2πl j y j ) x,y P N,z σ u,σ u {±} u σ u,σ u {±} u = 4 u σ u,σ u {±} u σ u(k u),σ u (lu) P,u N,z c φ (k j,l j )e 2πσjkjxj e 2πσ j ljyj c φ (k j,l j ) N 2 c φ (k j,l j ), e 2πiσu(ku) xu e 2πiσ u (lu) y u x P N,z y P N,z where we have used the notation σ u (k u ) = (σ j k j ). Noting that c φ (k,l) c /(k 2 l 2 ) for all k,l N, we have (e wor (P φ N,z ;H(Ksob = γ u (4π 4 ) u = 2,γ,s))) 2 k u,l u N u γ u (4π 4 ) u k u,l u P,u N,z ( ) c u γ u 4π 4 σ u,σ u {±} u σ u(k u),σ u (lu) P,u N,z c φ ( k j, l j ) k u,l u P,u N,z k 2 j l2 j c φ (k j,l j ) = γ u k u P,u N,z A further upper bound is obtained by using Jensen s inequality: (e wor (P φ N,z ;H(Ksob 2,γ,s )))2 γ u k u P,u N,z k 2 j 2 k 2 j 2.

16 6 Takashi Goda et al. γ /2 u k u P,u N,z k 2 j = (e wor (P N,z ;H(K kor,γ /2,s )))4, where the last equality is well known,see for instance [, ProofofTheorem 2]. This completes the proof of the theorem. Using the Fourier series of b 2 as shown in (2), the squared worst-case error (e wor (P N,z ;H(K kor,γ /2,s )))2 forarank-latticepointsetp N,z hasthefollowing computable formula (e wor (P N,z ;H(K kor,γ /2,s )))2 = = N = N γ /2 u N k u (Z\{}) u γ u /2 x P N,z x P N,z =u :s 2 x P N,z k j Z\{} (2c /2 ) u γ /2 u In particular, for product weights γ u = γ j, we have (e wor (P N,z ;H(K kor,γ /2,s )))2 = + N s x P N,z j= e 2πikjxj k 2 j b 2 (x j ). e 2πikjxj k 2 j [ ] +2c /2 γ /2 j b 2 (x j ). A worst-case error bound of a rank- lattice rule found by the so-called component-by-component (CBC) algorithm in the weighted Korobov space with general weights was proved, for instance, in [, Theorem 5], which is therefore also applicable in our setting. Corollary Let N be a prime. The CBC algorithm using the squared worstcase error (e wor (P N,z ;H(K kor,γ /2,s )))2 as a quality criterion can find a generating vector z {,...,N } s such that e wor (P φ N,z ;H(Ksob 2,γ,s )) (ewor (P N,z ;H(K kor,γ /2,s )))2 (N ) /λ =u :s γ λ/2 u ( 2 λ c λ/2 ζ(2λ) π 2λ holds for any λ (/2,], where ζ( ) denotes the zeta function. (5) ) u /λ Remark The number of points contained in a tent-transformed lattice point set P φ N,z is at most N. Under certain conditions on the weights, the worst-case error bound depends only polynomially on s or even becomes independent of s. We refer to [3,,4] for such results which also apply in our setting.

17 Lattice rules in non-periodic subspaces of Sobolev spaces 7 Remark 2 The fast CBC algorithm using the fast Fourier transform due to Nuyens and Cools [8] can be used to search for a good generating vector z. For product weights, we need O(sN logn) arithmetic operations and O(N) memory. Moreover, for product and order dependent (POD) weights, γ u = Γ u γ j, we need O(sN logn +s 2 N) arithmetic operations and O(sN) memory [5]. Remark 3 Calculating the Fourier series given in Lemma 3 is the crucial step to obtain the results of this section. We see that the decay of the Fourier coefficients c φ (k,l) O(/(k 2 l 2 )) leads to the almost optimal order of convergence O(N 2+ε ). In principle, a similar calculation can be done for a general α > 2. As can be expected from the proof in Appendix A, however, the decay of the Fourier coefficients remains c φ (k,l) O(/(k 2 l 2 )), which means that we cannot improve the order of convergence. 4.2 Symmetrized lattice rules in H(K sob(odd+α) ) For u : s and x = (x,...,x s ) [,] s, we define sym u (x) = (y,...,y s ) [,] s where y j := x j if j u and y j := x j otherwise. For a rank- lattice point set P N,z, the symmetrization of P N,z is defined as the multiset P sym N,z := sym u (P N,z ). u :s A QMC rule using P sym N,z is called a symmetrized (rank-) lattice rule. The worst-case error of symmetrized lattice rules in H(K sob(odd+α) ) is bounded as follows. We give the proof in Appendix B. Theorem 3 Let N,s N and z {,...,N } s. The squared worst-case error of a symmetrized rank- lattice rule in H(K sob(odd+α) ) for even α is bounded by 2 (e wor (P sym ))) 2 N,z ;H(Ksob(odd+α) γ u k u P,u N,z k α j (e wor (P N,z ;H(K kor α/2,γ /2,s )))4, where γ u = γ u (3/(2 2α+ π 2α )) u for u : s and γ /2 = (γ /2 u ) u :s. Following a similar argument as in Section 4., the squared worst-case error (e wor (P N,z ;H(K kor α/2,γ /2,s )))2 can be shown to have the computable formulas: (e wor (P N,z ;H(K kor α/2,γ /2,s )))2

18 8 Takashi Goda et al. = N x P N,z =u :s for general weights γ, and for product weights. ( 3 ( ) +α/2 2 ) u γ /2 u b α (x j ), (e wor (P N,z ;H(K kor α/2,γ /2,s [ )))2 = + ] s 3 + N ( ) +α/2 2 γ/2 j b α (x j ), x P N,z j= Corollary 2 Let N be a prime. The fast CBC algorithm using the squared worst-case error (e wor (P N,z ;H(K kor α/2,γ /2,s )))2 as a quality criterion can find a generating vector z {,...,N } s such that e wor (P sym N,z ;H(Ksob(odd+α) )) (e wor (P N,z ;H(K kor α/2,γ /2,s )))2 ( γ λ/2 3 λ/2 ) u ζ(αλ) (N ) /λ u 2 (2α+)λ/2 π αλ holds for any λ (/α,]. =u :s Remark 4 As proven in [8, Lemma 2], the number of points contained in a symmetrized lattice point set P sym N,z is 2s (N +) for odd N and 2 s N + forevenn.thus,intermsofthenumberofpoints,theexponentialdependence of the worst-case error bound on s cannot be avoided no matter how fast the weights decay. This is a major drawback of symmetrized lattice rules. /λ 5 Numerical results First let us consider a simple bi-variate test function f(x,y) = yexy e 2, whose integral over [,] 2 is. We approximate I(f) by using lattice rules, tent-transformed lattice rules and symmetrized lattice rules. For all of the three rules, Fibonacci lattice point sets are employed as the underlying nodes. Note that f H(Kα,γ,2 sob ) for any α N, but f / H(Ksob(odd+α) α,γ,2 ) for α > 2 since it does not satisfy the periodic condition on high order derivatives. Further since H(K sob(odd+α) 2,γ,2 ) = H(K2,γ,2 sob ) holds, both tent-transformed lattice rules and symmetrized lattice rules areexpected to achieveo(n 2+ε ) convergence.this can be confirmed in Figure. The integration error for symmetrized lattice

19 Lattice rules in non-periodic subspaces of Sobolev spaces 9-2 Absolute error log N 2 Fig. The results for f by using lattice rules ( ), tent-transformed lattice rules ( ) and symmetrized lattice rules (+). The dotted lines correspond to O(N ) and O(N 2 ), respectively. rules convergesexactly with order N 2, whereas the integration errorfor tenttransformed ones behaves periodically but still converges approximately with order N 2. Although a detailed analysis on the periodic behavior is beyond the scope of this paper, this periodicity coincides with whether the number of the underlying nodes is even or odd. Let us move onto the high-dimensional setting. We consider the following test integrands f (x) = f 2 (x) = s j= [+ω j ( x c j [ s + j= ω j +ω j x c2 j +c ], )], for the non-negative parameters c,c 2 and ω j. These smooth integrands were originally used in [6]. Note that the exact values of the integrals for f and for f 2 with the special cases c 2 = and c 2 = 2 are known: I(f ) = and s [+log(+ω j )] for c 2 =, j= I(f 2 ) = s [ + ωj tan ( ω j ) ] for c 2 = 2. j= We put s = 2 and ω j = /j 2 for the integrands. Since the problem is high-dimensional, we focus on tent-transformed lattice rules. We search for a generating vector using the fast CBC algorithm based on the quality criterion (5) with the choice of the weights γ j = /(4c j 2 ). We compare the performance of tent-transformed lattice rules with that of order 2 Sobol sequences

20 2 Takashi Goda et al Absolute error Absolute error log N log N Absolute error -6-7 Absolute error log N log N Absolute error Absolute error log N log N 2 Fig. 2 The results for f with c =.3 (top), f 2 with c 2 = (middle) and f 2 with c 2 = 2 (bottom) by using order 2 Sobol sequences (dashed) and tent-transformed lattice rules (solid). The left column shows the results for s = 2 and ω j = /j 2, while the right column does for s = and ω j = /j 3. In each graph, the dotted line corresponds to O(N 2 ) convergence. as implemented in [5]. Here higher order digital nets and sequences introduced in [4] are known to achieve high order convergence for non-periodic smooth integrands. A recent improvement even shows that those nets and sequences achieve the optimal order of the error convergence [], so that we employ order 2 Sobol sequences as a competitor for present numerical experiments. The absolute integration errors as functions of log 2 N are shown in the left column of Figure 2. For f 2 with c 2 =, order 2 Sobol sequences perform slightly better than tent-transformed lattice rules. For f and f 2 with c 2 = 2,

21 Lattice rules in non-periodic subspaces of Sobolev spaces 2 however, tent-transformed lattice rules outperform order 2 Sobol sequences for many N. For any case, the integration error converges approximately with order N 2, which supports our theoretical result. Finally we put s = and ω j = /j 3 for the same integrands. We search for a generating vector using the fast CBC algorithm based on (5) with γ j = /(4c j 3 ). The absolute integration errors as functions of log 2 N are shown in the right column of Figure 2, where we see that tent-transformed lattice rules compare well with order 2 Sobol sequences. Acknowledgements This work was supported by JSPS Grant-in-Aid for Young Scientists No. 5K2964 (T. G.), JSPS Grant-in-Aid for JSPS Fellows No. 7J466 (K. S.) and No. 7J265 (T. Y.), and JST CREST. A Proof of Lemma 3 In order to prove Lemma 3, we need the following result. Lemma 5 For k N, we have sin(2πkφ(x)) = 8 π k 4k 2 l 2 cos(2πlx). l= l: odd Proof Let l N. Since φ is given by { 2x if x [,/2], φ(x) = 2 2x otherwise, it is an easy exercise to check that and sin(2πkφ(x)) sin(2πlx) dx =, sin(2πkφ(x)) cos(2πlx) dx = sin(2πkx) cos(πlx) dx = (sin((2k +l)πx)+sin((2k l)πx)) dx 2 for even l, = 4 k π 4k 2 l 2 for odd l. Thus the Fourier series of sin(2πkφ( )) is given by sin(2πkφ(x)) = 2 cos(2πlx) sin(2πkφ(y)) cos(2πly) dy l= = 8 k π 4k 2 l 2 cos(2πlx). l= l: odd Hence we complete the proof.

22 22 Takashi Goda et al. Now we are ready to prove Lemma 3. Proof of Lemma 3 By definition, we have +K sob 2,, (φ(x),φ(y)) = b (φ(x))b (φ(y))+b 2 (φ(x))b 2 (φ(y)) b 4 (φ(x) φ(y)). (6) Using the Fourier series of a triangle wave provided in [2, Chapter,.444], we have b (φ(x)) = φ(x) 2 = 4 π 2 k: odd Thus the Fourier series of the first term of (6) is given by b (φ(x))b (φ(y)) = π 4 k,l= k,l: odd cos(2πkx) k 2. 6 k 2 l 2 cos(2πkx)cos(2πly). Using the Fourier series of b 2 as shown in (2) and the equality cos(2πkφ(x)) = cos(4πkx) which holds for any k N and x [,], we have b 2 (φ(x)) = 2π 2 cos(2πkφ(x)) k 2 = 2 π 2 Thus the Fourier series of the second term is given by b 2 (φ(x))b 2 (φ(y)) = π 4 k,l= k,l: even k: even cos(2πkx) k 2. 4 k 2 l 2 cos(2πkx)cos(2πly). Finally let us consider the third term of (6). Using the Fourier series of b 4 and the equality cos(2πkφ(x)) = cos(4πkx), we have b4 (φ(x) φ(y)) = 2 (2π) 4 = 2 π 4 k: even cos(2πk(φ(x) φ(y))) k 4 cos(2πkx) cos(2πky) k 4 8π 4 Using Lemma 5, the second term of the last expression is given by 8π 4 = 8 π 6 sin(2πkφ(x)) sin(2πkφ(y)) k 4 k 2 l,m= l,m: odd = 8 π 6 cos(2πlx) cos(2πmy) l,m= l,m: odd sin(2πkφ(x)) sin(2πkφ(y)) k 4. (4k 2 l 2 )(4k 2 m 2 ) cos(2πlx)cos(2πmy) k 2 (4k 2 l 2 )(4k 2 m 2 ). (7) Noting that the equalities 4k 2 l 2 = 2l 2 and (4k 2 l 2 ) 2 = ( π 2 4l ) l 2.

23 Lattice rules in non-periodic subspaces of Sobolev spaces 23 hold for any positive odd integer l, the inner sum of (7) can be evaluated as follows. In case of l = m, we have k 2 (4k 2 l 2 ) 2 = ( ( l 2 l 2 k 2 4 Otherwise if l m, we have = l 2 ( l 2 ( π 2 k 2 (4k 2 l 2 )(4k 2 m 2 ) = l 2 m 2 = l 2 m 2 ( π l 2 ( ( k m 2 l 2 m 2 ( m l 2 m 2 ) + 4k 2 l 2 )+ l ( π 2 2 2l 2 l2 2m 2 ) 4 (4k 2 l 2 ) 2 )) = l l 2 4k 2 l 2 l 2 )) 4k 2 m 2 )) = ( π 2 l 2 m l 2 2 m 2 ( 5 2 π2 4 ) l 2. By substituting these results on the Fourier series into (6), the result of the lemma follows. ). B Proof of Theorem 3 In what follows, for a function f defined over [,] 2, we define a function sym[f] by Since we have [,] s sym[f(x,y)] := K sob(odd+α) [,] s f(x,y)+f( x,y)+f(x, y)+f( x, y). 4 (x,y)dxdy = and for any x [,] s, it follows from (3) that (e wor (P sym N,z ;H(Ksob(odd+α) ))) 2 = + = + N 2 = N 2 (2 s N) 2 x,y P sym N,z 2 2s v,w :s γ u x,y P N,z x,y P N,z u :s K sob(odd+α) (x, y) As we assume that α is even, we have K sob(odd+α) (x,y)dy =, [,] s K sob(odd+α) (sym v (x),sym w (y)) [ +sym[k sob(odd+α) α,, (x j,y j )] ]. (8) +sym[k sob(odd+α) α,, (x, y)] α = sym b τ(x)b τ(y)+b α(x)b α(y) b 2α (x y) τ= τ : odd α ] = sym[b τ(x)b τ(y)] +sym[b α(x)b α(y)] sym[ b2α (x y). (9) τ= τ: odd

24 24 Takashi Goda et al. Since b τ(x) = b τ( x) for odd τ and b τ(x) = b τ( x) for even τ, we have sym[b τ(x)b τ(y)] = { for odd τ, b τ(x)b τ(y) for even τ. and the first term of (9) equals. Considering the Fourier series of b α for even α: b α(x) = ( )+α/2 (2π) α k Z\{} e 2πikx the Fourier series of the second term of (9) is given by sym[b α(x)b α(y)] = b α(x)b α(y) = k α 4 (2π) 2α = 2 ( )+α/2 (2π) α cos(2πkx) k α, k α l α cos(2πkx)cos(2πly). k,l= Finally, using the Fourier series of b 2α, the Fourier series of the third term of (9) is given by ] sym[ b2α (x y) = 2 (2π) 2α = 2 (2π) 2α sym[cos(2πk(x y))] k 2α cos(2πkx) cos(2πky) k 2α. By substituting these results on the Fourier series into (9), we have +sym[k sob(odd+α) α,, (x, y)] = (2π) 2α c sym (k,l)cos(2πkx)cos(2πly), k,l= where 6 c sym (k,l) = k 2α k = l, 4 k α l α k l. The rest of the proof follows exactly in the same manner as the proof of Theorem 2. References. N. Aronszajn, Theory of reproducing kernels, Trans. Am. Math. Soc. 68 (95), R. Cools, F. Y. Kuo, D. Nuyens and G. Suryanarayana, Tent-transformed lattice rules for integration and approximation of multivariate non-periodic functions, J. Complexity 36 (26), J. Dick, On the convergence rate of the component-by-component construction of good lattice rules, J. Complexity 2 (24), J. Dick, Walsh spaces containing smooth functions and quasi-monte Carlo rules of arbitrary high order, SIAM J. Numer. Anal. 46 (28), J. Dick, How to generate higher order Sobol points in Matlab and some numerical examples. (2). Accessed 24 July J. Dick, T. Goda and T. Yoshiki, Richardson extrapolation of polynomial lattice rules, ArXiv preprint, arxiv: J. Dick, F. Y. Kuo and I. H. Sloan, High-dimensional integration: the quasi-monte Carlo way, Acta Numer. 22 (23), J. Dick, D. Nuyens and F. Pillichshammer, Lattice rules for nonperiodic smooth integrands, Numer. Math. 26 (24),

25 Lattice rules in non-periodic subspaces of Sobolev spaces J. Dick and F. Pillichshammer, Digital Nets and Sequences: Discrepancy Theory and Quasi-Monte Carlo Integration, Cambridge University Press, 2.. J. Dick, I. H. Sloan, X. Wang and H. Woźniakowski, Good lattice rules in weighted Korobov spaces with general weights, Numer. Math. 3 (26), T. Goda, K. Suzuki and T. Yoshiki, Optimal order quadrature error bounds for infinitedimensional higher-order digital sequences, Found. Comput. Math. 8 (28), I. S. Gradshteyn and I. M. Ryzhik, Table of Integrals, Series, and Products, Academic Press, F. J. Hickernell, Obtaining O(N 2+ǫ ) convergence for lattice quadrature rules, In: K. T. Fang, F. J. Hickernell and H. Niederreiter, eds., Monte Carlo and Quasi-Monte Carlo Methods 2, Springer, Berlin, 22, pp F. Y. Kuo, Component-by-component constructions achieve the optimal rate of convergence for multivariate integration in weighted Korobov and Sobolev spaces, J. Complexity 9 (23), F. Y. Kuo, Ch. Schwab and I. H. Sloan, Quasi-Monte Carlo methods for highdimensional integration: the standard (weighted Hilbert space) setting and beyond, ANZIAM J. 53 (2), 37 (corrected in 54 (23) 26 29). 6. H. Niederreiter, Random Number Generation and Quasi-Monte Carlo Methods, CBMS- NSF Regional Conference Series in Applied Mathematics, Vol. 63, SIAM, Philadelphia, E. Novak and H. Woźniakowski, Tractability of multivariate problems. Volume I: Linear information, EMS Tracts in Mathematics, Vol. 6, European Mathematical Society (EMS), Zürich, D. Nuyens and R. Cools, Fast algorithms for component-by-component construction of rank- lattice rules in shift-invariant reproducing kernel Hilbert spaces, Math. Comp. 75 (26), I. H. Sloan and S. Joe, Lattice Methods for Multivariate Integration, Oxford University Press, I. H. Sloan and A. V. Reztsov, Component-by-component construction of good lattice rules, Math. Comp. 7 (22), I. H. Sloan and H. Woźniakowski, When are quasi-monte Carlo algorithms efficient for high-dimensional integrals? J. Complexity 4 (998), I. H. Sloan and H. Woźniakowski, Tractability of multivariate integration for weighted Korobov classes, J. Complexity 7 (2),