Tractability of Multivariate Approximation over a Weighted Unanchored Sobolev Space

We study d-variate L2-approximation for a weighted unanchored Sobolev space having smoothness m≥1. This space is equipped with an unusual norm which is, however, equivalent to the norm of the d-fold tensor product of the standard Sobolev space. One might hope that the problem should become easier as its smoothness increases. This is true for our problem if we are only concerned with asymptotic analysis: the nth minimal error is of order n−(m−δ) for any δ>0. However, it is unclear how long we need to wait before this asymptotic behavior kicks in. How does this waiting period depend on d and m? It is easy to prove that no matter how the weights are chosen, the waiting period is at least md, even if the error demand ε is arbitrarily close to 1. Hence, for m≥2, this waiting period is exponential in d, so that the problem suffers from the curse of dimensionality and is intractable. In other words, the fact that the asymptotic behavior improves with m is irrelevant when d is large. So we will be unable to vanquish the curse of dimensionality unless m=1, i.e., unless the smoothness is minimal. In this paper, we prove the more difficult fact that our problem can be tractable if m=1. That is, we can find an ε-approximation using polynomially-many (in d and ε−1) information operations, even if only function values are permitted. When m=1, it is even possible for the problem to be strongly tractable, i.e., we can find an ε-approximation using polynomially-many (in ε−1) information operations, independently of d. These positive results hold when the weights of the Sobolev space decay sufficiently quickly or are bounded finite-order weights, i.e., the d-variate functions we wish to approximate can be decomposed as sums of functions depending on at most ω variables, where ω is independent of d.

[1]  Henryk Wozniakowski,et al.  Information-based complexity , 1987, Nature.

[2]  Henryk Wozniakowski,et al.  Mathematik in den Naturwissenschaften Leipzig Generalized Tractability for Linear Functionals , 2007 .

[3]  R. Caflisch Monte Carlo and quasi-Monte Carlo methods , 1998, Acta Numerica.

[4]  Henryk Wozniakowski,et al.  Intractability Results for Integration and Discrepancy , 2001, J. Complex..

[5]  Henryk Wozniakowski,et al.  Weighted Tensor Product Algorithms for Linear Multivariate Problems , 1999, J. Complex..

[6]  Henryk Wozniakowski,et al.  Tractability of Integration in Non-periodic and Periodic Weighted Tensor Product Hilbert Spaces , 2002, J. Complex..

[7]  Henryk Wozniakowski,et al.  Polynomial-Time Algorithms for Multivariate Linear Problems with Finite-Order Weights: Worst Case Setting , 2005, Found. Comput. Math..

[8]  Henryk Wozniakowski,et al.  Explicit Cost Bounds of Algorithms for Multivariate Tensor Product Problems , 1995, J. Complex..

[9]  Henryk Wozniakowski,et al.  Finite-order weights imply tractability of linear multivariate problems , 2004, J. Approx. Theory.

[10]  Henryk Wozniakowski,et al.  Multivariate L∞ approximation in the worst case setting over reproducing kernel Hilbert spaces , 2008, J. Approx. Theory.

[11]  Henryk Wozniakowski,et al.  On the power of standard information for multivariate approximation in the worst case setting , 2009, J. Approx. Theory.

[12]  H. Woxniakowski Information-Based Complexity , 1988 .

[13]  Harald Niederreiter,et al.  Monte Carlo and Quasi-Monte Carlo Methods 2006 , 2007 .

[14]  Henryk Wozniakowski,et al.  Good Lattice Rules in Weighted Korobov Spaces with General Weights , 2006, Numerische Mathematik.

[15]  Christine Thomas-Agnan,et al.  Computing a family of reproducing kernels for statistical applications , 1996, Numerical Algorithms.

[16]  E. Novak,et al.  Tractability of Multivariate Problems , 2008 .

[17]  G. Wahba,et al.  Design Problems for Optimal Surface Interpolation. , 1979 .

[18]  Zvi Ziegler,et al.  Approximation theory and applications , 1983 .

[19]  Henryk Wozniakowski,et al.  On the Power of Standard Information for Weighted Approximation , 2001, Found. Comput. Math..