BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Friedrich Pillichshammer

Friedrich Pillichshammer appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
math.NAmath.NTNumerical Analysiseess.ASmath.FAmath.HO

Catalog footprint

What is connected

44works
6topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

44 published item(s)

preprint2023arXiv

Grid-Based Decimation for Wavelet Transforms with Stably Invertible Implementation

The constant center frequency to bandwidth ratio (Q-factor) of wavelet transforms provides a very natural representation for audio data. However, invertible wavelet transforms have either required non-uniform decimation -- leading to irregular data structures that are cumbersome to work with -- or require excessively high oversampling with unacceptable computational overhead. Here, we present a novel decimation strategy for wavelet transforms that leads to stable representations with oversampling rates close to one and uniform decimation. Specifically, we show that finite implementations of the resulting representation are energy-preserving in the sense of frame theory. The obtained wavelet coefficients can be stored in a timefrequency matrix with a natural interpretation of columns as time frames and rows as frequency channels. This matrix structure immediately grants access to a large number of algorithms that are successfully used in time-frequency audio processing, but could not previously be used jointly with wavelet transforms. We demonstrate the application of our method in processing based on nonnegative matrix factorization, in onset detection, and in phaseless reconstruction.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

A note on isotropic discrepancy and spectral test of lattice point sets

We show that the isotropic discrepancy of a lattice point set can be bounded from below and from above in terms of the spectral test of the corresponding integration lattice. From this we deduce that the isotropic discrepancy of any $N$-element lattice point set in $[0,1)^d$ is at least of order $N^{-1/d}$. This order of magnitude is best possible for lattice point sets in dimension $d$.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Tractability of approximation in the weighted Korobov space in the worst-case setting

In this paper we consider $L_p$-approximation, $p \in \{2,\infty\}$, of periodic functions from weighted Korobov spaces. In particular, we discuss tractability properties of such problems, which means that we aim to relate the dependence of the information complexity on the error demand $\varepsilon$ and the dimension $d$ to the decay rate of the weight sequence $(γ_j)_{j \ge 1}$ assigned to the Korobov space. Some results have been well known since the beginning of this millennium, others have been proven quite recently. We give a survey of these findings and will add some new results on the $L_\infty$-approximation problem. To conclude, we give a concise overview of results and collect a number of interesting open problems.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A note on the periodic $L_2$-discrepancy of Korobov's $p$-sets

We study the periodic $L_2$-discrepancy of point sets in the $d$-dimensional torus. This discrepancy is intimately connected with the root-mean-square $L_2$-discrepancy of shifted point sets, with the notion of diaphony, and with the worst case error of cubature formulas for the integration of periodic functions in Sobolev spaces of mixed smoothness. In discrepancy theory many results are based on averaging arguments. In order to make such results relevant for applications one requires explicit constructions of point sets with ``average'' discrepancy. In our main result we study Korobov's $p$-sets and show that this point sets have periodic $L_2$-discrepancy of average order. This result is related to an open question of Novak and Woźniakowski.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Exponential tractability of linear weighted tensor product problems in the worst-case setting for arbitrary linear functionals

We study the approximation of compact linear operators defined over certain weighted tensor product Hilbert spaces. The information complexity is defined as the minimal number of arbitrary linear functionals which is needed to obtain an $\varepsilon$-approximation for the $d$-variate problem. It is fully determined in terms of the weights and univariate singular values. Exponential tractability means that the information complexity is bounded by a certain function which depends polynomially on $d$ and logarithmically on $\varepsilon^{-1}$. The corresponding un-weighted problem was studied recently by Hickernell, Kritzer and Woźniakowski with many negative results for exponential tractability. The product weights studied in the present paper change the situation. Depending on the form of polynomial dependence on $d$ and logarithmic dependence on $\varepsilon^{-1}$, we study exponential strong polynomial, exponential polynomial, exponential quasi-polynomial, and exponential $(s,t)$-weak tractability with $\max(s,t)\ge1$. For all these notions of exponential tractability, we establish necessary and sufficient conditions on weights and univariate singular values for which it is indeed possible to achieve the corresponding notion of exponential tractability. The case of exponential $(s,t)$-weak tractability with $\max(s,t)<1$ is left for future study.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

On Quasi-Monte Carlo Methods in Weighted ANOVA Spaces

In the present paper we study quasi-Monte Carlo rules for approximating integrals over the $d$-dimensional unit cube for functions from weighted Sobolev spaces of regularity one. While the properties of these rules are well understood for anchored Sobolev spaces, this is not the case for the ANOVA spaces, which are another very important type of reference spaces for quasi-Monte Carlo rules. Using a direct approach we provide a formula for the worst case error of quasi-Monte Carlo rules for functions from weighted ANOVA spaces. As a consequence we bound the worst case error from above in terms of weighted discrepancy of the employed integration nodes. On the other hand we also obtain a general lower bound in terms of the number $n$ of used integration nodes. For the one-dimensional case our results lead to the optimal integration rule and also in the two-dimensional case we provide rules yielding optimal convergence rates.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Weighted integration over a cube based on digital nets and sequences

Quasi-Monte Carlo (QMC) methods are equal weight quadrature rules to approximate integrals over the unit cube with respect to the uniform measure. In this paper we discuss QMC integration with respect to general product measures defined on an arbitrary cube. We only require that the cumulative distribution function is invertible. We develop a worst-case error bound and study the dependence of the error on the number of points and the dimension for digital nets and sequences as well as polynomial lattice point sets, which are mapped to the domain using the inverse cumulative distribution function. We do not require any smoothness properties of the probability density function and the worst-case error does not depend on the particular choice of density function and its smoothness. The component-by-component construction of polynomial lattice rules is based on a criterion which depends only on the size of the cube but is otherwise independent of the product measure.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

$\boldsymbol{L}_{\infty}$-approximation in Korobov spaces with Exponential Weights

We study multivariate $\boldsymbol{L}_{\infty}$-approximation for a weighted Korobov space of periodic functions for which the Fourier coefficients decay exponentially fast. The weights are defined, in particular, in terms of two sequences $\boldsymbol{a}=\{a_j\}$ and $\boldsymbol{b}=\{b_j\}$ of positive real numbers bounded away from zero. We study the minimal worst-case error $e^{\boldsymbol{L}_{\infty}\mathrm{-app},Λ}(n,s)$ of all algorithms that use $n$ information evaluations from a class $Λ$ in the $s$-variate case. We consider two classes $Λ$ in this paper: the class $Λ^{\rm all}$ of all linear functionals and the class $Λ^{\rm std}$ of only function evaluations. We study exponential convergence of the minimal worst-case error, which means that $e^{\boldsymbol{L}_{\infty}\mathrm{-app},Λ}(n,s)$ converges to zero exponentially fast with increasing $n$. Furthermore, we consider how the error depends on the dimension $s$. To this end, we define the notions of $κ$-EC-weak, EC-polynomial and EC-strong polynomial tractability, where EC stands for "exponential convergence". In particular, EC-polynomial tractability means that we need a polynomial number of information evaluations in $s$ and $1+\log\,\varepsilon^{-1}$ to compute an $\varepsilon$-approximation. We derive necessary and sufficient conditions on the sequences $\boldsymbol{a}$ and $\boldsymbol{b}$ for obtaining exponential error convergence, and also for obtaining the various notions of tractability. The results are the same for both classes $Λ$.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Discrepancy of second order digital sequences in function spaces with dominating mixed smoothness

The discrepancy function measures the deviation of the empirical distribution of a point set in $[0,1]^d$ from the uniform distribution. In this paper, we study the classical discrepancy function with respect to the BMO and exponential Orlicz norms, as well as Sobolev, Besov and Triebel-Lizorkin norms with dominating mixed smoothness. We give sharp bounds for the discrepancy function under such norms with respect to infinite sequences.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

On Weyl products and uniform distribution modulo one

In the present paper we study the asymptotic behavior of trigonometric products of the form $\prod_{k=1}^N 2 \sin(πx_k)$ for $N \to \infty$, where the numbers $ω=(x_k)_{k=1}^N$ are evenly distributed in the unit interval $[0,1]$. The main result are matching lower and upper bounds for such products in terms of the star-discrepancy of the underlying points $ω$, thereby improving earlier results obtained by Hlawka in 1969. Furthermore, we consider the special cases when the points $ω$ are the initial segment of a Kronecker or van der Corput sequence. The paper concludes with some probabilistic analogues.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Truncation Dimension for Function Approximation

We consider approximation of functions of $s$ variables, where $s$ is very large or infinite, that belong to weighted anchored spaces. We study when such functions can be approximated by algorithms designed for functions with only very small number ${\rm dim^{trnc}}(\varepsilon)$ of variables. Here $\varepsilon$ is the error demand and we refer to ${\rm dim^{trnc}}(\varepsilon)$ as the $\varepsilon$-truncation dimension. We show that for sufficiently fast decaying product weights and modest error demand (up to about $\varepsilon \approx 10^{-5}$) the truncation dimension is surprisingly very small.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

$L_p$-discrepancy of the symmetrized van der Corput sequence

It is well known that the $L_p$-discrepancy for $p \in [1,\infty]$ of the van der Corput sequence is of exact order of magnitude $O((\log N)/N)$. This however is for $p \in (1,\infty)$ not best possible with respect to the lower bounds according to Roth and Proinov. For the case $p=2$ it is well known that the symmetrization trick due to Davenport leads to the optimal $L_2$-discrepancy rate $O(\sqrt{\log N}/N)$ for the symmetrized van der Corput sequence. In this note we show that this result holds for all $p \in (1,\infty)$. The proof is based on an estimate of the Haar coefficients of the corresponding local discrepancy and on the use of the Littlewood-Paley inequality.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Approximation in Hermite spaces of smooth functions

We consider $\mathbb{L}_2$-approximation of elements of a Hermite space of analytic functions over $\mathbb{R}^s$. The Hermite space is a weighted reproducing kernel Hilbert space of real valued functions for which the Hermite coefficients decay exponentially fast. The weights are defined in terms of two sequences $\boldsymbol{a} = \{a_j\}$ and $\boldsymbol{b} = \{b_j\}$ of positive real numbers. We study the $n$th minimal worst-case error $e(n,{\rm APP}_s;Λ^{\rm std})$ of all algorithms that use $n$ information evaluations from the class $Λ^{\rm std}$ which only allows function evaluations to be used. We study (uniform) exponential convergence of the $n$th minimal worst-case error, which means that $e(n,{\rm APP}_s; Λ^{\rm std})$ converges to zero exponentially fast with increasing $n$. Furthermore, we consider how the error depends on the dimension $s$. To this end, we study the minimal number of information evaluations needed to compute an $\varepsilon$-approximation by considering several notions of tractability which are defined with respect to $s$ and $\log \varepsilon^{-1}$. We derive necessary and sufficient conditions on the sequences $\boldsymbol{a}$ and $\boldsymbol{b}$ for obtaining exponential error convergence, and also for obtaining the various notions of tractability. It turns out that the conditions on the weight sequences are almost the same as for the information class $Λ^{\rm all}$ which uses all linear functionals. The results are also constructive as the considered algorithms are based on tensor products of Gauss-Hermite rules for multivariate integration. The obtained results are compared with the analogous results for integration in the same Hermite space. This allows us to give a new sufficient condition for EC-weak tractability for integration.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Component-by-component construction of shifted Halton sequences

We study quasi-Monte Carlo integration in a weighted anchored Sobolev space. As the underlying integration nodes we consider Halton sequences in prime bases $\boldsymbol{p}=(p_1,\ldots,p_s)$ which are shifted with a $\boldsymbol{p}$-adic shift based on $\boldsymbol{p}$-adic arithmetic. The error is studied in the worst-case setting. In a recent paper, Hellekalek together with the authors of this article proved optimal error bounds in the root mean square sense, where the mean was extended over the uncountable set of all possible $\boldsymbol{p}$-adic shifts. Here we show that candidates for good shifts can in fact be chosen from a finite set and can be found by a component-by-component algorithm.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Construction algorithms for plane nets in base $b$

The class of $(0,m,s)$-nets in base $b$ has been introduced by Niederreiter as examples of point sets in the $s$-dimensional unit cube with excellent uniform distribution properties. In particular such nets have been proved to have very low discrepancy. This property is essential for the use of nets in quasi-Monte Carlo rules for numerical integration. In this short note we propose two algorithms for the construction of plane $(0,m,2)$-nets in base~$b$.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Digital inversive vectors can achieve strong polynomial tractability for the weighted star discrepancy and for multivariate integration

We study high-dimensional numerical integration in the worst-case setting. The subject of tractability is concerned with the dependence of the worst-case integration error on the dimension. Roughly speaking, an integration problem is tractable if the worst-case error does not explode exponentially with the dimension. Many classical problems are known to be intractable. However, sometimes tractability can be shown. Often such proofs are based on randomly selected integration nodes. Of course, in applications true random numbers are not available and hence one mimics them with pseudorandom number generators. This motivates us to propose the use of pseudorandom vectors as underlying integration nodes in order to achieve tractability. In particular, we consider digital inverse vectors and present two examples of problems, the weighted star discrepancy and integration of Hölder continuous, absolute convergent Fourier- and cosine series, where the proposed method is successful.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

From van der Corput to modern constructions of sequences for quasi-Monte Carlo rules

In 1935 J.G. van der Corput introduced a sequence which has excellent uniform distribution properties modulo 1. This sequence is based on a very simple digital construction scheme with respect to the binary digit expansion. Nowadays the van der Corput sequence, as it was named later, is the prototype of many uniformly distributed sequences, also in the multi-dimensional case. Such sequences are required as sample nodes in quasi-Monte Carlo algorithms, which are deterministic variants of Monte Carlo rules for numerical integration. Since its introduction many people have studied the van der Corput sequence and generalizations thereof. This led to a huge number of results. On the occasion of the 125th birthday of J.G. van der Corput we survey many interesting results on van der Corput sequences and their generalizations. In this way we move from van der Corput's ideas to the most modern constructions of sequences for quasi-Monte Carlo rules, such as, e.g., generalized Halton sequences or Niederreiter's $(t,s)$-sequences.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Integration and approximation in cosine spaces of smooth functions

We study multivariate integration and approximation for functions belonging to a weighted reproducing kernel Hilbert space based on half-period cosine functions in the worst-case setting. The weights in the norm of the function space depend on two sequences of real numbers and decay exponentially. As a consequence the functions are infinitely often differentiable, and therefore it is natural to expect exponential convergence of the worst-case error. We give conditions on the weight sequences under which we have exponential convergence for the integration as well as the approximation problem. Furthermore, we investigate the dependence of the errors on the dimension by considering various notions of tractability. We prove sufficient and necessary conditions to achieve these tractability notions.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

On Equivalence of Anchored and ANOVA Spaces; Lower Bounds

We provide lower bounds for the norms of embeddings between $\boldsymbolγ$-weighted Anchored and ANOVA spaces of $s$-variate functions with mixed partial derivatives of order one bounded in $L_p$ norm ($p\in[1,\infty]$). In particular we show that the norms behave polynomially in $s$ for Finite Order Weights and Finite Diameter Weights if $p>1$, and increase faster than any polynomial in $s$ for Product Order-Dependent Weights and any $p$.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Tractability of Multivariate Approximation Defined over Hilbert Spaces with Exponential Weights

We study multivariate approximation defined over tensor product Hilbert spaces. The domain space is a weighted tensor product Hilbert space with exponential weights which depend on two sequences $\boldsymbol{a}=\{a_j\}_{j\in\mathbb{N}}$ and $\boldsymbol{b}=\{b_j\}_{j\in\mathbb{N}}$ of positive numbers, and on a bounded sequence of positive integers $\boldsymbol{m}=\{m_j\}_{j\in\mathbb{N}}$. The sequence $\boldsymbol{a}$ is non-decreasing and the sequence $\boldsymbol{b}$ is bounded from below by a positive number. We find necessary and sufficient conditions on $\boldsymbol{a},\boldsymbol{b}$ and $\boldsymbol{m}$ to achieve the standard and new notions of tractability in the worst case setting.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A reduced fast component-by-component construction of lattice points for integration in weighted spaces with fast decreasing weights

Lattice rules and polynomial lattice rules are quadrature rules for approximating integrals over the $s$-dimensional unit cube. Since no explicit constructions of such quadrature methods are known for dimensions $s > 2$, one usually has to resort to computer search algorithms. The fast component-by-component approach is a useful algorithm for finding suitable quadrature rules. We present a modification of the fast component-by-component algorithm which yields savings of the construction cost for (polynomial) lattice rules in weighted function spaces. The idea is to reduce the size of the search space for coordinates which are associated with small weights and are therefore of less importance to the overall error compared to coordinates associated with large weights. We analyze tractability conditions of the resulting QMC rules. Numerical results demonstrate the effectiveness of our method.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Discrepancy estimates for index-transformed uniformly distributed sequences

In this paper we show discrepancy bounds for index-transformed uniformly distributed sequences. From a general result we deduce very tight lower and upper bounds on the discrepancy of index-transformed van der Corput-, Halton-, and $(t,s)$-sequences indexed by the sum-of-digits function. We also analyze the discrepancy of sequences indexed by other functions, such as, e.g., $\lfloor n^α\rfloor$ with $0 < α< 1$.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Integration in Hermite spaces of analytic functions

We study integration in a class of Hilbert spaces of analytic functions defined on the $\mathbb{R}^s$. The functions are characterized by the property that their Hermite coefficients decay exponentially fast. We use Gauss-Hermite integration rules and show that the error of our algorithms decays exponentially fast. Furthermore, we give necessary and sufficient conditions under which we achieve exponential convergence with weak, polynomial, and strong polynomial tractability.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Numerical integration in $\log$-Korobov and $\log$-cosine spaces

QMC rules are equal weight quadrature rules for approximating integrals over $[0,1]^s$. One line of research studies the integration error of functions in the unit ball of so-called Korobov spaces, which are Hilbert spaces of periodic functions on $[0,1]^s$ with square integrable partial mixed derivatives of order $α$. Using Parseval's identity, this smoothness can be defined for all real numbers $α> 1/2$. This condition is necessary as otherwise the Korobov space contains discontinuous functions for which function evaluation is not well defined. This paper is concerned with more precise endpoint estimates of the integration error using QMC rules for Korobov spaces with $α$ arbitrarily close to $1/2$. To obtain such estimates we introduce a $\log$-scale for functions with smoothness close to $1/2$, which we call $\log$-Korobov spaces. We show that lattice rules can be used to obtain an integration error of order $\mathcal{O}(N^{-1/2} (\log N)^{-μ(1-λ)/2})$ for any $1/μ<λ\le 1$, where $μ>1$ is a power in the $\log$-scale. We also consider tractability of numerical integration for weighted Korobov spaces with product weights $(γ_j)_{j \in \mathbb{N}}$. It is known that if $\sum_{j=1}^\infty γ_j^τ< \infty$ for some $1/(2α) < τ\le 1$ one can obtain error bounds which are independent of the dimension. In this paper we give a more refined estimate for the case where $τ$ is close to $1/(2 α)$, namely we show dimension independent error bounds under the condition that $\sum_{j=1}^\infty γ_j \max\{1, \log γ_j^{-1}\}^{μ(1-λ)} < \infty$ for some $1/μ< λ\le 1$. The essential tool in our analysis is a $\log$-scale Jensen's inequality. The results described above also apply to integration in $\log$-cosine spaces using tent-transformed lattice rules.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Open type quasi-Monte Carlo integration based on Halton sequences in weighted Sobolev spaces

In this paper, we study quasi-Monte Carlo (QMC) integration in weighted Sobolev spaces. In contrast to many previous results the QMC algorithms considered here are of open type, i.e., they are extensible in the number of sample points without having to discard the samples already used. As the underlying integration nodes we consider randomized Halton sequences in prime bases $\boldsymbol{p}=(p_1,...,p_s)$ for which we study the root mean square (RMS) worst-case error. The randomization method is a $\boldsymbol{p}$-adic shift which is based on $\boldsymbol{p}$-adic arithmetic. The obtained error bounds are optimal in the order of magnitude of the number of sample nodes. Furthermore we obtain conditions on the coordinate weights under which the error bounds are independent of the dimension $s$. In terms of the field of Information-Based Complexity this means that the corresponding QMC rule achieves a strong polynomial tractability error bound. Our findings on the RMS worst-case error of randomized Halton sequences can be carried over to the RMS $L_2$-discrepancy. Except for the $\boldsymbol{p}$-adic shift our results are fully constructive and no search algorithms (such as the component-by-component algorithm) are required.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Optimal order of $L_p$-discrepancy of digit shifted Hammersley point sets in dimension 2

It is well known that the two-dimensional Hammersley point set consisting of $N=2^n$ elements (also known as Roth net) does not have optimal order of $L_p$-discrepancy for $p \in (1,\infty)$ in the sense of the lower bounds according to Roth (for $p \in [2,\infty)$) and Schmidt (for $p \in (1,2)$). On the other hand, it is also known that slight modifications of the Hammersley point set can lead to the optimal order $\sqrt{\log N}/N$ of $L_2$-discrepancy, where $N$ is the number of points. Among these are for example digit shifts or the symmetrization. In this paper we show that these modified Hammersley point sets also achieve optimal order of $L_p$-discrepancy for all $p \in (1,\infty)$.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Proof Techniques in Quasi-Monte Carlo Theory

In this survey paper we discuss some tools and methods which are of use in quasi-Monte Carlo (QMC) theory. We group them in chapters on Numerical Analysis, Harmonic Analysis, Algebra and Number Theory, and Probability Theory. We do not provide a comprehensive survey of all tools, but focus on a few of them, including reproducing and covariance kernels, Littlewood-Paley theory, Riesz products, Minkowski's fundamental theorem, exponential sums, diophantine approximation, Hoeffding's inequality and empirical processes, as well as other tools. We illustrate the use of these methods in QMC using examples.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The $\boldsymbol{p}$-adic diaphony of the Halton sequence

The $\boldsymbol{p}$-adic diaphony as introduced by Hellekalek is a quantitative measure for the irregularity of distribution of a sequence in the unit cube. In this paper we show how this notion of diaphony can be interpreted as worst-case integration error in a certain reproducing kernel Hilbert space. Our main result is an upper bound on the $\boldsymbol{p}$-adic diaphony of the Halton sequence.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The inverse of the star-discrepancy problem and the generation of pseudo-random numbers

The inverse of the star-discrepancy problem asks for point sets $P_{N,s}$ of size $N$ in the $s$-dimensional unit cube $[0,1]^s$ whose star-discrepancy $D^\ast(P_{N,s})$ satisfies $$D^\ast(P_{N,s}) \le C \sqrt{s/N},$$ where $C> 0$ is a constant independent of $N$ and $s$. The first existence results in this direction were shown by Heinrich, Novak, Wasilkowski, and Woźniakowski in 2001, and a number of improvements have been shown since then. Until now only proofs that such point sets exist are known. Since such point sets would be useful in applications, the big open problem is to find explicit constructions of suitable point sets $P_{N,s}$. We review the current state of the art on this problem and point out some connections to pseudo-random number generators.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The weighted star discrepancy of Korobov's $p$-sets

We analyze the weighted star discrepancy of so-called $p$-sets which go back to definitions due to Korobov in the 1950s and Hua and Wang in the 1970s. Since then, these sets have largely been ignored since a number of other constructions have been discovered which achieve a better convergence rate. However, it has recently been discovered that the $p$-sets perform well in terms of the dependence on the dimension. We prove bounds on the weighted star discrepancy of the $p$-sets which hold for any choice of weights. For product weights we give conditions under which the discrepancy bounds are independent of the dimension $s$. This implies strong polynomial tractability for the weighted star discrepancy. We also show that a very weak condition on the product weights suffices to achieve polynomial tractability.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Tractability of multivariate analytic problems

In the theory of tractability of multivariate problems one usually studies problems with finite smoothness. Then we want to know which $s$-variate problems can be approximated to within $\varepsilon$ by using, say, polynomially many in $s$ and $\varepsilon^{-1}$ function values or arbitrary linear functionals. There is a recent stream of work for multivariate analytic problems for which we want to answer the usual tractability questions with $\varepsilon^{-1}$ replaced by $1+\log \varepsilon^{-1}$. In this vein of research, multivariate integration and approximation have been studied over Korobov spaces with exponentially fast decaying Fourier coefficients. This is work of J. Dick, G. Larcher, and the authors. There is a natural need to analyze more general analytic problems defined over more general spaces and obtain tractability results in terms of $s$ and $1+\log \varepsilon^{-1}$. The goal of this paper is to survey the existing results, present some new results, and propose further questions for the study of tractability of multivariate analytic questions.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A metrical lower bound on the star discrepancy of digital sequences

In this paper we study uniform distribution properties of digital sequences over a finite field of prime order. In 1998 it was shown by Larcher that for almost all $s$-dimensional digital sequences the star discrepancy $D_N^\ast$ satisfies an upper bound of the form $D_N^\ast=O((\log N)^s (\log \log N)^{2+\varepsilon})$ for any $\varepsilon>0$. Generally speaking it is much more difficult to obtain good lower bounds for specific sequences than upper bounds. Here we show that Larchers result is best possible up to some $\log \log N$ term. More detailed, we prove that for almost all $s$-dimensional digital sequences the star discrepancy satisfies $D_N^\ast \ge c(q,s) (\log N)^s \log \log N$ for infinitely many $N \in \NN$, where $c(q,s)>0$ only depends on $q$ and $s$ but not on $N$.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Explicit constructions of point sets and sequences with low discrepancy

In this article we survey recent results on the explicit construction of finite point sets and infinite sequences with optimal order of $\mathcal{L}_q$ discrepancy. In 1954 Roth proved a lower bound for the $\mathcal{L}_2$ discrepancy of finite point sets in the unit cube of arbitrary dimension. Later various authors extended Roth's result to lower bounds also for the $\mathcal{L}_q$ discrepancy and for infinite sequences. While it was known already from the early 1980s on that Roth's lower bound is best possible in the order of magnitude, it was a longstanding open question to find explicit constructions of point sets and sequences with optimal order of $\mathcal{L}_2$ discrepancy. This problem was solved by Chen and Skriganov in 2002 for finite point sets and recently by the authors of this article for infinite sequences. These constructions can also be extended to give optimal order of the $\mathcal{L}_q$ discrepancy of finite point sets for $q \in (1,\infty)$. The main aim of this article is to give an overview of these constructions and related results.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Metrical lower bounds on the discrepancy of digital Kronecker-sequences

Digital Kronecker-sequences are a non-archimedean analog of classical Kronecker-sequences whose construction is based on Laurent series over a finite field. In this paper it is shown that for almost all digital Kronecker-sequences the star discrepancy satisfies $D_N^\ast \ge c(q,s) (\log N)^s \log \log N$ for infinitely many $N \in \NN$, where $c(q,s)>0$ only depends on the dimension $s$ and on the order $q$ of the underlying finite field, but not on $N$. This result shows that a corresponding metrical upper bound due to Larcher is up to some $\log \log N$ term best possible.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Optimal $\mathcal{L}_2$ discrepancy bounds for higher order digital sequences over the finite field $\mathbb{F}_2$

We show that the $\mathcal{L}_2$ discrepancy of the explicitly constructed infinite sequences of points $(\boldsymbol{x}_0,\boldsymbol{x}_1, \boldsymbol{x}_2,...)$ in $[0,1)^s$ over $\mathbb{F}_2$ introduced in [J. Dick, Walsh spaces containing smooth functions and quasi-Monte Carlo rules of arbitrary high order. SIAM J. Numer. Anal., {\bf 46}, 1519--1553, 2008] satisfy $$\mathcal{L}_{2,N}(\{\boldsymbol{x}_0,\boldsymbol{x}_1,..., \boldsymbol{x}_{N-1}\}) \le C_s N^{-1} (\log N)^{s/2} \quad {for all} N \ge 2,$$ and $$\mathcal{L}_{2,2^m}(\{\boldsymbol{x}_0,\boldsymbol{x}_1,..., \boldsymbol{x}_{2^m-1}\}) \le C_s 2^{-m} m^{(s-1)/2} \quad {for all} m \ge 1,$$ where $C_s > 0$ is a constant independent of $N$ and $m$. These results are best possible by lower bounds in [P.D. Proinov, On the $L^2$ discrepancy of some infinite sequences. Serdica, {\bf 11}, 3--12, 1985] and [K. F. Roth, On irregularities of distribution. Mathematika, {\bf 1}, 73--79, 1954]. Further, for every $N \ge 2$ we explicitly construct finite point sets $\{\boldsymbol{y}_0,..., \boldsymbol{y}_{N-1}\}$ in $[0,1)^s$ such that $$\mathcal{L}_{2,N}(\{\boldsymbol{y}_0,\boldsymbol{y}_1,..., \boldsymbol{y}_{N-1}\}) \le C_s N^{-1} (\log N)^{(s-1)/2}.$$ Another solution for finite point sets by a different construction was previously shown in [W. W. L. Chen and M. M. Skriganov, Explicit constructions in the classical mean squares problem in irregularity of point distribution. J. Reine Angew. Math., {\bf 545}, 67--95, 2002].

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Approximation of analytic functions in Korobov spaces

We study multivariate $L_2$-approximation for a weighted Korobov space of analytic periodic functions for which the Fourier coefficients decay exponentially fast. The weights are defined, in particular, in terms of two sequences $\boldsymbol{a} =\{a_j\}$ and $\boldsymbol{b} =\{b_j\}$ of numbers no less than one. Let $e^{L_2-\mathrm{app},Λ}(n,s)$ be the minimal worst-case error of all algorithms that use $n$ information functionals from the class $Λ$ in the $s$-variate case. We consider two classes $Λ$: the class $Λ^{\rm all}$ consists of all linear functionals and the class $Λ^{\rm std}$ consists of only function valuations. We study (EXP) exponential convergence. This means that $$ e^{L_2-\mathrm{app},Λ}(n,s) \le C(s)\,q^{\,(n/C_1(s))^{p(s)}}\quad{for all}\quad n, s \in \mathbb{N} $$ where $q\in(0,1)$, and $C,C_1,p:\mathbb{N} \rightarrow (0,\infty)$. If we can take $p(s)=p>0$ for all $s$ then we speak of (UEXP) uniform exponential convergence. We also study EXP and UEXP with (WT) weak, (PT) polynomial and (SPT) strong polynomial tractability. These concepts are defined as follows. Let $n(\e,s)$ be the minimal $n$ for which $e^{L_2-\mathrm{app},Λ}(n,s)\le \e$. Then WT holds iff $\lim_{s+\log\,\e^{-1}\to\infty}(\log n(\e,s))/(s+\log\,\e^{-1})=0$, PT holds iff there are $c,τ_1,τ_2$ such that $n(\e,s)\le cs^{τ_1}(1+\log\,\e^{-1})^{τ_2}$ for all $s$ and $\e\in(0,1)$, and finally SPT holds iff the last estimate holds for $τ_1=0$. The infimum of $τ_2$ for which SPT holds is called the exponent $τ^*$ of SPT. We prove that the results are the same for both classes $Λ$, and obtain conditions for WT, PT, SPT with and without EXP and UEXP.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Lattice rules for nonperiodic smooth integrands

The aim of this paper is to show that one can achieve convergence rates of $N^{-α+ δ}$ for $α> 1/2$ (and for $δ> 0$ arbitrarily small) for nonperiodic $α$-smooth cosine series using lattice rules without random shifting. The smoothness of the functions can be measured by the decay rate of the cosine coefficients. For a specific choice of the parameters the cosine series space coincides with the unanchored Sobolev space of smoothness 1. We study the embeddings of various reproducing kernel Hilbert spaces and numerical integration in the cosine series function space and show that by applying the so-called tent transformation to a lattice rule one can achieve the (almost) optimal rate of convergence of the integration error. The same holds true for symmetrized lattice rules for the tensor product of the direct sum of the Korobov space and cosine series space, but with a stronger dependence on the dimension in this case.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

On the existence of hyperplane sequences, with quality parameter and discrepancy bounds

It is well-known that digital $(t,m,s)$-nets and $(\Tfett,s)$-sequences over a finite field have excellent properties when they are used as underlying nodes in quasi-Monte Carlo integration rules. One very general sub-class of digital nets are hyperplane nets which can be viewed as a generalization of cyclic nets and of polynomial lattice point sets. In this paper we introduce infinite versions of hyperplane nets and call these sequences hyperplane sequences. Our construction is based on the recent duality theory for digital sequences according to Dick and Niederreiter. We then analyze the equidistribution properties of hyperplane sequences in terms of the quality function $\Tfett$ and the star discrepancy.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Efficient calculation of the worst-case error and (fast) component-by-component construction of higher order polynomial lattice rules

We show how to obtain a fast component-by-component construction algorithm for higher order polynomial lattice rules. Such rules are useful for multivariate quadrature of high-dimensional smooth functions over the unit cube as they achieve the near optimal order of convergence. The main problem addressed in this paper is to find an efficient way of computing the worst-case error. A general algorithm is presented and explicit expressions for base~2 are given. To obtain an efficient component-by-component construction algorithm we exploit the structure of the underlying cyclic group. We compare our new higher order multivariate quadrature rules to existing quadrature rules based on higher order digital nets by computing their worst-case error. These numerical results show that the higher order polynomial lattice rules improve upon the known constructions of quasi-Monte Carlo rules based on higher order digital nets.

0 citations0 reviewsNo code linkedOpen access