2023/09/14 by Pierre Gaillard, Gaillard, Pierre, Sébastien Gerchinovitz +3
Decision Sciences · Engineering · Mathematics · #Advanced Bandit Algorithms Research #Approximation Theory and Sequence Spaces #FOS: Computer and information sciences #FOS: Mathematics #Iterative Methods for Nonlinear Equations #Machine Learning (cs.LG) #Numerical Analysis (math.NA) #Reservoir Engineering and Simulation Methods #Risk and Portfolio Optimization
paper · pdf · doi:10.48550/arxiv.2309.07530
openalex publication_date 2023/09/14 · openalex created_date 2025/10/10 · openalex updated_date 2026/07/28
We study the classical problem of approximating a non-decreasing function f: X → Y in Lp(μ) norm by sequentially querying its values, for known compact real intervals X, Y and a known probability measure μ on \cX. For any function~f we characterize the minimum number of evaluations of f that algorithms need to guarantee an approximation f with an Lp(μ) error below ε after stopping. Unlike worst-case results that hold uniformly over all f, our complexity measure is dependent on each specific function f. To address this problem, we introduce GreedyBox, a generalization of an algorithm originally proposed by Novak (1992) for numerical integration. We prove that GreedyBox achieves an optimal sample complexity for any function f, up to logarithmic factors. Additionally, we uncover results regarding piecewise-smooth functions. Perhaps as expected, the Lp(μ) error of GreedyBox decreases much faster for piecewise-C2 functions than predicted by the algorithm (without any knowledge on the smoothness of f). A simple modification even achieves optimal minimax approximation rates for such functions, which we compute explicitly. In particular, our findings highlight multiple performance gaps between adaptive and non-adaptive algorithms, smooth and piecewise-smooth functions, as well as monotone or non-monotone functions. Finally, we provide numerical experiments to support our theoretical results.