Lower bounds on the $L_p$ discrepancy of digital NUT sequences
Ralph Kritzinger, Friedrich Pillichshammer

TL;DR
This paper establishes lower bounds on the $L_p$ discrepancy for digital NUT sequences, a significant subclass of digital sequences, contributing to the understanding of their uniformity properties.
Contribution
It provides the first known lower bounds for the $L_p$ discrepancy of specific subclasses of digital NUT sequences, advancing discrepancy theory.
Findings
Lower bounds for $L_p$ discrepancy of digital NUT sequences
Identification of subclasses with provable discrepancy limits
Enhancement of theoretical understanding of digital sequence uniformity
Abstract
We study the discrepancy of digital NUT sequences which are an important sub-class of digital -sequences in the sense of Niederreiter. The main result is a lower bound for certain sub-classes of digital NUT sequences.
Peer Reviews
No public reviews on file for this paper yet. If you reviewed it on a platform where reviews are public (OpenReview, ICLR, NeurIPS, ICML), you can paste yours below so the community can read it here.
Videos
No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.
Lower bounds on the discrepancy of digital NUT sequences
Ralph Kritzinger and Friedrich Pillichshammer The authors are supported by the Austrian Science Fund (FWF): Project F5509-N26, which is a part of the Special Research Program "Quasi-Monte Carlo Methods: Theory and Applications".
Abstract
We study the discrepancy of digital NUT sequences which are an important sub-class of digital -sequences in the sense of Niederreiter. The main result is a lower bound for certain sub-classes of digital NUT sequences.
Keywords: -discrepancy, van der Corput sequence, digital sequence
MSC 2000: 11K38, 11K31
1 Introduction
For a set of points in the (non-normalized) * discrepancy* for is defined as
[TABLE]
(with the usual modification if ), where
[TABLE]
is the (non-normalized) discrepancy function of .
We denote by the set of positive integers and define . Let be an infinite sequence in and, for , let denote the set consisting of the first elements of . It is well known that for all we have
[TABLE]
and
[TABLE]
(For functions , we write or , if there exists a positive constant that is independent of such that or , respectively.) The lower estimate for finite was first shown by Proĭnov [13] (see also [3]) based on famous results of Roth [14] and Schmidt [16] for finite point sets in dimension two. Using the method of Proĭnov in conjunction with a result of Halász [7] for finite point sets in dimension two the lower bound follows also for the -discrepancy. The estimate for was first shown by Schmidt [15] in 1972 (see also [1, 9, 17]).
In this paper we investigate the discrepancy of digital -sequences. Since we only deal with digital sequences over and in dimension 1 we restrict the necessary definitions to this case. For the general setting we refer to [2, 10, 11].
Let be the finite field of order 2, which we identify with the set equipped with arithmetic operations modulo 2. For the generation of a digital sequence over we require an infinite matrix over with the following property111A further technical condition which is sometimes required, see [11, p.72, (S6)], is that for each the sequence becomes eventually zero. Otherwise it could happen that one or more elements of the digital -sequence are 1 and therefore do not belong to .: for every the left upper submatrix has full rank. In order to construct the element for compute the base 2 expansion (which is actually finite), set and compute the matrix vector product
[TABLE]
over . Finally, set
[TABLE]
We denote the digital -sequence222In the general notation, the refers to the dimension and the [math] refers to the full rank condition of the generator matrix . constructed in this way by .
An important sub-class of digital -sequences which is studied in many papers (initiated by Faure [5]) are so-called digital NUT sequences whose generator matrices are of non-singular upper triangular (NUT) form
[TABLE]
For example, if is the identity matrix, then the corresponding digital NUT sequence is the van der Corput sequence in base 2. For information about digital NUT sequences and the van der Corput sequence see the survey [6] and the references therein.
For digital NUT sequences it is known (see [4, Theorem 1]), that
[TABLE]
and for general it is known (see [12, Theorem 2]) that
[TABLE]
Note that according to the lower bound of Schmidt the upper bound for the discrepancy in (2) is optimal in the order of magnitude in . This is not the case for finite .
Concerning lower bounds on the discrepancy of digital NUT sequences very few is known and only for very special cases. For the van der Corput sequence we have for all
[TABLE]
and hence for infinitely many ; see [12, Corollary 1].
For the so-called upper--sequence , which is generated by the matrix
[TABLE]
it is known that for every we have infinitely often; see [4].
In [4] the authors study the discrepancy of for special types of NUT matrices of the form
[TABLE]
with
[TABLE]
Note that these NUT sequences comprise the van der Corput sequence and the upper--sequence as special cases. For let denote the number of rows among the first rows of . For example, in case of the van der Corput sequence and in case of the sequence . Then it follows from [4, Lemma 4] that for every there exists an integer such that . This implies that if we have for every
[TABLE]
In general, however, it is a very difficult task to give precise lower bounds on the discrepancy of digital NUT sequences. We strongly conjecture the following:
Conjecture 1
For every digital NUT sequence we have
[TABLE]
Note that for every digital NUT sequence and for every we have
[TABLE]
where denotes the binary sum-of-digits function which is defined as whenever has binary expansion . The very last estimate in (11) follows from the proof of [8, Theorem 3.5 in Chapter 2].
Remark 1
The result in (11) can be generalized and improved in the following sense: For every and for every digital -sequence we have
[TABLE]
where
[TABLE]
We omit the proof.
The sum-of-digits function is very fluctuating. For example we have , but . In any case we have .
Remark 2
The inequalities in (11) shows that having only very few non-zero binary digits is a sufficient condition on which guarantees that has very low discrepancy. For example we have
[TABLE]
or
[TABLE]
or
[TABLE]
See Figure 1 for a comparison for the van der Corput sequence.
However, the condition on of having very few non-zero binary digits is not a necessary one for low discrepancy. For example, consider of the form . Then we have but: since the discrepancy of and of differ at most by 1 and since we obtain . Hence, while is very large, the discrepancy is low.
But in any case: the only possible candidates of that satisfy (10) are required to have .
In Section 2 we provide a lower bound for for special types of NUT matrices.
2 Lower bound on
We study two sub-classes of NUT matrices. The first class has a certain band structure. More detailled, the considered matrices are of the form where, for fixed ,
[TABLE]
For example, if , we obtain the identity matrix, i.e., .
Theorem 1
For all and we have
[TABLE]
The bound above is satisfied for of the form
[TABLE]
Remark 3
Following all the details in the proof the constant hidden in can be computed exactly. 2. 2.
For we have and hence the resulting NUT sequence is the van der Corput sequence. Theorem 1 gives
[TABLE]
This matches the corresponding value in (3).
We also study NUT matrices which have the same entries in each column above the diagonal; i.e. we deal with matrices of the form
[TABLE]
where is chosen arbitrarily. We set and . For let further be the minimal distance of consecutive zeroes and be the minimal distance of consecutive ones in the string , i.e. for we define
[TABLE]
Theorem 2
Let . For all and for N_{\boldsymbol{a}}=1+\sum_{i=1}^{m-1}2^{i}(1-a_{i})+2^{m}\ \ \ \mbox{for arbitrary m\in\mathbb{N}} we have
[TABLE]
if , and
[TABLE]
if .
Corollary 1
The first elements of a NUT-sequence generated by a matrix of the form satisfy
[TABLE]
for some constant if for some and for all or if for some and for all .
One example for a generator matrix satisfying the hypotheses of Corollary (1) is
[TABLE]
3 The proofs
The following auxiliary result will be the main tool of our proofs.
Lemma 1
For every NUT digital sequence and every of the form with and we have
[TABLE]
where the are given by the following matrix-vector product over :
[TABLE]
where the digits 1 in the latter vector are placed at positions for .
- Proof.
Let be the NUT digital sequence which is generated by the matrix . Let be of the form
[TABLE]
where . For consider
[TABLE]
where for we define . Every
[TABLE]
can be written as
[TABLE]
where and
[TABLE]
For fixed we decompose the matrix in the form
[TABLE]
where is the left upper sub-matrix of . To in (16) we associate
[TABLE]
where are the binary digits of and are the binary digits of . With this notation for in the range (16) we have
[TABLE]
This shows that the point set is a digitally shifted digital net with generating matrix and with digital shift vector
[TABLE]
Since is also a NUT matrix we find that the shift is of the form
[TABLE]
Note that the matrix has full rank, as is a NUT digital sequence. Hence the shifted digital net can be written as the set of points
[TABLE]
where for . Here and in the following denotes addition in .
We emphasize that do not depend on the ’s, whereas the components for may do so. Therefore we can also write
[TABLE]
where for and .
We have the following decomposition of :
[TABLE]
Therefore and from the fact that
[TABLE]
we obtain
[TABLE]
where
[TABLE]
and, for ,
[TABLE]
where for and all other ’s are zero. Note that .
We have
[TABLE]
Observe that
[TABLE]
and
[TABLE]
Hence
[TABLE]
This shows that we have
[TABLE]
From this we obtain for all that
[TABLE]
Hence
[TABLE]
For the very last double sum we have
[TABLE]
Hence
[TABLE]
This gives
[TABLE]
Now we give the proof of Theorem 1.
- Proof.
In order to simplify the notation we will write instead of in the following. For every and we have
[TABLE]
Now choose for , i.e.
[TABLE]
We have
[TABLE]
Therefore
[TABLE]
Hence, using Lemma 1, we get
[TABLE]
Now we have to determine the numbers . Observe that
[TABLE]
where the digits 1 in the latter vector are placed at positions for . From the structure of the matrix we find that
[TABLE]
and therefore
[TABLE]
Furthermore
[TABLE]
Putting all together we obtain
[TABLE]
Hence, using (18), we get
[TABLE]
In the following, we give the proof of Theorem 2.
- Proof.
Note that in the case the numbers appearing in Lemma 1 can also be understood in the following way: Let with , for and . Let
[TABLE]
Then we have for
[TABLE]
Now consider a matrix of the form and set , where . Then we have
[TABLE]
We observe that for and we have if and only if for some , and otherwise. Hence with Lemma 1 we find
[TABLE]
The fact that implies and further
[TABLE]
This completes the proof of the first claim (13). To derive (14) from (13), we show that changing the tuple which defines the matrix to does not change the integral of much. We use the following argument: Let for with
[TABLE]
It is not hard to show that where denotes the distance of a real number to its nearest integer. For as defined above we define the integer and prove . This is the case, since
[TABLE]
and therefore
[TABLE]
This implies inequality (14).
The reference list from the paper itself. Each links out to its DOI / PubMed record.
- 1[1] R. Béjian. Minoration de la discrépance d’une suite quelconque sur T 𝑇 T . Acta Arith. , 41: 185–202, 1982.
- 2[2] J. Dick and F. Pillichshammer. Digital Nets and Sequences: Discrepancy Theory and Quasi-Monte Carlo Integration . Cambridge University Press, Cambridge, 2010.
- 3[3] J. Dick and F. Pillichshammer. Explicit constructions of point sets and sequences with low discrepancy . In: P. Kritzer, H. Niederreiter, F. Pillichshammer, A. Winterhof (eds.), Uniform Distribution and Quasi-Monte Carlo Methods , pages 63–86. Radon Series on Computational and Applied Mathematics 15, De Gruyter, Berlin, 2014.
- 4[4] M. Drmota, G. Larcher and F. Pillichshammer. Precise distribution properties of the van der Corput sequence and related sequences. Manuscripta Math. , 118: 11–41, 2005.
- 5[5] H. Faure. Discrepancy and diaphony of digital ( 0 , 1 ) 0 1 (0,1) -sequences in prime bases. Acta Arith. , 117: 125–148, 2005.
- 6[6] H. Faure, P. Kritzer and F. Pillichshammer. From van der Corput to modern constructions of sequences for quasi-Monte Carlo rules. Indag. Math. , 26: 760–822, 2015.
- 7[7] G. Halász. On Roth’s method in the theory of irregularities of point distributions. In: Recent progress in analytic number theory , Vol. 2, pages 79–94. Academic Press, London-New York, 1981.
- 8[8] L. Kuipers and H. Niederreiter. Uniform distribution of sequences. John Wiley, New York, 1974. Reprint, Dover Publications, Mineola, NY, 2006.
