Dyadic diaphony of digital sequences

Size: px
Start display at page:

Download "Dyadic diaphony of digital sequences"

Transcription

1 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 give formulae for the dyadic diaphony of digital, s)-sequences over Z, s, These formulae show that for fixed s {, }, the dyadic diaphony has the same values for any digital, s)- sequence For s, it follows that the dyadic diaphony and the diaphony of special digital, )-sequences are up to a constant the same We give the exact asymptotic order of the dyadic diaphony of digital, s)-sequences and show that for s it satisfies a central limit theorem Introduction The diaphony F see [3] or [4, Definition 9] or [8, Exercise 57, p 6]) of the first elements of a sequence ω x n ) n in [, ) s is given by / F ω) e πi k,xn ρk), k Z s k where for k k,, k s ) Z s it is ρk) s i max, k i ) and, denotes the usual inner product in R s It is well known that the diaphony is a quantitative measure for the irregularity of distribution of the first points of a sequence In fact, a sequence ω is uniformly distributed modulo if and only if lim F ω) Throughout this paper we will call the diaphony the classical diaphony In [7] Hellekalek and Leeb introduced the notion of dyadic diaphony which is similar to the classical diaphony but with the trigonometric functions replaced by Walsh functions Before we give the exact definition of the dyadic diaphony recall that Walsh-functions in base can be defined as follows: for a non-negative integer k with base representation k κ m m + + κ + κ and a real x with canonical) base representation x x + x + the k-th Walsh function in base is defined as wal k x) : ) x κ +x κ + +x m+ κ m For dimension s, x,, x s [, ) and k,, k s we define s wal k,,k s x,, x s ) : wal kj x j ) The author is supported by the Austrian Science Foundation FWF), Project S969, that is part of the Austrian ational Research etwork Analytic Combinatorics and Probabilistic umber Theory j

2 For vectors k k,, k s ) s and x x,, x s ) R s we write wal k x) : wal k,,k s x,, x s ) ow we can give the definition of the dyadic diaphony see Hellekalek and Leeb [7]) Definition The dyadic diaphony F, of the first elements of a sequence ω x n ) n in [, ) s is defined by F, ω) 3 s k s k ψk) wal k x n ) where for k k,, k s ) it is ψk) s i ψk i) and for k, { if k, ψk) r if r k < r+ with r Throughout the paper we will write rk) r if r is the unique determined integer such that r k < r+ It is shown in [7, Theorem 3] that the dyadic diaphony is a quantitative measure for the irregularity of distribution of the first points of a sequence: a sequence ω is uniformly distributed modulo if and only if lim F, ω) Further it was shown in [] that the dyadic diaphony is up to a factor depending only on s the worst-case error for quasi-monte Carlo integration of functions from a certain Hilbert space We consider the dyadic diaphony of a special class of sequences in [, ) s, namely of so-called digital, s)-sequences over Z for s, Digital, s)-sequences or more generally digital t, s)-sequences were introduced by iederreiter [, ] and they provide at the moment the most efficient method to generate sequences with excellent distribution properties We remark that a digital, s)-sequence over Z only exists if s or s For higher dimensions s 3 the concept of digital t, s)-sequence over Z with t > has to be stressed see [] or []) Before we give the definition of digital, s)-sequences we introduce some notation: for a vector c γ, γ, ) Z and for m we denote the vector in Z m consisting of the first m components of c by cm), ie, cm) γ,, γ m ) Further for an matrix C over Z and for m we denote by Cm) the left upper m m submatrix of C Definition For s {, }, choose s matrices C,, C s over Z with the following property: for every m and every n m the vectors c ) m),, c n ) s) m), c m),, c m nm) s) are linearly independent in Z m j) Here c i is the i-th row vector of the matrix C j In particular for any m the matrix C j m) has full rank over Z for all j {,, s}) For n let n n + n + n + be the base representation of n For j {,, s} multiply the vector n n, n, ) with the matrix C j, C j n : x j n ), xj n ), ) Z /,

3 and set Finally set x n : x ) n,, x s) n ) x n j) : xj n ) + xj n ) + Every sequence x n ) n constructed in this way is called digital, s)-sequence over Z The matrices C,, C s are called the generator matrices of the sequence To guarantee that the points x n belong to [, ) s and not just to [, ] s ) and also for the analysis of the sequence we need the condition that for each n and j s, we have x j ni) for infinitely many i This condition is always satisfied if we assume that for each j s and r we have c j i,r for all sufficiently large i, where cj i,r are the entries of the matrix C j Throughout this paper we assume that the generator matrices fulfill this condition see [, p7] where this condition is called S6)) For example if s and if we choose as generator matrix the identity matrix, then the resulting digital, )-sequence over Z is the well known van der Corput sequence in base Hence the concept of digital, )-sequences over Z is a generalization of the construction principle of the van der Corput sequence For the classical diaphony it was proved by Faure [5] that F ω)) π u u, ) if ω is a digital, )-sequence over Z whose generator matrix C is a non-singular upper triangular matrix Faure and we shall do so as well) called these sequences UT-sequences Here denotes the distance to the nearest integer function, ie, x : minx x, x x )) See also [, 6, ] for further results concerning the classical diaphony of special -dimensional sequences The aim of this paper is to prove a similar formula for the dyadic diaphony of digital, s)-sequences over Z for s {, } see Theorems and ) These formulae show that for fixed s the dyadic diaphony is invariant for digital, s)-sequences over Z Further we find that the dyadic diaphony and the classical diaphony of UT-sequences s ) only differ by a multiplicative constant Corollary ) We obtain the exact asymptotic order of the dyadic diaphony of digital, s)-sequences over Z Corollary and 4) Moreover it follows from our formula that the squared dyadic diaphony of digital, )-sequences over Z satisfies a central limit theorem Corollary 3) For digital, )-sequences we will obtain a similar, but weaker result Corollary 5) The Results for s First we give the formula for the dyadic diaphony of digital, )-sequences over Z This formula shows that the dyadic diaphony is invariant for digital, )-sequences over Z Theorem Let ω be a digital, )-sequence over Z Then for any we have F, ω)) 3 u u 3

4 We defer the proof of this formula to Section 4 Remark In Theorem we have an infinite sum for the dyadic diaphony of a digital, )-sequence over Z This formula can easily be made computable since for m we have / u / u for u m + Therefore we have m ) F, ω)) 3 + u From Theorem we find the surprising result that the classical diaphony and the dyadic diaphony of a UT-sequence are essentially the same Corollary Let ω be a digital UT-sequence over Z Then for any we have u m F, ω) 3 π F ω) Proof This follows from Theorem together with Faures formula ) From Remark one can see immediately that the dyadic diaphony of a digital, )- sequence over Z is of order F, ω) O log /) But we can even be much more precise From a thorough analysis of the sum in Remark we obtain the exact dependence of the dyadic diaphony of digital, )-sequences over Z on log / Corollary Let ω be a digital, )-sequence over Z For m we have and F, ω)) m )m 9 m 9 m lim sup F, ω)) log 3 log The proof of this result will be given in Section 5 We just remark that the result for the lim sup follows also from a result of Chaix and Faure [, Théoréme 43] for the classical diaphony of the van der Corput sequence together with Corollary and Theorem In [3] the authors showed that the star discrepancy and all L p -discrepancies of the van der Corput sequence in base satisfy a central limit theorem The same arguments as in their proofs can now be used to obtain the subsequent result Corollary 3 Let ω be a digital, )-sequence over Z Then for every real y we have { M # < M : F, ω)) 4 log + y } 4 log Φy) + o), 3 where Φy) π y e t dt denotes the normal distribution function and log denotes the logarithm to the base Ie, the squared dyadic diaphony of a digital, )-sequence over Z satisfies a central limit theorem 4

5 Remark Together with Corollary we also obtain a central limit theorem for the square of the classical diaphony of UT-sequences Proof As already mentioned, the proof follows exactly the lines of the proof of [3, Theorem ] One only has to compute the expectation and the variance of the random variable S m m X w, w where X is uniformly distributed on [, ) By tedious but straightforward calculations we obtain ES m m/ and VarS m m/ m )/6 3 The Results for s We give the formula for the dyadic diaphony of digital, )-sequences over Z which shows that the dyadic diaphony is invariant for digital, )-sequences as well Theorem Let ω be a digital, )-sequence over Z Then for any we have F, ω)) 9 4 u u u We defer the proof of this formula to Section 4 Remark 3 In Theorem we have an infinite sum for the dyadic diaphony of a digital, )-sequence over Z Again this formula can easily be made computable compare with Remark ) For m we have F, ω)) 9 m ) 4 + 3m 4 u u + m 4 u From Remark 3 one can see immediately that the dyadic diaphony of a digital, )- sequence over Z is of order F, ω) Olog /) Also here we obtain from a thorough analysis of the sum in Remark 3 the exact dependence of the dyadic diaphony of digital, )-sequences over Z on log / Corollary 4 Let ω be a digital, )-sequence over Z Then for any m we have and F, ω)) m 8 + 7m O m m ) lim sup F, ω)) log ) 8log ) The proof of this result will be given in Section 5 Following this proof the Om/ m )- term in the above bound can easily be made explicit Unfortunately we could not show that the squared dyadic diaphony of a digital, )- sequence over Z satisfies a central limit theorem However, we were able to prove the following result 5

6 Corollary 5 Let ω be a digital, )-sequence over Z Then for any ε > we have { ) } lim < m 3 : m 3 ε < F, ω) < 3 log 3 + ε Proof By tedious but straightforward calculations using Theorem we obtain and F, ω)) 3 3 m m + Om m ) m m From this the result immediately follows F, ω)) m4 m + Om 3 m ) 4 The Proofs of Theorems and For the proofs of Theorems and we need the subsequent lemma This result was implicitly proved in [3] For the sake of completeness we provide the short proof Lemma Let the non-negative integer U have binary expansion U U + U + + U m m For any non-negative integer n U let n n + n + + n m m be the binary representation of n For p m let Up) : U + + U p p Let b, b,, b m be arbitrary elements of Z, not all zero Then U ) b n + +b m n m ) b w+u w+ + +b m U m w+ U, where w is minimal such that b w Proof From splitting up the sum we obtain U ) b n + +b m n m w+ U w+ + +U m m w ) + Uw) ) nw ) b w+u w+ + +b m U m ) nw ) b w+n w+ + +b m n m Uw) + ) b w+u w+ + +b m U m ) nw w+ { ) b w+u w+ + +b m U m Uw) if Uw) < w, w+ Uw) if Uw) w, { ) b w+u w+ + +b m U m Uw)/ w+ w+ if Uw)/ w+ < /, Uw)/ w+ if Uw)/ w+ /, ) b w+u w+ + +b m U m w+ Uw) w+ 6

7 Since Uw) U w+ the result follows w+ ow we can give the Proof of Theorem Let r k < r+ Then k k + k + + k r r with k i {, }, i < r and k r Let, denote the usual inner product in Z and let c i Z be the i-th row vector of the generator matrix C of the digital, )-sequence for short we write C instead of C here) Since the i-th digit x n i) of the point x n, i, n, is given by c i, n see Definition ) we have wal k x n ) ) k c, n + +k r c r+, n ) k c + +k r c r+, n ) Let C c i,j ) i,j For k, k k + k + + k r r, k i {, }, i < r and k r define uk) : min{l : k c,l + +k r c r+,l } ote that since C generates a digital, )-sequence over Z we obviously have uk) r + For fixed k, r k < r+ let b b, b, ) : k c + + k r c r+ Let m m If uk) m we obtain from ) together with Lemma, wal k x n ) ) b, n ) n b + +n m b m ) n uk) + ) uk)b uk) + uk) uk) But if uk) > m the above equality is trivially true Therefore we have F, ω)) k r r uk) rk) r+ uk) r k r r uk) r+ u u u u u u ru uk) ) r+ k r uk)u r+ r k r uk)u 7

8 ow we have to evaluate the sum r+ for r u and u This is the number of k r uk)u vectors k,, k r ) Z r such that Cr + ) k k r x u+ x r+ Z r+ 3) for arbitrary x u+,, x r+ Z Recall that for an integer m we denote by Cm) the left upper m m submatrix of the matrix C, see Section ) We consider two cases: Assume that r u Then system 3) becomes k Cr + ) k r Since the r + ) r + ) matrix Cr + ) is regular over Z it is clear that there exists a vector k k,, k r ) Z r+ such that Cr + ) k,,, ) Assume that k r, then we have Cr) k,, k r ),, ) Again we know that Cr) is regular over Z and therefore we obtain k k r Hence k, the zero vector in Z r+ This is now a contradiction since k is a solution of the system Cr + ) k,,, ) Therefore we have u k u uk)u Assume that r u Since Cr) is regular over Z it is clear that Dr) : Cr) ) is regular over Z Hence for any vector k Z r there is a vector l Z r such that k Dr) l Therefore system 3) can be rewritten as Cr + ) Dr) l ) x u+ x r+ 8

9 with l Z r ow we use the definition of the matrix Dr) and find that the above system is equivalent to the system l + c r+ r + ), l r x u+ d d d r d r x r+ where d,, d r ) : c,r+,, c r,r+ )Dr) ow one can easily see that for arbitrary x u+,, x r there exists exactly one solution l l,, l r ) Z r such that the first r lines of the above system are fulfilled Further there is exactly one possible choice of x r+ Z such that this vector l is a solution of the above system Therefore we obtain r+ k r uk)u r u ow we have F, ω)) The result follows u u ) u + u ) r u 6 r ru u u Proof of Theorem Let ω x n ) n be a digital, )-sequence over Z Let x n x n, y n ) for n Clearly the sequences x n ) n and y n ) n are digital, )-sequences over Z We have F, ω)) 8 k k ψk) wal k x n ) wal 8 rk) k x n ) + wal 8 rl) l y n ) k l + wal 8 rk)+rl) k x n )wal l y n ) k,l 3 u + wal 8 rk)+rl) k x n )wal l y n ), 4) u k,l where for the last equality we used Theorem We have to consider Σ : wal rk)+rl) k x n )wal l y n ) k,l 9

10 Assume that r k < r+ and t l < t+ Then k k + k + + k r r with k i {, }, i < r and k r and l l + l + + l t t with l j {, }, j < t and l t Let c i Z be the i-th row vector of the generator matrix C and let d i Z be the i-th row vector of the generator matrix C, i Since the i-th digit x n i) of x n is given by c i, n and the i-th digit y n i) of y n is given by d i, n see Definition ) we have wal k x n )wal l y n ) ) k c, n + +k r c r+, n +l d, n + +l t d t+, n ) k c + +k r c r+ +l d + +l tdt+, n Let C c i,j ) i,j and C d i,j ) i,j Define uk, l) : min{j : k c,j + + k r c r+,j + l d,j + + l s d s+,j } Since C, C generate a digital, )-sequence over Z we obviously have uk, l) r+s+ As in the proof of Theorem we now apply Lemma and obtain wal k x n )wal l y n ) uk,l) Therefore we have Σ k,l r,t r,t rk)+rl) uk,l) r+t We have to evaluate the double-sum r+ k r u t+ uk,l) l t uk,l) r+t+ u r+t u r+ k r t+ l t } {{ } uk,l)u number of k,, k r, l,, l t Z such that c, c r+, d, d t+, c, c r+, d, d t+, c,r+t+ c r+,r+t+ d,r+t+ d t+,r+t+ uk,l) r+ k r uk,l) t+ l t } {{ } uk,l)u for u r + t + This is the k k r l l t x u+ x r+t+ 5)

11 for arbitrary x u+,, x r+t+ Z ote that the above matrix we denote it by Cr, t) is an r + t + ) r + t + ) matrix Since C, C generate a digital, )-sequence over Z, it follows that Cr, t) is regular ow we consider three cases: Assume that u r + t + Then system 5) becomes Cr, t) h 6) Since Cr, t) is regular there exists a vector h k,, k r, l,, l t ) Z r+z+, h, such that Cr, t) h,,, ) Assume that l t Then But then c, c r+, d, d t, c, c r+, d, d t, c,r+t+ c r+,r+t+ d,r+t+ d t,r+t+ c, c r+, d, d t, c, c r+, d, d t, c,r+t+ c r+,r+t+ d,r+t+ d t,r+t+ k k r l l t k k r l l t Since the above matrix is again regular we obtain that k,, k r, l,, l t ),, ) and therefore h, a contradiction Hence l t and in the same way one can show that k r We have shown that system 6) has exactly one solution and therefore we have r+ k r s+ l s } {{ } uk,l)u Assume that u r + t + Let x Z Since Cr, t) is regular there exists exactly one vector h Z r+t+ such that Cr, t) h,,,, x) Assume that h is of the form h k,, k r,, l,, l t, ) particular we have h c, c r+, d, d t+, c, c r+, d, d t+, c,r+t c r+,r+t d,r+t d t+,r+t Z r+t+ In 7)

12 Since c, c r, d, d t, c, c r, d, d t, c,r+t c r,r+t d,r+t d t,r+t is regular we find that h is the unique solution of 7) Hence h is exactly the same vector as in the first case where u r + t + But then h cannot be a solution of Cr, t) h,,,, x) Therefore we obtain r+ k r s+ l s } {{ } uk,l)u 3 Assume that u r + t We rewrite system 5) in the form c, c r, d, d t, c, c r, d, d t, c,r+t c r,r+t d,r+t d t,r+t c,r+t+ c r,r+t+ d,r+t+ d t,r+t+ c,r+t+ c r,r+t+ d,r+t+ d t,r+t+ k k r l l t x u+ x r+t x r+t+ x r+t+ + y r,t where y r,t c r+, + d t+,,, c r+,r+t+ + d t+,r+t+ ) Z r+t+ Since the upper r + t) r + t) sub-matrix of the above matrix is regular we find for arbitrary x u+,, x r+t exactly one solution of the first r + t rows of the above system This solution can be made a solution of the whole system by an adequate choice of x r+t+ and x r+t+ Therefore we have r+ k r s+ l s } {{ } uk,l)u r+t u

13 ow we have Σ r+t u r+t u r,t r,t r+t u r+t u u u u u r,t r+t u u u u wu u u + 4) + 6 8u ) u u u u + 6 r+t u + r+t+) r,t r+t+ + 6 r+t u w + w + 6 u r,t r+tu u u ) u u ) u The result follows by inserting this expression into 4) u r+t+ ) 5 The Proofs of Corollaries and 4 We will say that a real β in [, ) is m-bit if β b + + bm m with b i {, } Ie, an m-bit number is of the form k/ m with k {,,, m } The essential technical tool for the proof of Corollary is provided by Lemma Assume that β, b b this here and in the following always means base representation) has two equal consecutive digits b i b i+ with i m and let i be minimal with this property, ie, Replace β by Then m u u γ m u u β + β, b i+ or β, b i+ or β, b i+ or β, b i+ γ, b i+ resp γ, b i+ resp γ, b i+ resp γ, b i+ 9 9 ) )i τ) in the first two cases, i ) )i τ in the last two cases, i 3

14 where τ :, b i+ b i+3 Remark 4 In any case we have m u u γ m u u β with equality iff τ in the first two cases and iff τ in the last two cases Proof of Lemma This is simple calculation We just handle the first case here m i u γ u β ) γ i β + u γ) u β ) 8) u Here γ 3 u γ) 3 u ) + τ and i β τ Further, for even u we have i i+ ) ) u+ u+ i + τ i u and u β 3 i + τ i+ u, and for odd u we have u γ) ) + u+ τ and u β ) + u+ τ 3 i i u 3 i i+ u Inserting this into 8) we obtain m u γ u β ) 9 u The other cases are calculated in the same way + i ) τ) From Lemma we obtain the subsequent result concerning the maximum of m u u β over all β We remark that in [9] the authors considered the same problem without the square at the -function Lemma 3 Consider β R with the canonical base representation ie, with infinitely many digits equal to zero) Then there exists m max β u u β m )m 7 m 7 m and it is attained if and only if β is of the form β with β ) ) m+ or β ) m ) 3 3 Remark 5 ote that 3 ) ) m+ {, if m is odd,, if m is even, and ) m ) {, if m is odd, 3, if m is even 4

15 Proof of Lemma 3 For any γ, c c c m c m+ with fixed c,, c m the sum m u u γ obviously becomes maximal if c m and c m+ c m+, or if c m and c m+ c m+ Hence by Lemma the u β m sup β u only can be attained, respectively approached by if m is even, and by if m is odd ow we check easily that m u for k,, 6 and the result follows β, or b m is the last zero) β, or β 3, b m is the last one) β 4, or b m is the last zero) β 5, or β 6, b m is the last one) u β k m )m 7 m 7 m We give the Proof of Corollary We have max m m u u m max u β m β m bit )m 7 m 7 m by Lemma 3 The result follows now together with Remark u For the proof of Corollary 4 we can in principle proceed as for the proof of Corollary However, in this case the detailed computations are by far more involved than above First we have Lemma 4 Assume that β, b b has two equal consecutive digits b i b i+ with i m and let i be minimal with this property, ie, β, b i+ or β, b i+ or β, b i+ or β, b i+ 5

16 Replace β by γ, b i+ resp γ, b i+ resp γ, b i+ resp γ, b i+ Then m u u γ m u) + m 9 m 9 where τ :, b i+ b i+3 m u u β m u) ) )i i + 4 ) i)) τ) in the first two cases, i 9 7 i i ) )i i + 4 ) i)) τ in the last two cases, i 9 7 i i Remark 6 In any case, for m > 3, we have m u u γ m u) m u u β m u) Proof of Lemma 4 We have m u γ u β )m u) u i m γ m i) i β + u γ) m u ) u β m u) ) u The result now follows as in the proof of Lemma by some tedious but straightforward algebra With Lemma 4 we obtain Lemma 5 We have m max β m bit u { m u β m u) + m ) ) m 7 m + m 3 8 m + m m + ) )) m 7 7 m + 4 m 3 For even m the maximum is attained if and only if {, β 3 + ), 3 m m ) For odd m the maximum is attained if and only if { ), β 3 + m ), 3 m+ 8 m ) and and if m is even, if m is odd 6

17 Proof For short we write f m β) : m u u β m u) Let m be even It follows from Lemma 4 that the m-bit number β which maximizes our sum has to be of the form β, b m b m or β, b m b m First we deal with β, b m b m ow we consider four cases corresponding to the possible choices for b m and b m If b m, b m ), ), then f m β ) m 8 + m If b m, b m ), ), then f m β ) m 8 + m If b m, b m ), ), then f m β ) m 8 + m If b m, b m ), ), then f m β ) m 8 + m m 6 m 6 m 7 m 8 64 m 8 m m 5 m + m 7 m 8 m 8 m m 4 m m m 8 6 m 8 m m m m 7 4 m 8 4 m 8 m Therefore we find that the choice b m b m ), ) gives the maximal value, ie, β, For β, b m b m we find in the same way that b m, b m ), ) gives the maximal value, ie, β, Since f m β ) f m β ) m 8 + m m 4 7 m ) m + 3 ) 8 ) 8 m the result follows for even m For odd m 3 the result can be proved analogously We give the Proof of Corollary 4 We have max m m u u m u max u β m u) m β m bit 8 + m m ) 8 + O m u by Lemma 5 The result follows now together with Remark 3 Acknowledgement The author wishes to thank Peter Kritzer for his many comments and suggestions 7

18 References [] Chaix, H, Faure, H, Discrépance et diaphonie en dimension un Acta Arith 63: 3-4, 993 [] Dick, J, Pillichshammer, F, Diaphony, discrepancy, spectral test and worst-case error Math Comput Simulation 7: 59 7, 5 [3] Drmota, M, Larcher, G, Pillichshammer, F, Precise distribution properties of the van der Corput sequence and related sequences Manuscripta Math 8: 4, 5 [4] Drmota, M, Tichy, RF, Sequences, Discrepancies and Applications Lecture otes in Mathematics 65, Springer-Verlag, Berlin, 997 [5] Faure, H, Discrepancy and Diaphony of digital, )-sequences in prime base Acta Arith 7: 5 48, 4 [6] Grozdanov, V, On the diaphony of one class of one-dimensional sequences Internat J Math Math Sci 9: 5 4, 996 [7] Hellekalek, P, Leeb, H, Dyadic diaphony Acta Arith 8: 87 96, 997 [8] Kuipers, L, iederreiter, H, Uniform Distribution of Sequences John Wiley, ew York, 974 [9] Larcher, G, Pillichshammer, F, Sums of distances to the nearest integer and the discrepancy of digital nets Acta Arith 6: , 3 [] iederreiter, H, Point sets and sequences with small discrepancy Monatsh Math 4: , 987 [] iederreiter, H, Random umber Generation and Quasi-Monte Carlo Methods o 63 in CBMS-SF Series in Applied Mathematics SIAM, Philadelphia, 99 [] Proinov, PD, Grozdanov, VS, On the diaphony of the van der Corput-Halton sequence J umber Theory 3: 94 4, 988 [3] Zinterhof, P, Über einige Abschätzungen bei der Approximation von Funktionen mit Gleichverteilungsmethoden Sitzungsber Österr Akkad Wiss Math-atur Kl II 85: 3, 976 Author s Address: Friedrich Pillichshammer, Institut für Finanzmathematik, Universität Linz, Altenbergstraße 69, A-44 Linz, Austria friedrichpillichshammer@jkuat 8

Uniform distribution of sequences connected with the weighted sum-of-digits function

Uniform 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 information

f(x n ) [0,1[ s f(x) dx.

f(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 information

uniform distribution theory

uniform 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 information

APPLIED MATHEMATICS REPORT AMR04/2 DIAPHONY, DISCREPANCY, SPECTRAL TEST AND WORST-CASE ERROR. J. Dick and F. Pillichshammer

APPLIED 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 information

A NEW UPPER BOUND ON THE STAR DISCREPANCY OF (0,1)-SEQUENCES. Peter Kritzer 1. Abstract

A 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 information

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 ϱ α. 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 information

L 2 Discrepancy of Two-Dimensional Digitally Shifted Hammersley Point Sets in Base b

L 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 information

Construction algorithms for plane nets in base b

Construction 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 information

On the inverse of the discrepancy for infinite dimensional infinite sequences

On the inverse of the discrepancy for infinite dimensional infinite sequences On the inverse of the discrepancy for infinite dimensional infinite sequences Christoph Aistleitner Abstract In 2001 Heinrich, ovak, Wasilkowski and Woźniakowski proved the upper bound s,ε) c abs sε 2

More information

Discrepancy Bounds for Nets and Sequences

Discrepancy 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 information

The Halton sequence and its discrepancy in the Cantor expansion

The Halton sequence and its discrepancy in the Cantor expansion Period Math Hung 2017 75:128 141 DOI 10.1007/s10998-016-0169-5 The Halton sequence and its discrepancy in the Cantor expansion Alena Haddley 1 Poj Lertchoosakul 2 Radhakrishnan air 1 Published online:

More information

arxiv: v1 [math.co] 8 Feb 2013

arxiv: v1 [math.co] 8 Feb 2013 ormal numbers and normality measure Christoph Aistleitner arxiv:302.99v [math.co] 8 Feb 203 Abstract The normality measure has been introduced by Mauduit and Sárközy in order to describe the pseudorandomness

More information

High Dimensional Integration: New Weapons Fighting the Curse of Dimensionality

High Dimensional Integration: New Weapons Fighting the Curse of Dimensionality High Dimensional Integration: ew Weapons Fighting the Curse of Dimensionality Peter Zinterhof Technical Report 26-4 October 26 Department of Computer Sciences Jakob-Haringer-Straße 2 52 Salzburg Austria

More information

DIGITAL EXPANSION OF EXPONENTIAL SEQUENCES

DIGITAL EXPANSION OF EXPONENTIAL SEQUENCES DIGITAL EXPASIO OF EXPOETIAL SEQUECES MICHAEL FUCHS Abstract. We consider the q-ary digital expansion of the first terms of an exponential sequence a n. Using a result due to Kiss und Tichy [8], we prove

More information

Star discrepancy of generalized two-dimensional Hammersley point sets

Star discrepancy of generalized two-dimensional Hammersley point sets Star discrepancy of generalized two-dimensional Hammersley point sets Henri Faure Astract We generalize to aritrary ases recent results on the star discrepancy of digitally shifted two-dimensional Hammersley

More information

A good permutation for one-dimensional diaphony

A good permutation for one-dimensional diaphony Monte Carlo Methods Appl. 16 (2010), 307 322 DOI 10.1515/MCMA.2010.015 de Gruyter 2010 A good permutation for one-dimensional diaphony Florian Pausinger and Wolfgang Ch. Schmid Astract. In this article

More information

Constructions of digital nets using global function fields

Constructions 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 information

arxiv: v2 [math.nt] 27 Mar 2018

arxiv: 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 information

Low-discrepancy sequences obtained from algebraic function fields over finite fields

Low-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 information

THE DISCREPANCY OF GENERALIZED VAN-DER-CORPUT-HALTON SEQUENCES

THE DISCREPANCY OF GENERALIZED VAN-DER-CORPUT-HALTON SEQUENCES THE DISCREPANCY OF GENERALIZED VAN-DER-CORPUT-HALTON SEQUENCES MICHAEL DRMOTA Abstract. The purpose of this paper is to provide upper bounds for the discrepancy of generalized Van-der-Corput-Halton sequences

More information

Tutorial on quasi-monte Carlo methods

Tutorial 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 information

Incomplete exponential sums over finite fields and their applications to new inversive pseudorandom number generators

Incomplete 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 information

Exponential Convergence and Tractability of Multivariate Integration for Korobov Spaces

Exponential 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 information

On the maximal density of sum-free sets

On the maximal density of sum-free sets ACTA ARITHMETICA XCV.3 (2000) On the maximal density of sum-free sets by Tomasz Luczak (Poznań) and Tomasz Schoen (Kiel and Poznań) 1. Introduction. For a set A N, let A(n) = A {1,..., n} and { } P(A)

More information

LS-sequences of points in the unit square

LS-sequences of points in the unit square LS-sequences of points in the unit square arxiv:2.294v [math.t] 3 ov 202 Ingrid Carbone, Maria Rita Iacò, Aljoša Volčič Abstract We define a countable family of sequences of points in the unit square:

More information

On the Classification of LS-Sequences

On the Classification of LS-Sequences On the Classification of LS-Sequences Christian Weiß arxiv:1706.08949v3 [math.nt] 4 Dec 2017 February 9, 2018 Abstract This paper addresses the question whether the LS-sequences constructed in [Car12]

More information

On the discrepancy estimate of normal numbers

On 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 information

Construction Algorithms for Higher Order Polynomial Lattice Rules

Construction Algorithms for Higher Order Polynomial Lattice Rules 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

More information

On families of anticommuting matrices

On families of anticommuting matrices On families of anticommuting matrices Pavel Hrubeš December 18, 214 Abstract Let e 1,..., e k be complex n n matrices such that e ie j = e je i whenever i j. We conjecture that rk(e 2 1) + rk(e 2 2) +

More information

Uniformly distributed sequences of partitions

Uniformly distributed sequences of partitions Uniformly distributed sequences of partitions Aljoša Volčič Università della Calabria Uniform distribution and QMC methods A. Volčič (Università della Calabria) U.d. sequences of partitions Linz, October

More information

Construction of digital nets from BCH-codes

Construction of digital nets from BCH-codes Construction of digital nets from BCH-codes Yves Edel Jürgen Bierbrauer Abstract We establish a link between the theory of error-correcting codes and the theory of (t, m, s)-nets. This leads to the fundamental

More information

A LeVeque-type lower bound for discrepancy

A LeVeque-type lower bound for discrepancy reprinted from Monte Carlo and Quasi-Monte Carlo Methods 998, H. Niederreiter and J. Spanier, eds., Springer-Verlag, 000, pp. 448-458. A LeVeque-type lower bound for discrepancy Francis Edward Su Department

More information

Van der Corput sets with respect to compact groups

Van der Corput sets with respect to compact groups Van der Corput sets with respect to compact groups Michael Kelly and Thái Hoàng Lê Abstract. We study the notion of van der Corput sets with respect to general compact groups. Mathematics Subject Classification

More information

Improved Discrepancy Bounds for Hybrid Sequences. Harald Niederreiter. RICAM Linz and University of Salzburg

Improved Discrepancy Bounds for Hybrid Sequences. Harald Niederreiter. RICAM Linz and University of Salzburg Improved Discrepancy Bounds for Hybrid Sequences Harald Niederreiter RICAM Linz and University of Salzburg MC vs. QMC methods Hybrid sequences The basic sequences Deterministic discrepancy bounds The proof

More information

A class of trees and its Wiener index.

A class of trees and its Wiener index. A class of trees and its Wiener index. Stephan G. Wagner Department of Mathematics Graz University of Technology Steyrergasse 3, A-81 Graz, Austria wagner@finanz.math.tu-graz.ac.at Abstract In this paper,

More information

The Distribution of Generalized Sum-of-Digits Functions in Residue Classes

The Distribution of Generalized Sum-of-Digits Functions in Residue Classes Journal of umber Theory 79, 9426 (999) Article ID jnth.999.2424, available online at httpwww.idealibrary.com on The Distribution of Generalized Sum-of-Digits Functions in Residue Classes Abigail Hoit Department

More information

ON THE CLASSIFICATION OF LS-SEQUENCES

ON THE CLASSIFICATION OF LS-SEQUENCES uniform distribution theory DOI: 10.2478/udt-2018 0012 Unif. Distrib. Theory 13 2018, no.2, 83 92 ON THE CLASSIFICATION OF LS-SEQUENCES Christian Weiß Hochschule Ruhr West, Mülheim an der Ruhr, GERMANY

More information

OPTIMAL SCALING FOR P -NORMS AND COMPONENTWISE DISTANCE TO SINGULARITY

OPTIMAL SCALING FOR P -NORMS AND COMPONENTWISE DISTANCE TO SINGULARITY published in IMA Journal of Numerical Analysis (IMAJNA), Vol. 23, 1-9, 23. OPTIMAL SCALING FOR P -NORMS AND COMPONENTWISE DISTANCE TO SINGULARITY SIEGFRIED M. RUMP Abstract. In this note we give lower

More information

NORMAL NUMBERS AND UNIFORM DISTRIBUTION (WEEKS 1-3) OPEN PROBLEMS IN NUMBER THEORY SPRING 2018, TEL AVIV UNIVERSITY

NORMAL NUMBERS AND UNIFORM DISTRIBUTION (WEEKS 1-3) OPEN PROBLEMS IN NUMBER THEORY SPRING 2018, TEL AVIV UNIVERSITY ORMAL UMBERS AD UIFORM DISTRIBUTIO WEEKS -3 OPE PROBLEMS I UMBER THEORY SPRIG 28, TEL AVIV UIVERSITY Contents.. ormal numbers.2. Un-natural examples 2.3. ormality and uniform distribution 2.4. Weyl s criterion

More information

Approximation of Functions Using Digital Nets

Approximation 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 information

arxiv:math/ v1 [math.nt] 27 Feb 2004

arxiv:math/ v1 [math.nt] 27 Feb 2004 arxiv:math/0402458v1 [math.nt] 27 Feb 2004 On simultaneous binary expansions of n and n 2 Giuseppe Melfi Université de Neuchâtel Groupe de Statistique Espace de l Europe 4, CH 2002 Neuchâtel, Switzerland

More information

c i r i i=1 r 1 = [1, 2] r 2 = [0, 1] r 3 = [3, 4].

c i r i i=1 r 1 = [1, 2] r 2 = [0, 1] r 3 = [3, 4]. Lecture Notes: Rank of a Matrix Yufei Tao Department of Computer Science and Engineering Chinese University of Hong Kong taoyf@cse.cuhk.edu.hk 1 Linear Independence Definition 1. Let r 1, r 2,..., r m

More information

THE SUM OF DIGITS OF PRIMES

THE SUM OF DIGITS OF PRIMES THE SUM OF DIGITS OF PRIMES Michael Drmota joint work with Christian Mauduit and Joël Rivat Institute of Discrete Mathematics and Geometry Vienna University of Technology michael.drmota@tuwien.ac.at www.dmg.tuwien.ac.at/drmota/

More information

Approximation of High-Dimensional Rank One Tensors

Approximation of High-Dimensional Rank One Tensors Approximation of High-Dimensional Rank One Tensors Markus Bachmayr, Wolfgang Dahmen, Ronald DeVore, and Lars Grasedyck March 14, 2013 Abstract Many real world problems are high-dimensional in that their

More information

Monte Carlo Method for Numerical Integration based on Sobol s Sequences

Monte Carlo Method for Numerical Integration based on Sobol s Sequences Monte Carlo Method for Numerical Integration based on Sobol s Sequences Ivan Dimov and Rayna Georgieva Institute of Information and Communication Technologies, Bulgarian Academy of Sciences, 1113 Sofia,

More information

Min-Rank Conjecture for Log-Depth Circuits

Min-Rank Conjecture for Log-Depth Circuits Min-Rank Conjecture for Log-Depth Circuits Stasys Jukna a,,1, Georg Schnitger b,1 a Institute of Mathematics and Computer Science, Akademijos 4, LT-80663 Vilnius, Lithuania b University of Frankfurt, Institut

More information

7 Asymptotics for Meromorphic Functions

7 Asymptotics for Meromorphic Functions Lecture G jacques@ucsd.edu 7 Asymptotics for Meromorphic Functions Hadamard s Theorem gives a broad description of the exponential growth of coefficients in power series, but the notion of exponential

More information

Solving a linear equation in a set of integers II

Solving a linear equation in a set of integers II ACTA ARITHMETICA LXXII.4 (1995) Solving a linear equation in a set of integers II by Imre Z. Ruzsa (Budapest) 1. Introduction. We continue the study of linear equations started in Part I of this paper.

More information

Lattices and Hermite normal form

Lattices and Hermite normal form Integer Points in Polyhedra Lattices and Hermite normal form Gennady Shmonin February 17, 2009 1 Lattices Let B = { } b 1,b 2,...,b k be a set of linearly independent vectors in n-dimensional Euclidean

More information

Estimates for probabilities of independent events and infinite series

Estimates for probabilities of independent events and infinite series Estimates for probabilities of independent events and infinite series Jürgen Grahl and Shahar evo September 9, 06 arxiv:609.0894v [math.pr] 8 Sep 06 Abstract This paper deals with finite or infinite sequences

More information

uniform distribution theory

uniform distribution theory Uniform Distribution Theory 6 (2011), no.1, 185 200 uniform distribution theory CONSTRUCTIONS OF UNIFORMLY DISTRIBUTED SEQUENCES USING THE b-adic METHOD Peter Hellekalek Harald Niederreiter ABSTRACT. For

More information

METRIC HEIGHTS ON AN ABELIAN GROUP

METRIC 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 information

Random and Quasi-Random Point Sets

Random and Quasi-Random Point Sets Peter Hellekalek Gerhard Larcher (Editors) Random and Quasi-Random Point Sets With contributions by JözsefBeck Peter Hellekalek Fred J. Hickemell Gerhard Larcher Pierre L'Ecuyer Harald Niederreiter Shu

More information

Notes on Equidistribution

Notes on Equidistribution otes on Equidistribution Jacques Verstraëte Department of Mathematics University of California, San Diego La Jolla, CA, 92093. E-mail: jacques@ucsd.edu. Introduction For a real number a we write {a} for

More information

THE SUM OF DIGITS OF n AND n 2

THE SUM OF DIGITS OF n AND n 2 THE SUM OF DIGITS OF n AND n 2 KEVIN G. HARE, SHANTA LAISHRAM, AND THOMAS STOLL Abstract. Let s q (n) denote the sum of the digits in the q-ary expansion of an integer n. In 2005, Melfi examined the structure

More information

A NOTE ON NORMAL NUMBERS IN MATRIX NUMBER SYSTEMS

A NOTE ON NORMAL NUMBERS IN MATRIX NUMBER SYSTEMS A NOTE ON NORMAL NUMBERS IN MATRIX NUMBER SYSTEMS MANFRED G. MADRITSCH Abstract. We consider a generalization of normal numbers to matrix number systems. In particular we show that the analogue of the

More information

Some results on the existence of t-all-or-nothing transforms over arbitrary alphabets

Some results on the existence of t-all-or-nothing transforms over arbitrary alphabets Some results on the existence of t-all-or-nothing transforms over arbitrary alphabets Navid Nasr Esfahani, Ian Goldberg and Douglas R. Stinson David R. Cheriton School of Computer Science University of

More information

Optimal compression of approximate Euclidean distances

Optimal compression of approximate Euclidean distances Optimal compression of approximate Euclidean distances Noga Alon 1 Bo az Klartag 2 Abstract Let X be a set of n points of norm at most 1 in the Euclidean space R k, and suppose ε > 0. An ε-distance sketch

More information

On Systems of Diagonal Forms II

On 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 information

MINIMAL GENERATING SETS OF GROUPS, RINGS, AND FIELDS

MINIMAL GENERATING SETS OF GROUPS, RINGS, AND FIELDS MINIMAL GENERATING SETS OF GROUPS, RINGS, AND FIELDS LORENZ HALBEISEN, MARTIN HAMILTON, AND PAVEL RŮŽIČKA Abstract. A subset X of a group (or a ring, or a field) is called generating, if the smallest subgroup

More information

On Riesz-Fischer sequences and lower frame bounds

On Riesz-Fischer sequences and lower frame bounds On Riesz-Fischer sequences and lower frame bounds P. Casazza, O. Christensen, S. Li, A. Lindner Abstract We investigate the consequences of the lower frame condition and the lower Riesz basis condition

More information

Updated Tables of Parameters of (T, M, S)-Nets

Updated Tables of Parameters of (T, M, S)-Nets Updated Tables of Parameters of (T, M, S)-Nets Andrew T. Clayman, 1 K. Mark Lawrence, 2 Gary L. Mullen, 3 Harald Niederreiter, 4 N. J. A. Sloane 5 1 Raytheon Systems Company, State College Operations,

More information

Fermat numbers and integers of the form a k + a l + p α

Fermat numbers and integers of the form a k + a l + p α ACTA ARITHMETICA * (200*) Fermat numbers and integers of the form a k + a l + p α by Yong-Gao Chen (Nanjing), Rui Feng (Nanjing) and Nicolas Templier (Montpellier) 1. Introduction. In 1849, A. de Polignac

More information

A proof of the Jordan normal form theorem

A proof of the Jordan normal form theorem A proof of the Jordan normal form theorem Jordan normal form theorem states that any matrix is similar to a blockdiagonal matrix with Jordan blocks on the diagonal. To prove it, we first reformulate it

More information

CS286.2 Lecture 8: A variant of QPCP for multiplayer entangled games

CS286.2 Lecture 8: A variant of QPCP for multiplayer entangled games CS286.2 Lecture 8: A variant of QPCP for multiplayer entangled games Scribe: Zeyu Guo In the first lecture, we saw three equivalent variants of the classical PCP theorems in terms of CSP, proof checking,

More information

Matrix approximation and Tusnády s problem

Matrix approximation and Tusnády s problem European Journal of Combinatorics 28 (2007) 990 995 www.elsevier.com/locate/ejc Matrix approximation and Tusnády s problem Benjamin Doerr Max-Planck-Institut für Informatik, Stuhlsatzenhausweg 85, 66123

More information

Definition 2.3. We define addition and multiplication of matrices as follows.

Definition 2.3. We define addition and multiplication of matrices as follows. 14 Chapter 2 Matrices In this chapter, we review matrix algebra from Linear Algebra I, consider row and column operations on matrices, and define the rank of a matrix. Along the way prove that the row

More information

AN ELEMENTARY PROOF OF THE SPECTRAL RADIUS FORMULA FOR MATRICES

AN ELEMENTARY PROOF OF THE SPECTRAL RADIUS FORMULA FOR MATRICES AN ELEMENTARY PROOF OF THE SPECTRAL RADIUS FORMULA FOR MATRICES JOEL A. TROPP Abstract. We present an elementary proof that the spectral radius of a matrix A may be obtained using the formula ρ(a) lim

More information

Digital lattice rules

Digital 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 information

QUALITATIVE CONTROLLABILITY AND UNCONTROLLABILITY BY A SINGLE ENTRY

QUALITATIVE CONTROLLABILITY AND UNCONTROLLABILITY BY A SINGLE ENTRY QUALITATIVE CONTROLLABILITY AND UNCONTROLLABILITY BY A SINGLE ENTRY D.D. Olesky 1 Department of Computer Science University of Victoria Victoria, B.C. V8W 3P6 Michael Tsatsomeros Department of Mathematics

More information

Numerical Analysis: Solving Systems of Linear Equations

Numerical Analysis: Solving Systems of Linear Equations Numerical Analysis: Solving Systems of Linear Equations Mirko Navara http://cmpfelkcvutcz/ navara/ Center for Machine Perception, Department of Cybernetics, FEE, CTU Karlovo náměstí, building G, office

More information

Discreteness of Transmission Eigenvalues via Upper Triangular Compact Operators

Discreteness of Transmission Eigenvalues via Upper Triangular Compact Operators Discreteness of Transmission Eigenvalues via Upper Triangular Compact Operators John Sylvester Department of Mathematics University of Washington Seattle, Washington 98195 U.S.A. June 3, 2011 This research

More information

Fibonacci sets and symmetrization in discrepancy theory

Fibonacci sets and symmetrization in discrepancy theory Fibonacci sets and symmetrization in discrepancy theory Dmitriy Bilyk, V.N. Temlyakov, Rui Yu Department of Mathematics, University of South Carolina, Columbia, SC, 29208, U.S.A. Abstract We study the

More information

Linear Dependence of Stationary Distributions in Ergodic Markov Decision Processes

Linear Dependence of Stationary Distributions in Ergodic Markov Decision Processes Linear ependence of Stationary istributions in Ergodic Markov ecision Processes Ronald Ortner epartment Mathematik und Informationstechnologie, Montanuniversität Leoben Abstract In ergodic MPs we consider

More information

USING GRAPHS AND GAMES TO GENERATE CAP SET BOUNDS

USING GRAPHS AND GAMES TO GENERATE CAP SET BOUNDS USING GRAPHS AND GAMES TO GENERATE CAP SET BOUNDS JOSH ABBOTT AND TREVOR MCGUIRE Abstract. Let F 3 be the field with 3 elements and consider the k- dimensional affine space, F k 3, over F 3. A line of

More information

M. VAN BAREL Department of Computing Science, K.U.Leuven, Celestijnenlaan 200A, B-3001 Heverlee, Belgium

M. 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

On the discrepancy of circular sequences of reals

On the discrepancy of circular sequences of reals On the discrepancy of circular sequences of reals Fan Chung Ron Graham Abstract In this paper we study a refined measure of the discrepancy of sequences of real numbers in [0, ] on a circle C of circumference.

More information

The cocycle lattice of binary matroids

The cocycle lattice of binary matroids Published in: Europ. J. Comb. 14 (1993), 241 250. The cocycle lattice of binary matroids László Lovász Eötvös University, Budapest, Hungary, H-1088 Princeton University, Princeton, NJ 08544 Ákos Seress*

More information

Discrete Mathematics. Coquet-type formulas for the rarefied weighted Thue Morse sequence

Discrete Mathematics. Coquet-type formulas for the rarefied weighted Thue Morse sequence Discrete Mathematics 11 (2011) 1724 174 Contents lists available at ScienceDirect Discrete Mathematics journal homepage: www.elsevier.com/locate/disc Coquet-type formulas for the rarefied weighted Thue

More information

Discrete Math, Fourteenth Problem Set (July 18)

Discrete Math, Fourteenth Problem Set (July 18) Discrete Math, Fourteenth Problem Set (July 18) REU 2003 Instructor: László Babai Scribe: Ivona Bezakova 0.1 Repeated Squaring For the primality test we need to compute a X 1 (mod X). There are two problems

More information

Construction of Equidistributed Generators Based on Linear Recurrences Modulo 2

Construction of Equidistributed Generators Based on Linear Recurrences Modulo 2 Construction of Equidistributed Generators Based on Linear Recurrences Modulo 2 Pierre L Ecuyer and François Panneton Département d informatique et de recherche opérationnelle Université de Montréal C.P.

More information

Linear Algebra 2 Spectral Notes

Linear Algebra 2 Spectral Notes Linear Algebra 2 Spectral Notes In what follows, V is an inner product vector space over F, where F = R or C. We will use results seen so far; in particular that every linear operator T L(V ) has a complex

More information

2 K. ENTACHER 2 Generalized Haar function systems In the following we x an arbitrary integer base b 2. For the notations and denitions of generalized

2 K. ENTACHER 2 Generalized Haar function systems In the following we x an arbitrary integer base b 2. For the notations and denitions of generalized BIT 38 :2 (998), 283{292. QUASI-MONTE CARLO METHODS FOR NUMERICAL INTEGRATION OF MULTIVARIATE HAAR SERIES II KARL ENTACHER y Deartment of Mathematics, University of Salzburg, Hellbrunnerstr. 34 A-52 Salzburg,

More information

Fundamentals of Linear Algebra. Marcel B. Finan Arkansas Tech University c All Rights Reserved

Fundamentals of Linear Algebra. Marcel B. Finan Arkansas Tech University c All Rights Reserved Fundamentals of Linear Algebra Marcel B. Finan Arkansas Tech University c All Rights Reserved 2 PREFACE Linear algebra has evolved as a branch of mathematics with wide range of applications to the natural

More information

23 Elements of analytic ODE theory. Bessel s functions

23 Elements of analytic ODE theory. Bessel s functions 23 Elements of analytic ODE theory. Bessel s functions Recall I am changing the variables) that we need to solve the so-called Bessel s equation 23. Elements of analytic ODE theory Let x 2 u + xu + x 2

More information

4. Ergodicity and mixing

4. Ergodicity and mixing 4. Ergodicity and mixing 4. Introduction In the previous lecture we defined what is meant by an invariant measure. In this lecture, we define what is meant by an ergodic measure. The primary motivation

More information

Approximation numbers of Sobolev embeddings - Sharp constants and tractability

Approximation numbers of Sobolev embeddings - Sharp constants and tractability Approximation numbers of Sobolev embeddings - Sharp constants and tractability Thomas Kühn Universität Leipzig, Germany Workshop Uniform Distribution Theory and Applications Oberwolfach, 29 September -

More information

(x 1 +x 2 )(x 1 x 2 )+(x 2 +x 3 )(x 2 x 3 )+(x 3 +x 1 )(x 3 x 1 ).

(x 1 +x 2 )(x 1 x 2 )+(x 2 +x 3 )(x 2 x 3 )+(x 3 +x 1 )(x 3 x 1 ). CMPSCI611: Verifying Polynomial Identities Lecture 13 Here is a problem that has a polynomial-time randomized solution, but so far no poly-time deterministic solution. Let F be any field and let Q(x 1,...,

More information

The 4-periodic spiral determinant

The 4-periodic spiral determinant The 4-periodic spiral determinant Darij Grinberg rough draft, October 3, 2018 Contents 001 Acknowledgments 1 1 The determinant 1 2 The proof 4 *** The purpose of this note is to generalize the determinant

More information

Brownian motion. Samy Tindel. Purdue University. Probability Theory 2 - MA 539

Brownian motion. Samy Tindel. Purdue University. Probability Theory 2 - MA 539 Brownian motion Samy Tindel Purdue University Probability Theory 2 - MA 539 Mostly taken from Brownian Motion and Stochastic Calculus by I. Karatzas and S. Shreve Samy T. Brownian motion Probability Theory

More information

Boundedly complete weak-cauchy basic sequences in Banach spaces with the PCP

Boundedly complete weak-cauchy basic sequences in Banach spaces with the PCP Journal of Functional Analysis 253 (2007) 772 781 www.elsevier.com/locate/jfa Note Boundedly complete weak-cauchy basic sequences in Banach spaces with the PCP Haskell Rosenthal Department of Mathematics,

More information

Additive irreducibles in α-expansions

Additive irreducibles in α-expansions Publ. Math. Debrecen Manuscript Additive irreducibles in α-expansions By Peter J. Grabner and Helmut Prodinger * Abstract. The Bergman number system uses the base α = + 5 2, the digits 0 and, and the condition

More information

QUASI-RANDOM NUMBER GENERATION IN THE VIEW OF PARALLEL IMPLEMENTATION. 1. Introduction

QUASI-RANDOM NUMBER GENERATION IN THE VIEW OF PARALLEL IMPLEMENTATION. 1. Introduction Trends in Mathematics Information Center for Mathematical Sciences Volume 6, Number 2, December, 2003, Pages 155 162 QUASI-RANDOM NUMBER GENERATION IN THE VIEW OF PARALLEL IMPLEMENTATION CHI-OK HWANG Abstract.

More information

On the law of the iterated logarithm for the discrepancy of lacunary sequences

On the law of the iterated logarithm for the discrepancy of lacunary sequences On the law of the iterated logarithm for the discrepancy of lacunary sequences Christoph Aistleitner Abstract A classical result of Philipp (1975) states that for any sequence (n k ) k 1 of integers satisfying

More information

APPLIED 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. 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 information

The Free Central Limit Theorem: A Combinatorial Approach

The Free Central Limit Theorem: A Combinatorial Approach The Free Central Limit Theorem: A Combinatorial Approach by Dennis Stauffer A project submitted to the Department of Mathematical Sciences in conformity with the requirements for Math 4301 (Honour s Seminar)

More information

Linear Algebra and its Applications

Linear Algebra and its Applications Linear Algebra and its Applications 435 (011) 1845 1856 Contents lists available at ScienceDirect Linear Algebra and its Applications journal homepage: wwwelseviercom/locate/laa Hurwitz rational functions

More information

IMPROVEMENT ON THE DISCREPANCY. 1. Introduction

IMPROVEMENT ON THE DISCREPANCY. 1. Introduction t m Mathematical Publications DOI: 10.2478/tmmp-2014-0016 Tatra Mt. Math. Publ. 59 (2014 27 38 IMPROVEMENT ON THE DISCREPANCY OF (t es-sequences Shu Tezuka ABSTRACT. Recently a notion of (t es-sequences

More information

Counting Matrices Over a Finite Field With All Eigenvalues in the Field

Counting Matrices Over a Finite Field With All Eigenvalues in the Field Counting Matrices Over a Finite Field With All Eigenvalues in the Field Lisa Kaylor David Offner Department of Mathematics and Computer Science Westminster College, Pennsylvania, USA kaylorlm@wclive.westminster.edu

More information

Econ Slides from Lecture 7

Econ Slides from Lecture 7 Econ 205 Sobel Econ 205 - Slides from Lecture 7 Joel Sobel August 31, 2010 Linear Algebra: Main Theory A linear combination of a collection of vectors {x 1,..., x k } is a vector of the form k λ ix i for

More information