Construction Algorithms for Higher Order Polynomial Lattice Rules
|
|
- August Norton
- 6 years ago
- Views:
Transcription
1 Construction Algorithms for Higher Order Polynomial Lattice Rules Jan Baldeaux, Josef Dick, Julia Greslehner and Friedrich Pillichshammer Dedicated to Gerhard Larcher on the occasion of his 50th birthday Abstract Higher order polynomial lattice point sets are special types of digital higher order nets which are known to achieve almost optimal convergence rates when used in a quasi-monte Carlo algorithm to approximate high-dimensional integrals over the unit cube. Recently it has been shown that higher order polynomial lattice point sets of good quality must exist. However, it was not shown how to construct such point sets avoiding an exhaustive search. This is the contribution of the present paper. We use a component-by-component approach to construct higher order polynomial lattice rules achieving optimal convergence rates for functions of arbitrarily high smoothness and at the same time under certain conditions on the weights strong polynomial tractability. In addition, we show how to combine a sieve-type algorithm with the component-by-component approach to construct higher order polynomial lattice rules adjusting themselves to the smoothness of the The support of the Australian Research Council under its Centre of Excellence Program is gratefully acknowledged. The support of the Australian Research Council under its Centre of Excellence Program is gratefully acknowledged. The author is supported by an Australian Research Council Queen Elizabeth II Research Fellowship. F.P. is supported by the Austrian Science Foundation FWF, Project S9609, that is part of the Austrian National Research Network Analytic Combinatorics and Probabilistic Number Theory.
2 2 integrand up to a certain given degree. Analogous results for higher order Korobov polynomial lattice point sets are presented as well. Introduction Quasi-Monte Carlo rules are equal weight integration formulas used to approximate integrals over the unit cube [0, ] s, where the dimension s is typically large. In particular, one approximates an integral [0,] s fx dx by Q N,s f = N N n=0 fx n where x 0,..., x N [0, s. Popular choices for the underlying integration nodes x 0,..., x N [0, s are either lattice point sets see [4, 5] or digital t, m, s-nets see [2, 4]; in this paper, we focus on digital nets. Recently, digital higher order nets were introduced by Dick [3] which include digital t, m, s-nets as special cases and have the appealing property that they can exploit the smoothness of the integrand under consideration. This is not possible with ordinary digital t, m, s-nets. To be more precise, if the integrand under consideration has square integrable mixed partial derivatives of order α in each variable, then digital higher order nets consisting of N points can produce integration errors converging at a rate of N α+ε with arbitrary small ε > 0. Having established the desirable properties of digital higher order nets, an important question is how to construct them. One possible answer to this question was given in [3, Section 4.4], where it was shown how to obtain digital higher order nets from classical digital nets. The aim of the present paper is to provide constructions of digital higher order nets independent of classical digital nets. As alternative construction we consider polynomial lattice rules which first have been introduced by Niederreiter [3, 4] as special cases of digital nets and later generalized in [8] as special cases of digital higher order nets. Quasi-Monte Carlo rules using such point sets as integration nodes are nowadays known as higher order polynomial lattice rules. We refer to [9, 4] for more information on polynomial lattice rules. In [8] the existence of higher order polynomial lattice rules achieving optimal convergence rates was established with an averaging argument. In addition, these rules can at the same time achieve strong polynomial tractability. However,
3 3 this approach is purely probabilistic and gives no hint of how to construct such point sets. A further nonconstructive existence result for good higher order polynomial lattice rules is presented in [6]. This approach uses the concept of a figure of merit. The results presented in these papers were very encouraging. In particular, for large values of α and s, the higher order polynomial lattice rules improved on the construction based on classical digital nets from [3, Section 4.4]. In this paper we find a construction procedure producing higher order polynomial lattice rules which perform well when applied to numerical integration. In particular, we consider two cases: Firstly, we use a componentby-component CBC approach an idea first used in [6] to produce higher order polynomial lattice rules achieving the optimal rate of convergence for functions having higher order mixed partial derivatives, see Algorithm and Theorem 3.. Furthermore, by combining the CBC approach with a sieve - type algorithm as used in [0] we can even construct higher order polynomial lattice rules which automatically adjust themselves to the smoothness of the integrand in terms of the convergence of the integration error within a certain arbitrarily high range; see Algorithm 2 and Theorem 4.2. We point out already here, that an analogous result for lattice rules is not known. In addition, we study a special case of higher order polynomial lattice rules, namely higher order Korobov polynomial lattice rules as introduced in [8]; see also [] for the classical definition. We find that analogous results as for the higher order polynomial lattice rules constructed using a CBC and CBC sieve approach can be established. The approaches to constructing digital higher order nets presented in this paper are direct, that is, we avoid employing classical digital nets, which earlier constructions see [3, Section 4.4] relied on. Furthermore, the higher order polynomial lattice rules are constructed to perform well when applied to numerical integration and without reference to the quality parameter t of the resulting digital higher order net. The structure of the paper is as follows: In Section 2 we recall higher order polynomial lattice rules, discuss the function space under consideration and present a result on numerical integration in this function space employing higher order polynomial lattice rules. In Section 3 we use a CBC approach to construct higher order polynomial lattice rules achieving optimal rates of convergence for functions of a given smoothness and in Section 4 we show how to construct higher order polynomial lattice rules achieving optimal convergence rates for a given range of smoothness parameters using a CBC
4 4 sieve algorithm. Finally, in Section 5, analogous results for higher order Korobov polynomial lattice rules are established. 2 Preliminaries In this section we introduce higher order polynomial lattice rules which can achieve arbitrarily high convergence rates, the function space under consideration, and a result on numerical integration in this function space when using higher order polynomial lattice rules. 2. Polynomial lattice rules for arbitrarily smooth functions For a prime b let Z b be the finite field with b elements and let Z b x be the field of formal Laurent series over Z b. Elements of Z b x are formal Laurent series, L = t l x l, l=w where w is an arbitrary integer and all t l Z b. Note that Z b x contains the field of rational functions over Z b as a subfield. Further let Z b [x] be the set of all polynomials over Z b. For an integer n let v n be the map from Z b x to the interval [0, defined by n v n t l x l = t l b l. l=w l=max,w The following definition of higher order polynomial lattice rules given in [8] is a slight generalization of the definition from [3], see also [4]. Definition 2. Let b be prime and let m n be integers. For a given dimension s, choose px Z b [x] with degpx = n and let q x,..., q s x Z b [x]. For 0 h < b m let h = h 0 + h b + + h m b m be the b-adic expansion of h. With each such h we associate the polynomial hx = m r=0 h r x r Z b [x].
5 5 Then S p,m,n q, where q = q x,..., q s x, is the point set consisting of the b m points hxq x hxqs x x h = v n,..., v n [0, s, px px for 0 h < b m. A quasi-monte Carlo rule using the point set S p,m,n q is called a polynomial lattice rule. Remark 2. Using similar arguments as for the classical case n = m, see [3, 4], it can be shown that the point set S p,m,n q is a digital net in the sense of [3] which can be seen as a generalisation of the classical definition of digital nets according to Niederreiter [2, 3, 4]. The generating matrices C,..., C s Z n m b of this digital net can be obtained in the following way: For j s consider the expansions q j x px = u j l x l Z b x, l=w j where w j Z. Then the elements c j l,r of the n m matrix C j over Z b are given by c j l,r = uj r+l Z b, for j s, l n, 0 r m. For the rest of the paper, we make use of the following notation: We write h for vectors over Zb and h for vectors over Z or R. Polynomials over Z b are denoted by hx and vectors of polynomials by hx. Furthermore, given an integer h with b-adic expansion h = r=0 h rb r, we denote the associated polynomial by hx, which is given by n hx = h r x r. r=0 For arbitrary kx = k x,..., k s x Z b [x] s and qx = q x,..., q s x Z b [x] s, we define the inner product kx qx = s k j xq j x Z b [x],
6 6 and we write qx 0 mod px if px divides qx in Z b [x]. We remark here that for our results only the degree of the polynomial px is important and not the specific choice of px itself we assume though that px is irreducible, but this assumption could be removed by a more complicated analysis. 2.2 Walsh functions and the function space W α,s,γ We now define the space of functions we are going to study. This function space is based on Walsh functions whose definition is recalled in the following. Let N 0 denote the set of nonnegative and N the set of positive integers. Each k N has a unique b-adic representation k = a i=0 κ ib i with digits κ i {0,..., b } for 0 i a, where κ a 0. For k = 0 we have a = 0 and κ 0 = 0. Similarly, each x [0, has a b-adic representation x = i= ξ ib i with digits ξ i {0,..., b } for i. This representation is unique in the sense that infinitely many of the ξ i must differ from b. We define the kth Walsh function in base b, wal k : [0, C by wal k x := exp2πiξ κ ξ a+ κ a /b. For dimension s 2 and vectors k = k,..., k s N s 0 and x = x,..., x s [0, s we define wal k : [0, s C by wal k x := wal kj x j. It follows from the definition above that Walsh functions are piecewise constant functions. For more information on Walsh functions, see, e.g., [, 8] or [9, Appendix A]. When studying integration errors resulting from the approximation of an integral based on a digital net or digital higher order net or a higher order polynomial lattice rule, it is convenient to consider the Walsh series of the integrand f. In particular, for f L 2 [0, ] s, the Walsh series of f is given by fx fkwal k x, k N s 0 where the Walsh coefficients fk are given by fk = fxwal k x dx. [0,] s
7 7 In general, the Walsh series given in equation need not converge to f, however, for the space of Walsh series W α,s,γ, which we define in the following, it does, see also [3]. For more details on the convergence of Walsh series, we refer to [3] or [9]. Throughout the paper we assume that b is a fixed prime, all polynomials are over Z b [x] and all Walsh functions are also considered in the same base b. The function space under consideration in this paper is the space W α,s,γ L 2 [0, ] s as introduced in [3]. Here γ = γ j is a sequence of positive, nonincreasing weights, which are introduced to model the importance of different variables for our approximation problem, see [7]. For s N let [s] := {,..., s} and for u [s] let γ u := j u γ j be the weight associated with the projection onto components whose index is contained in u. Given a positive integer k with base b expansion k = κ b a + κ 2 b a2 + +κ v b av, a v < < a, v, we define µ α k := a + +a minv,α. Furthermore we put µ α 0 := 0. For k N 0 and a weight γ > 0, we define a function { if k = 0, r α γ, k := γb µαk otherwise. If we consider a vector k N s 0 of the form k = k,..., k s, we set r α γ, k := r α γ j, k j. Definition 2.2 The space W α,s,γ L 2 [0, ] s consists of all Walsh series f = k N fkwal s k for which the norm 0 f Wα,s,γ := sup fk k N s r 0 α γ, k. 2 is finite. For α 2, the following property was shown in [3]: Let f : [0, ] s R be such that all mixed partial derivatives up to order α in each variable are square integrable, then f W α,s,γ. Furthermore, an inequality using a Sobolev type norm and the norm in equation 2 was shown in [3], see also [2, 4]. Consequently, the results we are going to establish in the following
8 8 for functions in W α,s,γ also apply automatically to smooth functions. The assumption α > is needed to ensure that the sum of the absolute values of the Walsh coefficients converges. For the case α =, which requires a different analysis, we refer to [7] or to [9]. 2.3 Numerical Integration in W α,s,γ We are interested in the worst-case error of multivariate integration in W α,s,γ using a quasi-monte Carlo rule Q b m,s, which is given by eq b m,s, W α,s,γ = The initial error is given by eq 0,s, W α,s,γ = sup I s f Q b m,sf. 3 f W α,s,γ f Wα,s,γ sup I s f = I s. f W α,s,γ f Wα,s,γ We denote the quasi-monte Carlo rule based on a polynomial lattice rule S p,m,n q by Q b m,sq and the associated worst-case integration error by e b m,αq, p. The next proposition gives information on this quantity. Proposition 2. Let b be a prime and α 2 an integer. Then the worstcase integration error for multivariate integration in W α,s,γ using the polynomial lattice rule S p,m,n q is given by where e b m,αq, p = k D pq r α γ, k, D p q := { k N s 0 \ {0} : kx qx ax mod px with degax < n m}. 4 Proof. Combine [3, equation 5.2] with the determination of the dual net D of a polynomial lattice from [8, Section 4].
9 3 Component-by-component construction of polynomial lattice rules We propose the following algorithm to construct a polynomial lattice rule that achieves higher order convergence. We remark that unlike the results presented in Section 4, we only deal with a fixed α in this section. For ease of notation, we proceed as follows: We use q = qx Z b [x], p = px Z b [x] and a = ax Z b [x]; also, if we consider the polynomial associated with an integer k, we use k = kx Z b [x]. We put G b,n := {q Z b [x] : degq < n}. We also make use of the following lemma, which appeared in a weaker and nonexplicit form as [8, Lemma 4.2]. Lemma 3. Let α 2 be an integer. Then for every /α < λ we have where rαγ, λ l γ λ C b,α,λ, l= α C b,α,λ := C b α b,α,λ + b λα b b λi, i= { α if λ =, C b,α,λ = b b α b λ α if λ <. b b λ b λ α Furthermore, the series l= rλ αγ, l diverges to as λ goes to /α from the right. Proof. Let l = λ b a + + λ v b av where v, 0 < a v < < a and λ i {,..., b }. We divide the sum over all l N into two parts, namely firstly where v α and secondly where v > α. For the first part we have α b v v= b λa + +a v 0<a v< <a 9
10 0 = α b v v= α b b λ v= =: Cb,α,λ. a b λa a =v a 2 =v v = a v b λa 2 b λav a v= { α if λ =, b b α b λ α if λ <, b b λ b λ α For the second part we have b α = = = = baα b λa + +a α 0<a α< <a b α b b α b b α b b α b λα b a b λa a =α a 2 =α a α= α i= α i= b aα b λaα b λi b λi. Hence, we have shown that a α =a α+ a α= a α b λa 2 a α= b λa α b aα b λaα b λα aα b aα b λaα b λa 2 a 2 =a 3 + b λa a =a 2 + As b α b λα b α λ b α γ b λα b i= γ λ C b,α,λ + α i= b λi rαγ, λ l b α b λα b l= α i= =: γ λ C b λi b,α,λ. whenever λ /α from the right we also obtain b λi the second assertion. Now we show that a component-by-component approach can be used to construct a polynomial lattice rule that achieves higher order convergence, where for d s, we set q d = q,..., q d. Note that we consider this vector instead of, q 2,..., q s, c.f. [5, Algorithm 4.3], as otherwise the
11 projection onto the first coordinate does not achieve a convergence rate of b αm, see also [8, Remark 2.3]. The component-by-component algorithm for a fixed α is summarised in Algorithm. Algorithm CBC algorithm for fixed α Require: b a prime, s, m N and weights γ = γ j j. : Choose an irreducible polynomial p Z b [x], with degp = n. 2: for d = to s do 3: find q d G b,n by minimising e b m,αq,..., q d, p as a function of q d. 4: end for 5: return q = q,..., q s. Theorem 3. Let b be prime, let s, n, m, α N, m n and let α 2. Let p Z b [x] be irreducible with degp = n. Suppose q,..., q s G s b,n is constructed using Algorithm. Then for all d =,..., s we have: e b m,αq,..., q d, p b minτm,n d + 3γ /τ j C b,α,/τ τ τ < α. Proof. We firstly show the result for d =. By Proposition 2., e b m,αq, p = r α γ, k. k D pq The algorithm chooses q as to minimise the worst-case error, so we have Hence for all /α < λ we have e b m,αq, p e b m,αq, p, q G b,n. e b m,αq, p λ e b n b m,αq, p λ. q G b,n Using an argument very similar to the one used in the proof of [8, Proposition 4.3], it can be shown that for all /α < λ e b m,αq, p λ e b n b m,αq, p λ γ λ C b,α,λ b m + b λn. q G b,n
12 2 Consequently, setting τ = /λ we obtain e b m,αq, p + 2γ λ C b,α,λ /λ b minm/λ,n + 3γ /τ C b,α,/τ τ b minmτ,n. We now assume that for some d < s we have q d Gd b,n such that d e b m,αq d, p b minτm,n + 3γ /τ j C b,α,/τ τ. We consider e b m,αq d, q d+, p = = k,k d+ D pq d,q d+ k D pq d r α γ, k + r α γ, kr α γ d+, k d+ k d+ = = e b m,αq d, p + θq d, q d+, where we set θq d, q d+ := k d+ = r α γ d+, k d+ r α γ d+, k d+ k N d 0 k,k d+ D pq d,q d+ k N d 0 k,k d+ D pq d,q d+ r α γ, k. r α γ, k We see from Algorithm that qd+ is chosen in such a way that the worst-case error e b m,αq d, q d+, p is minimised. Since the only dependence on q d+ is in θq d, q d+ we have θq d, q d+ θq d, q d+ for all q d+ G b,n. This implies that for all /α < λ we have θq d, qd+ λ θq b d, q n d+ λ q d+ G b,n = b n q d+ G b,n k d+ = r α γ d+, k d+ k N d 0 k,k d+ D pq d,q d+ r α γ, k λ
13 3 b n q d+ G b,n k d+ = k d+ = p k d+ + b n rαγ λ d+, k d+ k d+ = p k d+ rαγ λ d+, k d+ r λ αγ d+, k d+ k N d 0 k q d a mod p dega<n m q d+ G b,n k N d 0 k,k d+ D pq d,q d+ r α γ, k k N d 0 k q d +k d+q d+ a mod p dega<n m λ r α γ, k λ r λ αγ, k, where we used Jensen s inequality, which states that for a sequence a k of nonnegative reals we have a k λ a λ k for any 0 < λ. Now we have k d+ = p k d+ rαγ λ d+, k d+ γλ d+ C b,α,λ b λn where we used the following: rαγ, λ k = k= p k + k N d 0 k q d a mod p dega<n m k N d 0 \{0} k q d a mod p dega<n m γλ d+ C b,α,λ b λn + eb m,αq d, p λ, rαγ, λ b n l + l= l=0 b n k= p k r α γ, k λ r α γ, k λ r λ αγ, k + b n l.
14 4 For l > 0 we have r α γ, b n l b n r α γ, l. Further for k < b n the polynomial p never divides k since degp = n. Hence rαγ, λ k = k= p k rαγ, λ b n l b λn l= l= r λ αγ, l γλ C b,α,λ b λn. Next we consider the case where k d+ is not a multiple of p. Here we have b n q d+ G b,n b n Now we have Hence k d+ = p k d+ k d+ = p k d+ q d+ G b,n = k N d 0 rαγ λ d+, k d+ r λ αγ d+, k d+ q d+ G b,n k N d 0 \{0} k q d +k d+q d+ a mod p dega<n m k N d 0 r λ αγ, k r λ αγ, kb n m d = b n m + C b,α,λ γj λ. k,k d+ D pq d,q d+ r α γ, k k N d 0 k q d +k d+q d+ a mod p dega<n m r λ αγ, k a Z b [x] q d+ G b,n dega<n m k q d +k d+ q d+ a mod p λ r λ αγ, k. b n k d+ = p k d+ r λ αγ d+, k d+ q d+ G b,n k N d 0 \{0} k q d +k d+q d+ a mod p dega<n m r λ αγ, k
15 5 Consequently, b n k d+ = b m C b,α,λγ λ d+ r λ αγ d+, k d+ b n m d + C b,α,λ γj λ d + C b,α,λ γj λ. θq d, q d+ γ λ d+ C b,α,λ b λn + e b m,αq d, p λ + b m C b,α,λγ λ d+ γ d+ C /λ b,α,λ [ d + C b,α,λ γj λ /λ b λn + e b m,αq d, p λ + b m /λ d + C b,α,λ γj ] λ. We now set τ = /λ and use the induction hypothesis to obtain θq d, qd+ γ d+ Cb,α,/τ τ b + e n/τ b m,αq d, p /τ + d + C b m b,α,/τ γ /τ τ γ d+ Cb,α,/τ τ 3 d + 3γ /τ b minm,n/τ j C b,α,/τ = Finally, we have 3 τ b minτm,n γ d+c τ b,α,/τ d + 3γ /τ j C b,α,/τ τ. e b m,αq d+, p = e b m,αq d, p + θq d, qd+ d + 3γ /τ b minτm,n j C b,α,/τ τ = 3 τ + b γ d+c τ minτm,n b,α,/τ d + 3γ /τ C b,α,/τ τ b minτm,n + 3τ γ d+ C τ b,α,/τ j d + 3γ /τ C b,α,/τ τ j j τ
16 6 d+ + 3γ /τ b minτm,n j C b,α,/τ τ, where we again used Jensen s inequality. From Theorem 3. we obtain the following corollary. Corollary 3. Let b be prime, let s, n, m, α N, m n and α 2. Let p Z b [x] be irreducible with degp = n and suppose q G s b,n is constructed using Algorithm. We have where e b m,αq, p γ α δ j c s,α,γ,δ := c s,α,γ,δ b minα δm,n 0 < δ α, α δ + 3γ α δ j C b,α,. α δ Suppose <, then c s,α,γ,δ c,α,γ,δ < and we have e b m,αq, p c,α,γ,δ 0 < δ α. b minα δm,n Thus the worst-case error is bounded independently of the dimension. Under the assumption A := lim sup s s γ j/log s < we obtain c s,α,γ,α c η s 2C b,α,a+η and therefore e b m,αq, p c ηs 2C b,α,a+η η > 0, b m where c η depends only on η. Thus the worst-case error satisfies a bound which depends only polynomially on the dimension. Proof. The first part follows from Theorem 3. by setting τ = α δ. The second and the third part follow from the first part in exactly the same way as in the proof of [5, Corollary 4.5]. The above result shows that higher order polynomial lattice rules can achieve a worst-case error satisfying at the same time the almost optimal convergence rate and a bound which depends only polynomially or even does not depend on the dimension s the technical term for such a behavior is strong polynomial tractability. Until now it is not known whether this is possible for ordinary lattice rules.
17 4 Optimal convergence rates for a range of smoothness parameters In this section, we construct polynomial lattices which are optimal for a range of smoothness parameters; we use α and τ α to denote the smoothness, where 2 α β, τ α < α. We set A m,n,s,α,p λ := e λ b sn b,αq m s, p. q s G s b,n Proposition 4. For α 2 and /α < λ we have 2 A m,n,s,α,p λ + + γ λ b minm,λn j C b,α,λ. Proof. Using Proposition 2. and Jensen s inequality we obtain, A m,n,s,α,p λ b sn = q G s b,n k D pq rαγ, λ k b sn k N s 0 \{0} r λ αγ, k q G s b,n k q a mod p dega<n m 7. 5 In the case where all components of k are multiples of p every q satisfies the equation k q 0 mod p and hence we have b sn q G s b,n k q a mod p dega<n m = and the sum over all k which satisfy this condition is therefore bounded by k N s 0 \{0} k 0 mod p r λ αγ, k = + rαγ λ j, k. k=0 p k
18 8 Now we have rαγ λ j, k = k=0 p k rαγ λ j, b n l + l=0 l=0 b n k= p k r λ αγ j, k + b n l. For l > 0 we have r α γ j, b n l b n r α γ j, l and further for k < b n the polynomial p never divides k since degp = n. Hence rαγ λ j, k = + k=0 p k l= r λ αγ j, b n l + b λn rαγ λ j, l. l= Therefore, k N s 0 \{0} k 0 mod p r λ αγ, k + = =u [s] b λn + b λn γj λ C b,α,λ b u λn γ λ u C u b,α,λ + + γj λ C b,α,λ. In the case where there is at least one component of k which is not a multiple of p we have = b sn b m q G s b,n k q a mod p dega<n m and therefore this part of equation 5 is bounded by b m k N s 0 \{0} k 0 mod p r λ αγ, k b m k N s 0 \{0} r λ αγ, k + b m + γj λ C b,α,λ.
19 9 Altogether we now obtain that A m,n,s,α,p λ b γ λ m b λn j C b,α,λ γ λ b minm,λn j C b,α,λ as required. G s b,βm where Let α β and set n = βm. Let ν denote the equiprobable measure on. For c and τ < α β the following set is introduced: C b,α c, τ := { q G s b,βm : e b m,αq, p E b,α,γ,s,m c, τ }, 6 Furthermore, let C b,α c := E b,α,γ,s,m c, τ := 2τ c τ τ<α b τm C b,α c, τ + τ + γ /τ j C b,α,/τ. = { q G s b,βm : e b m,αq, p E b,α,γ,s,m c, τ τ < α }. 7 Note that the intersection τ<α C b,αc, τ can be understood as an intersection of finitely many sets since C b,α c, τ has only finitely many elements. Lemma 4. Let c and τ < α β, then we have νc b,α c, τ > c. Proof. We denote C b,α c, τ := G s b,βm \ C b,αc, τ. Then for all τ < α we have A m,βm,s,α,p /τ = e /τ b sβm b m,αq, p q G s b,βm > νc b,α c, τ 2c b m + + γ /τ j C b,α,/τ. Now using Proposition 4. we obtain νc b,α c, τ < c and the result follows.
20 20 Lemma 4.2 Let c, then we have Proof. Let τ < α be such that νc b,α c > c. E b,α,γ,s,m c, τ = inf E b,α,γ,s,mc, τ τ<α note that by Lemma 3. we have E b,α,γ,s,m c, τ whenever τ α and hence we can find τ with the demanded property. Then we have C b,α c, τ C b,α c, τ = C b,α c τ<α and hence the result follows from Lemma 4.. If we choose c = β in Lemma 4.2, then we obtain νc b,α β > β and consequently we have β β β ν C b,α β = ν C b,α β νc b,α β > 0. α=2 α=2 Hence we obtain the following theorem which establishes the existence of a q G s b,βm which achieves the optimal convergence rate for a range of α s. Theorem 4. Let β, m, s N, β 2 and let p Z b [x] with degp = βm. Then there exists a q G s b,βm such that τα e b m,αq, p 2τα β τα + + γ /τα j C b,α,/τα 8 b ταm for all 2 α β and for all τ α < α. The proof of Theorem 4. suggests that in principle we can find q which satisfies equation 8 for all 2 α β and all τ α < α by using a so-called sieve algorithm which will be explained in the following. Use a computer search to find β b βms + of the b βms vectors q in G s b,βm which satisfy α=2 e b m,2q, p E b,2,γ,s,m β, τ 2 τ 2 < 2,
21 2 and label this set T 2. By Lemma 4.2 we know that at least such a number of vectors exists. Then proceed by using a computer search to find 2β b βms + vectors q in T 2 which satisfy e b m,3q, p E b,3,γ,s,m β, τ 3 τ 3 < 3 and label this set T 3. Since 3 3 ν C b,α β = ν C b,α β α=2 α=2 3 νc b,α β > 2 β, we know that there are at least 2β b βms + values in T 2 to populate the set T 3. In the same way we proceed to construct the sets T 4,..., T β. Theorem 4. guarantees that T β is not empty and we may select q to be any vector from T β. This vector satisfies equation 8 for all 2 α β and all τ α < α. However, in practice such a search algorithm would not be applicable since it is much too time consuming. For this reason we show in the following how the sieve algorithm may be combined with the component-by-component CBC algorithm which makes the search algorithm applicable. Such an algorithm, we call it CBC sieve algorithm, is presented in Algorithm 2. For its statement we use the following notation: For 2 α β and p Z b [x] with degp = βm we define the following: for d = 0 and q G b,βm we set θ α 0, q := e b m,αq, p, and for d N, q d G d b,βm and q d+ G b,βm we set α=2 θ α q d, q d+ := e b m,αq d, q d+, p e b m,αq d, p. Furthermore, for short we use the notation M d,α,γ τ := b m d + 3βγ /τ j C b,α,/τ. Now we prove the following result.
22 22 Algorithm 2 CBC sieve algorithm for 2 α β Require: b a prime, s, m, β N, β 2, and p Z b [x] with degp = βm. : Set T,d := G b,βm for all d s and q 0 = 0. 2: for d = 0 to s do 3: for α = 2 to β do 4: perform a computer search to find α β b βm + elements q in T α,d+ to populate the set T α,d+, which is a subset of 5: if d = 0 then 6: { q T α,d+ : θ α 0, q } τα + 3γ /τα b ταm C b,α,/τα τα < α 7: else 8: { q T α,d+ : θ α q d, q 3βγ /τα d+ C τα } b,α,/τ α M d,α,γ τ α τα < α 9: end if 0: end for : Select q d+ T β,d+. 2: Set q d+ = q d, q d+. 3: end for 4: return q = q s. Theorem 4.2 Let s, m, β N, β 2, then Algorithm 2 constructs a vector q d Gd b,βm such that e b m,αq d, p b ταm d + 3βγ /τα C b,α,/τα τα for all τ α < α and for all 2 α β. To prove Theorem 4.2 we introduce the following set: for q d G d b,βm let F b,α c, q d be the set of all q d+ G b,βm such that θ α q d, q d+ j 3cγ /τα d+ C b,α,/τ α M d,α,γ τ α τα 9
23 23 for all τ α < α. Lemma 4.3 Let 2 α β and let c. q d G d b,βm such that Assume that there exists a for all τ α < α and for all 2 α β. Then e b m,αq d, p M d,α,γ τ α τα 0 ν F b,α c, q d > c. Proof. From the proof of Theorem 3. and using assumption 0 for all /α < λ we have b βm q d+ G b,βm θ α q d, q d+ λ γd+c λ b,α,λ b + e λαm b m,αq d, p λ + 3γ λ d+c b,α,λ M d,α,γ /λ. b m d + γj λ C b,α,λ From this the result follows in the same way as in the proof of Lemmas 4. and 4.2. Now we give the proof of Theorem 4.2. Proof. We proceed by induction on d and firstly show the result for d =. Having fixed d =, we proceed by induction on α. For q C b,α β we have see equation 7 e b m,αq, p + 3βγ/τα bταm C b,α,/τα τα τ α < α, for 2 α β. According to Lemma 4.2, νc b,α β > β, 2 α β, hence there are β b βm + elements to populate T 2,. Assume now that for 2 α < β, there are α β b βm + elements to populate T α,, hence νt α, > α β. We want to show that ν {q T α, : e b m,α+q, p M,α+,γ τ α+ τ α+ τ α+ < α + } > αβ,
24 24 which implies that there are αβ b βm + elements to populate T α+,. But {q T α, : e b m,α+q, p M,α+,γ τ α+ τ α+ τ α+ < α + } = T α, {q G b,βm : e b m,α+q, p M,α+,γ τ α+ τ α+ τ α+ < α + }, hence we get equation from the induction assumption and from Lemma 4.2. Thus we have proven the assertion for d =. We now assume that for d < s, the algorithm has found q d so that e b m,αq d, p M d,α,γ τ α τα 2 for all τ α < α and for all 2 α β and again we proceed by induction. According to Lemma 4.3, under the assumption 2, we have ν F b,α β, q d > β 2 α β, hence there are β b βm + elements to populate T 2,d+. We now assume that for 2 α < β, there are α β b βm + elements to populate T α,d+, hence νt α,d+ > α β. Since { q T α,d+ : θq d, q 3βγ /τ τα+ α+ d+ C b,α+,/τα+ M d+,α+,γ τ α+ } τ α+ < α + = T α,d+ F b,α+ β, q d we obtain from the inductive hypothesis and from Lemma 4.3 that { ν q T α,d+ : θq d, q 3βγ /τ α+ d+ C b,α+,/τα+ M d+,α+,γ τ α+ } τ α+ < α + > αβ, τα+ which implies that there are αβ b βm + elements to populate T α+,d+. Therefore Algorithm 2 finds a q d+ G b,βm such that θ α q d, q d+ 3βγ /τα d+ C τα b,α,/τ α M d,α,γ τ α
25 25 for all τ α < α and for all 2 α β. Using equation 2 we obtain, e b m,αq d, q d+, p = e b m,αq d, p + θ α q d, q d+ M d,α,γ τ α τα + 3βγ /τα d+ C b,α,/τ α τα M d+,α,γ τ α τα, for all τ α < α and for all 2 α β. 5 Optimal convergence rates for a range of smoothness parameters using Korobov polynomial lattice rules In this section we study a special case of polynomial lattice rules, namely Korobov polynomial lattice rules. We show the existence of higher order Korobov polynomial lattice rules which achieve optimal rates of convergence for a range of smoothness parameters and present an algorithm which shows how to construct such higher order Korobov polynomial lattice rules. This algorithm is the same as the sieve algorithm presented in Section 4, but due to the structure of Korobov polynomial lattice rules, the cost of such an algorithm is feasible. We now present the results which are used to establish the existence of a higher order Korobov polynomial lattice rule achieving optimal rates of convergence for a range of smoothness parameters and its construction. For the remainder of this section, we use φq := q, q 2,..., q s mod p, q G b,βm to denote the generating vector of the higher order Korobov polynomial lattice rule S p,m,βm φq and e b m,αφq, p to denote the corresponding worst-case error, 2 α β; we recall that α β and n = βm. As in Section 3 we point out that we use generating vectors φq := q, q 2,..., q s mod p instead of, q,..., q s see e.g. [5, Algorithm 4.6], as otherwise the projection onto the first coordinate does not achieve a convergence rate of b αm. We start with the following proposition, which is analogous to Proposition 4., where we set à m,n,s,α,p λ := e λ b n b,αφq, p. m q G b,n
26 26 Proposition 5. For α 2 and /α < λ we have [ ] Ã m,n,s,α,p λ s γ λ b minm,λn j C b,α,λ. Proof. Using Proposition 2. and Jensen s inequality we obtain, Ã m,n,s,α,p λ r b αγ, λ k n = q G b,βm k D pφq k N s 0 \{0} r λ αγ, k b n q G b,βm k φq a mod p dega<n m. 3 In the case where all components of k are multiples of p every q G b,βm satisfies the equation k φq 0 mod p and hence we have = b n k N s 0 \{0} k 0 mod p q G b,βm k φq a mod p dega<n m and the sum over all k which satisfy this condition is therefore bounded by rαγ, λ k + + γ λ b λn j C b,α,λ, see the proof of Proposition 4.. In the case where there is at least one component of k which is not a multiple of p we have b n q G b,βm k φq a mod p dega<n m sb m, because for any choice a there are at most s solutions q to k φq a mod p and there are b n m possible choices for a. Therefore this part of equation 3 is bounded by s b m + + γj λ C b,α,λ.
27 27 Altogether we now obtain that à m,n,s,α,p λ s b m + b λn s + b minm,λn γj λ C b,α,λ + γj λ C b,α,λ as required. Let ν denote the equiprobable measure on G b,βm. For c and τ < α β the following set is introduced: { } C b,α c, τ := q G b,βm : e b m,αφq, p Ẽb,α,γ,s,mc, τ, 4 where Ẽ b,α,γ,s,m c, τ := cτ s + τ b τm τ + + γ /τ j C b,α,/τ. Furthermore, let C b,α c := = τ<α C b,α c, τ 5 { } q G b,βm : e b m,αφq, p Ẽb,α,γ,s,mc, τ τ < α. Lemma 5. Let c and τ < α β. Then we have ν C b,α c, τ > c. Proof. The proof follows exactly along the lines of the proof of Lemma 4.. Lemma 5.2 Let c. Then we have ν C b,α c > c. Proof. The proof follows exactly along the lines of the proof of Lemma 4.2.
28 28 Algorithm 3 Korobov sieve algorithm Require: b a prime, s, m, β N, β 2, and p Z b [x] with degp = βm. : Set T := G b,βm. 2: for α = 2 to β do 3: perform a computer search to find α β b βm + elements q in T α to populate the set T α, which is a subset of { } q T α : e b m,αφq, p Ẽb,α,γ,s,mβ, τ α τ α < α 4: end for 5: Select q T β 6: return q. As in Section 4 we now introduce a sieve algorithm see Algorithm 3 which shows how to obtain a generating vector for a higher order Korobov polynomial lattice rule, which achieves optimal convergence rates for a range of smoothness parameters. The next theorem shows that Algorithm 3 does indeed produce such a vector. Theorem 5. Let s, m, β N, β 2. Then Algorithm 3 finds an element q G b,βm such that e b m,αφq, p s + τα τα β τα + + γ /τα j C b,α,/τα, b ταm for all τ α < α, 2 α β. Proof. We prove the result by induction on α. For α = 2, by Lemma 5.2, ν Cb,α β > β, 2 α β, so there are at least β b βm + vectors to populate the set T 2. We now assume that there are α β b βm + elements in the set T α, where 2 α < β, hence νt α > α β. We want to show that { } ν q T α : e b m,α+φq, p Ẽb,α+,γ,s,mβ, τ α+ τ α+ < α + > αβ, 6
29 29 which implies that there are αβ b βm + elements to populate the set T α+. Since { } q T α : e b m,α+φq, p Ẽb,α+,γ,s,mβ, τ α+ τ α+ < α + = T α C b,α+ β we obtain equation 6 from the induction assumption and from Lemma 5.2. References [] H.E. Chrestenson, A class of generalized Walsh functions, Pacific J. Math [2] J. Dick, Explicit constructions of quasi-monte Carlo rules for the numerical integration of high-dimensional periodic functions, SIAM J. Numer. Anal [3] J. Dick, Walsh spaces containing smooth functions and quasi-monte Carlo rules of arbitrary high order, SIAM J. Numer. Anal [4] J. Dick, The decay of the Walsh coefficients of smooth functions, Bull. Austral. Math. Soc [5] J. Dick, F. Kuo, F. Pillichshammer, I. Sloan, Construction algorithms for polynomial lattice rules for multivariate integration, Math. Comp [6] J. Dick, P. Kritzer, F. Pillichshammer, W. Ch. Schmid, On the existence of higher order polynomial lattices based on a generalized figure of merit, J. Complexity [7] J. Dick, F. Pillichshammer, Multivariate integration in weighted Hilbert spaces based on Walsh functions and weighted Sobolev spaces, J. Complexity [8] J. Dick, F. Pillichshammer, Strong tractability of multivariate integration of arbitrary high order using digitally shifted polynomial lattice rules, J. Complexity
30 30 [9] J. Dick, F. Pillichshammer, Digital Nets and Sequences, Cambridge University Press, Cambridge, 200 to appear. [0] J. Dick, F. Pillichshammer, B. Waterhouse, The construction of good extensible rank- lattices, Math. Comp [] N.M. Korobov, Properties and calculation of optimal coefficients, Dokl. Akad. Nauk SSSR In Russian. [2] H. Niederreiter, Point sets and sequences with small discrepancy, Monatsh. Math [3] H. Niederreiter, Low-discrepancy point sets obtained by digital constructions over finite fields, Czech. Math. J [4] H. Niederreiter, Random Number Generation and Quasi-Monte Carlo Methods, CBMS NSF Series in Applied Mathematics 63, SIAM, Philadelphia, 992. [5] I. H. Sloan, S. Joe, Lattice methods for multiple integration, Oxford Science Publications, The Clarendon Press Oxford University Press, New York, 994. [6] I. Sloan, A. Reztsov, Component-by-component construction of good lattice rules, Math. Comp [7] I.H. Sloan, H. Woźniakowski, When are quasi-monte Carlo algorithms efficient for high dimensional integrals?, J. Complexity [8] J.L. Walsh, A closed set of normal orthogonal functions, Amer. J. Math Author s Addresses: Jan Baldeaux and Josef Dick, School of Mathematics and Statistics, The University of New South Wales, Sydney, NSW 2052, Australia. janbaldeaux@gmail.com, josef.dick@unsw.edu.au Julia Greslehner and Friedrich Pillichshammer, Institut für Finanzmathematik, Universität Linz, Altenbergstraße 69, A-4040 Linz, Austria. julia.greslehner@gmx.at, friedrich.pillichshammer@jku.at
On the existence of higher order polynomial. lattices with large figure of merit ϱ α.
On the existence of higher order polynomial lattices with large figure of merit Josef Dick a, Peter Kritzer b, Friedrich Pillichshammer c, Wolfgang Ch. Schmid b, a School of Mathematics, University of
More informationAPPLIED MATHEMATICS REPORT AMR04/2 DIAPHONY, DISCREPANCY, SPECTRAL TEST AND WORST-CASE ERROR. J. Dick and F. Pillichshammer
APPLIED MATHEMATICS REPORT AMR4/2 DIAPHONY, DISCREPANCY, SPECTRAL TEST AND WORST-CASE ERROR J. Dick and F. Pillichshammer January, 24 Diaphony, discrepancy, spectral test and worst-case error Josef Dick
More informationApproximation of Functions Using Digital Nets
Approximation of Functions Using Digital Nets Josef Dick 1, Peter Kritzer 2, and Frances Y. Kuo 3 1 UNSW Asia, 1 Kay Siang Road, Singapore 248922. E-mail: j.dick@unswasia.edu.sg 2 Fachbereich Mathematik,
More informationarxiv: v2 [math.na] 18 May 2011
Efficient calculation of the worst-case error and (fast component-by-component construction of higher order polynomial lattice rules arxiv:1105.2599v2 [math.na] 18 May 2011 Jan Baldeaux Josef Dick Gunther
More informationConstruction algorithms for plane nets in base b
Construction algorithms for plane nets in base b Gunther Leobacher, Friedrich Pillichshammer and Thomas Schell arxiv:1506.03201v1 [math.nt] 10 Jun 2015 Abstract The class of (0, m, s)-nets in base b has
More informationExponential Convergence and Tractability of Multivariate Integration for Korobov Spaces
Exponential Convergence and Tractability of Multivariate Integration for Korobov Spaces Josef Dick, Gerhard Larcher, Friedrich Pillichshammer and Henryk Woźniakowski March 12, 2010 Abstract In this paper
More informationDyadic diaphony of digital sequences
Dyadic diaphony of digital sequences Friedrich Pillichshammer Abstract The dyadic diaphony is a quantitative measure for the irregularity of distribution of a sequence in the unit cube In this paper we
More informationDigital lattice rules
Digital lattice rules Multivariate integration and discrepancy estimates A thesis presented to The University of New South Wales in fulfilment of the thesis requirement for the degree of Doctor of Philosophy
More informationA NEW UPPER BOUND ON THE STAR DISCREPANCY OF (0,1)-SEQUENCES. Peter Kritzer 1. Abstract
INTEGERS: ELECTRONIC JOURNAL OF COMBINATORIAL NUMBER THEORY 5(3 (2005, #A11 A NEW UPPER BOUND ON THE STAR DISCREPANCY OF (0,1-SEQUENCES Peter Kritzer 1 Department of Mathematics, University of Salzburg,
More informationTutorial on quasi-monte Carlo methods
Tutorial on quasi-monte Carlo methods Josef Dick School of Mathematics and Statistics, UNSW, Sydney, Australia josef.dick@unsw.edu.au Comparison: MCMC, MC, QMC Roughly speaking: Markov chain Monte Carlo
More informationAPPLIED MATHEMATICS REPORT AMR04/16 FINITE-ORDER WEIGHTS IMPLY TRACTABILITY OF MULTIVARIATE INTEGRATION. I.H. Sloan, X. Wang and H.
APPLIED MATHEMATICS REPORT AMR04/16 FINITE-ORDER WEIGHTS IMPLY TRACTABILITY OF MULTIVARIATE INTEGRATION I.H. Sloan, X. Wang and H. Wozniakowski Published in Journal of Complexity, Volume 20, Number 1,
More informationNew Multilevel Algorithms Based on Quasi-Monte Carlo Point Sets
New Multilevel Based on Quasi-Monte Carlo Point Sets Michael Gnewuch Institut für Informatik Christian-Albrechts-Universität Kiel 1 February 13, 2012 Based on Joint Work with Jan Baldeaux (UTS Sydney)
More informationApplication of QMC methods to PDEs with random coefficients
Application of QMC methods to PDEs with random coefficients a survey of analysis and implementation Frances Kuo f.kuo@unsw.edu.au University of New South Wales, Sydney, Australia joint work with Ivan Graham
More informationarxiv: v1 [math.na] 13 Jul 2017
Richardson extrapolation of polynomial lattice rules Josef Dick, Takashi Goda, Takehito Yoshiki arxiv:707.03989v [math.na] 3 Jul 207 July 4, 207 Abstract We study multivariate numerical integration of
More informationLattice rules and sequences for non periodic smooth integrands
Lattice rules and sequences for non periodic smooth integrands dirk.nuyens@cs.kuleuven.be Department of Computer Science KU Leuven, Belgium Joint work with Josef Dick (UNSW) and Friedrich Pillichshammer
More informationInfinite-Dimensional Integration on Weighted Hilbert Spaces
Infinite-Dimensional Integration on Weighted Hilbert Spaces Michael Gnewuch Department of Computer Science, Columbia University, 1214 Amsterdam Avenue, New Yor, NY 10027, USA Abstract We study the numerical
More informationLow-discrepancy sequences obtained from algebraic function fields over finite fields
ACTA ARITHMETICA LXXII.3 (1995) Low-discrepancy sequences obtained from algebraic function fields over finite fields by Harald Niederreiter (Wien) and Chaoping Xing (Hefei) 1. Introduction. We present
More informationL 2 Discrepancy of Two-Dimensional Digitally Shifted Hammersley Point Sets in Base b
L Discrepancy of Two-Dimensional Digitally Shifted Hammersley Point Sets in Base Henri Faure and Friedrich Pillichshammer Astract We give an exact formula for the L discrepancy of two-dimensional digitally
More informationOptimal Randomized Algorithms for Integration on Function Spaces with underlying ANOVA decomposition
Optimal Randomized on Function Spaces with underlying ANOVA decomposition Michael Gnewuch 1 University of Kaiserslautern, Germany October 16, 2013 Based on Joint Work with Jan Baldeaux (UTS Sydney) & Josef
More informationOn the discrepancy estimate of normal numbers
ACTA ARITHMETICA LXXXVIII.2 (1999 On the discrepancy estimate of normal numbers 1. Introduction by M. B. Levin (Tel-Aviv Dedicated to Professor N. M. Korobov on the occasion of his 80th birthday 1.1. A
More informationSolutions to Tutorial 8 (Week 9)
The University of Sydney School of Mathematics and Statistics Solutions to Tutorial 8 (Week 9) MATH3961: Metric Spaces (Advanced) Semester 1, 2018 Web Page: http://www.maths.usyd.edu.au/u/ug/sm/math3961/
More informationIncomplete exponential sums over finite fields and their applications to new inversive pseudorandom number generators
ACTA ARITHMETICA XCIII.4 (2000 Incomplete exponential sums over finite fields and their applications to new inversive pseudorandom number generators by Harald Niederreiter and Arne Winterhof (Wien 1. Introduction.
More informationConstructions of digital nets using global function fields
ACTA ARITHMETICA 105.3 (2002) Constructions of digital nets using global function fields by Harald Niederreiter (Singapore) and Ferruh Özbudak (Ankara) 1. Introduction. The theory of (t, m, s)-nets and
More informationPart III. Quasi Monte Carlo methods 146/349
Part III Quasi Monte Carlo methods 46/349 Outline Quasi Monte Carlo methods 47/349 Quasi Monte Carlo methods Let y = (y,...,y N ) be a vector of independent uniform random variables in Γ = [0,] N, u(y)
More informationLow Discrepancy Sequences in High Dimensions: How Well Are Their Projections Distributed?
Low Discrepancy Sequences in High Dimensions: How Well Are Their Projections Distributed? Xiaoqun Wang 1,2 and Ian H. Sloan 2 1 Department of Mathematical Sciences, Tsinghua University, Beijing 100084,
More informationarxiv: v2 [math.nt] 27 Mar 2018
arxiv:1710.09313v [math.nt] 7 Mar 018 The Champernowne constant is not Poissonian Ísabel Pirsic, Wolfgang Stockinger Abstract We say that a sequence (x n n N in [0,1 has Poissonian pair correlations if
More informationConstruction algorithms for good extensible lattice rules
Construction algorithms for good extensible lattice rules Harald Niederreiter and Friedrich Pillichshammer Dedicated to Ian H. Sloan on the occasion of his 70th birthday Abstract Extensible olynomial lattice
More informationIan Sloan and Lattice Rules
www.oeaw.ac.at Ian Sloan and Lattice Rules P. Kritzer, H. Niederreiter, F. Pillichshammer RICAM-Report 2016-34 www.ricam.oeaw.ac.at Ian Sloan and Lattice Rules Peter Kritzer, Harald Niederreiter, Friedrich
More informationf(x n ) [0,1[ s f(x) dx.
ACTA ARITHMETICA LXXX.2 (1997) Dyadic diaphony by Peter Hellekalek and Hannes Leeb (Salzburg) 1. Introduction. Diaphony (see Zinterhof [13] and Kuipers and Niederreiter [6, Exercise 5.27, p. 162]) is a
More informationAppendix B for The Evolution of Strategic Sophistication (Intended for Online Publication)
Appendix B for The Evolution of Strategic Sophistication (Intended for Online Publication) Nikolaus Robalino and Arthur Robson Appendix B: Proof of Theorem 2 This appendix contains the proof of Theorem
More informationJoshua N. Cooper Department of Mathematics, Courant Institute of Mathematics, New York University, New York, NY 10012, USA.
CONTINUED FRACTIONS WITH PARTIAL QUOTIENTS BOUNDED IN AVERAGE Joshua N. Cooper Department of Mathematics, Courant Institute of Mathematics, New York University, New York, NY 10012, USA cooper@cims.nyu.edu
More informationIN QUASI-MONTE CARLO INTEGRATION
MATHEMATICS OF COMPUTATION VOLUME, NUMBER 0 JULY, PAGES - GOOD PARAMETERS FOR A CLASS OF NODE SETS IN QUASI-MONTE CARLO INTEGRATION TOM HANSEN, GARY L. MULLEN, AND HARALD NIEDERREITER Dedicated to the
More informationDiscrepancy Bounds for Nets and Sequences
Discrepancy Bounds for Nets and Sequences Peter Kritzer Johann Radon Institute for Computational and Applied Mathematics (RICAM) Austrian Academy of Sciences Linz, Austria Joint work with Henri Faure (Marseille)
More informationOn the Behavior of the Weighted Star Discrepancy Bounds for Shifted Lattice Rules
On the Behavior of the Weighted Star Discrepancy Bounds for Shifted Lattice Rules Vasile Sinescu and Pierre L Ecuyer Abstract We examine the question of constructing shifted lattice rules of rank one with
More informationSTRONG TRACTABILITY OF MULTIVARIATE INTEGRATION USING QUASI MONTE CARLO ALGORITHMS
MATHEMATICS OF COMPUTATION Volume 72, Number 242, Pages 823 838 S 0025-5718(02)01440-0 Article electronically published on June 25, 2002 STRONG TRACTABILITY OF MULTIVARIATE INTEGRATION USING QUASI MONTE
More information2014:05 Incremental Greedy Algorithm and its Applications in Numerical Integration. V. Temlyakov
INTERDISCIPLINARY MATHEMATICS INSTITUTE 2014:05 Incremental Greedy Algorithm and its Applications in Numerical Integration V. Temlyakov IMI PREPRINT SERIES COLLEGE OF ARTS AND SCIENCES UNIVERSITY OF SOUTH
More informationGeometrical Constructions for Ordered Orthogonal Arrays and (T, M, S)-Nets
Geometrical Constructions for Ordered Orthogonal Arrays and (T, M, S)-Nets Ryoh Fuji-Hara and Ying Miao Institute of Policy and Planning Sciences University of Tsukuba Tsukuba 305-8573, Japan fujihara@sk.tsukuba.ac.jp
More informationApproximation exponents for algebraic functions in positive characteristic
ACTA ARITHMETICA LX.4 (1992) Approximation exponents for algebraic functions in positive characteristic by Bernard de Mathan (Talence) In this paper, we study rational approximations for algebraic functions
More informationMA257: INTRODUCTION TO NUMBER THEORY LECTURE NOTES
MA257: INTRODUCTION TO NUMBER THEORY LECTURE NOTES 2018 57 5. p-adic Numbers 5.1. Motivating examples. We all know that 2 is irrational, so that 2 is not a square in the rational field Q, but that we can
More informationON SUMS OF PRIMES FROM BEATTY SEQUENCES. Angel V. Kumchev 1 Department of Mathematics, Towson University, Towson, MD , U.S.A.
INTEGERS: ELECTRONIC JOURNAL OF COMBINATORIAL NUMBER THEORY 8 (2008), #A08 ON SUMS OF PRIMES FROM BEATTY SEQUENCES Angel V. Kumchev 1 Department of Mathematics, Towson University, Towson, MD 21252-0001,
More informationQuasi-Monte Carlo integration over the Euclidean space and applications
Quasi-Monte Carlo integration over the Euclidean space and applications f.kuo@unsw.edu.au University of New South Wales, Sydney, Australia joint work with James Nichols (UNSW) Journal of Complexity 30
More informationKatholieke Universiteit Leuven Department of Computer Science
Extensions of Fibonacci lattice rules Ronald Cools Dirk Nuyens Report TW 545, August 2009 Katholieke Universiteit Leuven Department of Computer Science Celestijnenlaan 200A B-3001 Heverlee (Belgium Extensions
More informationQuasi-Monte Carlo Methods for Applications in Statistics
Quasi-Monte Carlo Methods for Applications in Statistics Weights for QMC in Statistics Vasile Sinescu (UNSW) Weights for QMC in Statistics MCQMC February 2012 1 / 24 Quasi-Monte Carlo Methods for Applications
More informationANOVA decomposition of convex piecewise linear functions
ANOVA decomposition of convex piecewise linear functions W. Römisch Abstract Piecewise linear convex functions arise as integrands in stochastic programs. They are Lipschitz continuous on their domain,
More informationIntegration of permutation-invariant. functions
permutation-invariant Markus Weimar Philipps-University Marburg Joint work with Dirk Nuyens and Gowri Suryanarayana (KU Leuven, Belgium) MCQMC2014, Leuven April 06-11, 2014 un Outline Permutation-invariant
More informationConvergence rate analysis for averaged fixed point iterations in the presence of Hölder regularity
Convergence rate analysis for averaged fixed point iterations in the presence of Hölder regularity Jonathan M. Borwein Guoyin Li Matthew K. Tam October 23, 205 Abstract In this paper, we establish sublinear
More informationBasic Algebra. Final Version, August, 2006 For Publication by Birkhäuser Boston Along with a Companion Volume Advanced Algebra In the Series
Basic Algebra Final Version, August, 2006 For Publication by Birkhäuser Boston Along with a Companion Volume Advanced Algebra In the Series Cornerstones Selected Pages from Chapter I: pp. 1 15 Anthony
More informationA strongly polynomial algorithm for linear systems having a binary solution
A strongly polynomial algorithm for linear systems having a binary solution Sergei Chubanov Institute of Information Systems at the University of Siegen, Germany e-mail: sergei.chubanov@uni-siegen.de 7th
More informationMACMAHON S PARTITION ANALYSIS IX: k-gon PARTITIONS
MACMAHON S PARTITION ANALYSIS IX: -GON PARTITIONS GEORGE E. ANDREWS, PETER PAULE, AND AXEL RIESE Dedicated to George Szeeres on the occasion of his 90th birthday Abstract. MacMahon devoted a significant
More informationQMC designs and covering of spheres by spherical caps
QMC designs and covering of spheres by spherical caps Ian H. Sloan University of New South Wales, Sydney, Australia Joint work with Johann Brauchart, Josef Dick, Ed Saff (Vanderbilt), Rob Womersley and
More informationExplicit constructions of quasi-monte Carlo rules for the numerical integration of high dimensional periodic functions
Explicit constructions of quasi-monte Carlo rules for the numerical integration of high dimensional periodic functions By J. Dick at Sydney Astract In this paper we give first explicit constructions of
More informationNUMBER FIELDS WITHOUT SMALL GENERATORS
NUMBER FIELDS WITHOUT SMALL GENERATORS JEFFREY D. VAALER AND MARTIN WIDMER Abstract. Let D > be an integer, and let b = b(d) > be its smallest divisor. We show that there are infinitely many number fields
More informationAbsolute Values and Completions
Absolute Values and Completions B.Sury This article is in the nature of a survey of the theory of complete fields. It is not exhaustive but serves the purpose of familiarising the readers with the basic
More informationSeries of Error Terms for Rational Approximations of Irrational Numbers
2 3 47 6 23 Journal of Integer Sequences, Vol. 4 20, Article..4 Series of Error Terms for Rational Approximations of Irrational Numbers Carsten Elsner Fachhochschule für die Wirtschaft Hannover Freundallee
More informationQMC methods in quantitative finance. and perspectives
, tradition and perspectives Johannes Kepler University Linz (JKU) WU Research Seminar What is the content of the talk Valuation of financial derivatives in stochastic market models using (QMC-)simulation
More informationOn Gauss sums and the evaluation of Stechkin s constant
On Gauss sums and the evaluation of Stechkin s constant William D. Banks Department of Mathematics University of Missouri Columbia, MO 65211 USA bankswd@missouri.edu Igor E. Shparlinski Department of Pure
More informationNOTES ON DIOPHANTINE APPROXIMATION
NOTES ON DIOPHANTINE APPROXIMATION Jan-Hendrik Evertse January 29, 200 9 p-adic Numbers Literature: N. Koblitz, p-adic Numbers, p-adic Analysis, and Zeta-Functions, 2nd edition, Graduate Texts in Mathematics
More informationuniform distribution theory
Uniform Distribution Theory 2 (2007), no.1, 1 10 uniform distribution theory UNIFORM DISTRIBUTION OF SEQUENCES CONNECTED WITH THE WEIGHTED SUM-OF-DIGITS FUNCTION Friedrich Pillichshammer ABSTRACT. In this
More informationDIOPHANTINE APPROXIMATION AND CONTINUED FRACTIONS IN POWER SERIES FIELDS
DIOPHANTINE APPROXIMATION AND CONTINUED FRACTIONS IN POWER SERIES FIELDS A. LASJAUNIAS Avec admiration pour Klaus Roth 1. Introduction and Notation 1.1. The Fields of Power Series F(q) or F(K). Let p be
More informationProgress in high-dimensional numerical integration and its application to stochastic optimization
Progress in high-dimensional numerical integration and its application to stochastic optimization W. Römisch Humboldt-University Berlin Department of Mathematics www.math.hu-berlin.de/~romisch Short course,
More informationUniform distribution of sequences connected with the weighted sum-of-digits function
Uniform distribution of seuences connected with the weighted sum-of-digits function Friedrich Pillichshammer Abstract In this paper we consider seuences which are connected with the so-called weighted
More informationTo Professor W. M. Schmidt on his 60th birthday
ACTA ARITHMETICA LXVII.3 (1994) On the irreducibility of neighbouring polynomials by K. Győry (Debrecen) To Professor W. M. Schmidt on his 60th birthday 1. Introduction. Denote by P the length of a polynomial
More informationMATH FINAL EXAM REVIEW HINTS
MATH 109 - FINAL EXAM REVIEW HINTS Answer: Answer: 1. Cardinality (1) Let a < b be two real numbers and define f : (0, 1) (a, b) by f(t) = (1 t)a + tb. (a) Prove that f is a bijection. (b) Prove that any
More informationReal Analysis Math 131AH Rudin, Chapter #1. Dominique Abdi
Real Analysis Math 3AH Rudin, Chapter # Dominique Abdi.. If r is rational (r 0) and x is irrational, prove that r + x and rx are irrational. Solution. Assume the contrary, that r+x and rx are rational.
More informationA Belgian view on lattice rules
A Belgian view on lattice rules Ronald Cools Dept. of Computer Science, KU Leuven Linz, Austria, October 14 18, 2013 Atomium Brussels built in 1958 height 103m figure = 2e coin 5 10 6 in circulation Body
More informationHomework #2 Solutions Due: September 5, for all n N n 3 = n2 (n + 1) 2 4
Do the following exercises from the text: Chapter (Section 3):, 1, 17(a)-(b), 3 Prove that 1 3 + 3 + + n 3 n (n + 1) for all n N Proof The proof is by induction on n For n N, let S(n) be the statement
More informationFAST QMC MATRIX-VECTOR MULTIPLICATION
FAST QMC MATRIX-VECTOR MULTIPLICATION JOSEF DICK, FRANCES Y. KUO, QUOC T. LE GIA, CHRISTOPH SCHWAB Abstract. Quasi-Monte Carlo (QMC) rules 1/N N 1 n=0 f(y na) can be used to approximate integrals of the
More informationChapter 8. P-adic numbers. 8.1 Absolute values
Chapter 8 P-adic numbers Literature: N. Koblitz, p-adic Numbers, p-adic Analysis, and Zeta-Functions, 2nd edition, Graduate Texts in Mathematics 58, Springer Verlag 1984, corrected 2nd printing 1996, Chap.
More informationFactorization of integer-valued polynomials with square-free denominator
accepted by Comm. Algebra (2013) Factorization of integer-valued polynomials with square-free denominator Giulio Peruginelli September 9, 2013 Dedicated to Marco Fontana on the occasion of his 65th birthday
More informationAn intractability result for multiple integration. I.H.Sloan y and H.Wozniakowski yy. University of New South Wales. Australia.
An intractability result for multiple integration I.H.Sloan y and H.Wozniakowski yy Columbia University Computer Science Technical Report CUCS-019-96 y School of Mathematics University of ew South Wales
More informationWhat s new in high-dimensional integration? designing for applications
What s new in high-dimensional integration? designing for applications Ian H. Sloan i.sloan@unsw.edu.au The University of New South Wales UNSW Australia ANZIAM NSW/ACT, November, 2015 The theme High dimensional
More informationarxiv: v1 [math.gr] 3 Feb 2019
Galois groups of symmetric sextic trinomials arxiv:1902.00965v1 [math.gr] Feb 2019 Alberto Cavallo Max Planck Institute for Mathematics, Bonn 5111, Germany cavallo@mpim-bonn.mpg.de Abstract We compute
More informationIterative algorithms based on the hybrid steepest descent method for the split feasibility problem
Available online at www.tjnsa.com J. Nonlinear Sci. Appl. 9 (206), 424 4225 Research Article Iterative algorithms based on the hybrid steepest descent method for the split feasibility problem Jong Soo
More information1 Math 241A-B Homework Problem List for F2015 and W2016
1 Math 241A-B Homework Problem List for F2015 W2016 1.1 Homework 1. Due Wednesday, October 7, 2015 Notation 1.1 Let U be any set, g be a positive function on U, Y be a normed space. For any f : U Y let
More informationarxiv: v1 [math.co] 7 Jul 2014
Sum-ratio estimates over arbitrary finite fields Oliver Roche-Newton arxiv:1407.1654v1 [math.co] 7 Jul 2014 July 11, 2018 Abstract The aim of this note is to record a proof that the estimate max{ A+A,
More informationCourse 212: Academic Year Section 1: Metric Spaces
Course 212: Academic Year 1991-2 Section 1: Metric Spaces D. R. Wilkins Contents 1 Metric Spaces 3 1.1 Distance Functions and Metric Spaces............. 3 1.2 Convergence and Continuity in Metric Spaces.........
More information2. The Concept of Convergence: Ultrafilters and Nets
2. The Concept of Convergence: Ultrafilters and Nets NOTE: AS OF 2008, SOME OF THIS STUFF IS A BIT OUT- DATED AND HAS A FEW TYPOS. I WILL REVISE THIS MATE- RIAL SOMETIME. In this lecture we discuss two
More informationFoSP. FoSP. WARING S PROBLEM WITH DIGITAL RESTRICTIONS IN F q [X] Manfred Madritsch. Institut für Analysis und Computational Number Theory (Math A)
FoSP Algorithmen & mathematische Modellierung FoSP Forschungsschwerpunkt Algorithmen und mathematische Modellierung WARING S PROBLM WITH DIGITAL RSTRICTIONS IN F q [X] Manfred Madritsch Institut für Analysis
More informationOn Systems of Diagonal Forms II
On Systems of Diagonal Forms II Michael P Knapp 1 Introduction In a recent paper [8], we considered the system F of homogeneous additive forms F 1 (x) = a 11 x k 1 1 + + a 1s x k 1 s F R (x) = a R1 x k
More informationReliable Error Estimation for Quasi-Monte Carlo Methods
Reliable Error Estimation for Quasi-Monte Carlo Methods Department of Applied Mathematics, Joint work with Fred J. Hickernell July 8, 2014 Outline 1 Motivation Reasons Problem 2 New approach Set-up of
More informationMETRIC HEIGHTS ON AN ABELIAN GROUP
ROCKY MOUNTAIN JOURNAL OF MATHEMATICS Volume 44, Number 6, 2014 METRIC HEIGHTS ON AN ABELIAN GROUP CHARLES L. SAMUELS ABSTRACT. Suppose mα) denotes the Mahler measure of the non-zero algebraic number α.
More informationSMOOTHNESS PROPERTIES OF SOLUTIONS OF CAPUTO- TYPE FRACTIONAL DIFFERENTIAL EQUATIONS. Kai Diethelm. Abstract
SMOOTHNESS PROPERTIES OF SOLUTIONS OF CAPUTO- TYPE FRACTIONAL DIFFERENTIAL EQUATIONS Kai Diethelm Abstract Dedicated to Prof. Michele Caputo on the occasion of his 8th birthday We consider ordinary fractional
More information#A69 INTEGERS 13 (2013) OPTIMAL PRIMITIVE SETS WITH RESTRICTED PRIMES
#A69 INTEGERS 3 (203) OPTIMAL PRIMITIVE SETS WITH RESTRICTED PRIMES William D. Banks Department of Mathematics, University of Missouri, Columbia, Missouri bankswd@missouri.edu Greg Martin Department of
More informationOn components of vectorial permutations of F n q
On components of vectorial permutations of F n q Nurdagül Anbar 1, Canan Kaşıkcı 2, Alev Topuzoğlu 2 1 Johannes Kepler University, Altenbergerstrasse 69, 4040-Linz, Austria Email: nurdagulanbar2@gmail.com
More informationMATH 205C: STATIONARY PHASE LEMMA
MATH 205C: STATIONARY PHASE LEMMA For ω, consider an integral of the form I(ω) = e iωf(x) u(x) dx, where u Cc (R n ) complex valued, with support in a compact set K, and f C (R n ) real valued. Thus, I(ω)
More informationELLIPTIC CURVES SEMINAR: SIEGEL S THEOREM
ELLIPTIC CURVES SEMINAR: SIEGEL S THEOREM EVAN WARNER 1. Siegel s Theorem over Q 1.1. Statement of theorems. Siegel s theorem, in its simplest form, is the fact that a nonsingular elliptic curve contains
More informationM. VAN BAREL Department of Computing Science, K.U.Leuven, Celestijnenlaan 200A, B-3001 Heverlee, Belgium
MATRIX RATIONAL INTERPOLATION WITH POLES AS INTERPOLATION POINTS M. VAN BAREL Department of Computing Science, K.U.Leuven, Celestijnenlaan 200A, B-3001 Heverlee, Belgium B. BECKERMANN Institut für Angewandte
More information= 1 2x. x 2 a ) 0 (mod p n ), (x 2 + 2a + a2. x a ) 2
8. p-adic numbers 8.1. Motivation: Solving x 2 a (mod p n ). Take an odd prime p, and ( an) integer a coprime to p. Then, as we know, x 2 a (mod p) has a solution x Z iff = 1. In this case we can suppose
More informationGeneralizing the Tolerance Function for Guaranteed Algorithms
Generalizing the Tolerance Function for Guaranteed Algorithms Lluís Antoni Jiménez Rugama Joint work with Professor Fred J. Hickernell Room 120, Bldg E1, Department of Applied Mathematics Email: lluisantoni@gmail.com
More informationCourse 311: Michaelmas Term 2005 Part III: Topics in Commutative Algebra
Course 311: Michaelmas Term 2005 Part III: Topics in Commutative Algebra D. R. Wilkins Contents 3 Topics in Commutative Algebra 2 3.1 Rings and Fields......................... 2 3.2 Ideals...............................
More informationSYMMETRY AND SPECIALIZABILITY IN THE CONTINUED FRACTION EXPANSIONS OF SOME INFINITE PRODUCTS
SYMMETRY AND SPECIALIZABILITY IN THE CONTINUED FRACTION EXPANSIONS OF SOME INFINITE PRODUCTS J MC LAUGHLIN Abstract Let fx Z[x] Set f 0x = x and for n 1 define f nx = ff n 1x We describe several infinite
More informationNOTES ON ZHANG S PRIME GAPS PAPER
NOTES ON ZHANG S PRIME GAPS PAPER TERENCE TAO. Zhang s results For any natural number H, let P (H) denote the assertion that there are infinitely many pairs of distinct primes p, q with p q H; thus for
More information7: FOURIER SERIES STEVEN HEILMAN
7: FOURIER SERIES STEVE HEILMA Contents 1. Review 1 2. Introduction 1 3. Periodic Functions 2 4. Inner Products on Periodic Functions 3 5. Trigonometric Polynomials 5 6. Periodic Convolutions 7 7. Fourier
More informationCover Page. The handle holds various files of this Leiden University dissertation
Cover Page The handle http://hdl.handle.net/1887/32076 holds various files of this Leiden University dissertation Author: Junjiang Liu Title: On p-adic decomposable form inequalities Issue Date: 2015-03-05
More informationGEOMETRIC CONSTRUCTIONS AND ALGEBRAIC FIELD EXTENSIONS
GEOMETRIC CONSTRUCTIONS AND ALGEBRAIC FIELD EXTENSIONS JENNY WANG Abstract. In this paper, we study field extensions obtained by polynomial rings and maximal ideals in order to determine whether solutions
More informationMeasurable functions are approximately nice, even if look terrible.
Tel Aviv University, 2015 Functions of real variables 74 7 Approximation 7a A terrible integrable function........... 74 7b Approximation of sets................ 76 7c Approximation of functions............
More informationAdapting quasi-monte Carlo methods to simulation problems in weighted Korobov spaces
Adapting quasi-monte Carlo methods to simulation problems in weighted Korobov spaces Christian Irrgeher joint work with G. Leobacher RICAM Special Semester Workshop 1 Uniform distribution and quasi-monte
More informationHot new directions for QMC research in step with applications
Hot new directions for QMC research in step with applications Frances Kuo f.kuo@unsw.edu.au University of New South Wales, Sydney, Australia Frances Kuo @ UNSW Australia p.1 Outline Motivating example
More informationEXPLICIT EVALUATIONS OF SOME WEIL SUMS. 1. Introduction In this article we will explicitly evaluate exponential sums of the form
EXPLICIT EVALUATIONS OF SOME WEIL SUMS ROBERT S. COULTER 1. Introduction In this article we will explicitly evaluate exponential sums of the form χax p α +1 ) where χ is a non-trivial additive character
More informationA proof of a partition conjecture of Bateman and Erdős
proof of a partition conjecture of Bateman and Erdős Jason P. Bell Department of Mathematics University of California, San Diego La Jolla C, 92093-0112. US jbell@math.ucsd.edu 1 Proposed Running Head:
More information