The cost of deterministic, adaptive, automatic algorithms: Cones, not balls

• Published in: Journal of Complexity Vol. 30; no. 1; pp. 21 - 45
• Main Authors: Clancy, Nicholas, Ding, Yuhan, Hamilton, Caleb, Hickernell, Fred J., Zhang, Yizhi
• Format: Journal Article
• Language: English
• Published: Elsevier Inc 01.02.2014
• Subjects: Function recovery; Guarantee; Integration; Quadrature; Quadrature; Integration; Guarantee; Function recovery
• ISSN: 0885-064X; 1090-2708
• DOI: 10.1016/j.jco.2013.09.002

Automatic numerical algorithms attempt to provide approximate solutions that differ from exact solutions by no more than a user-specified error tolerance. The computational cost is often determined adaptively by the algorithm based on the function values sampled. While adaptive, automatic algorithms are widely used in practice, most lack guarantees, i.e., conditions on input functions that ensure that the error tolerance is met. This article establishes a framework for guaranteed, adaptive, automatic algorithms. Sufficient conditions for success and two-sided bounds on the computational cost are provided in Theorems 2 and 3. Lower bounds on the complexity of the problem are given in Theorem 6, and conditions under which the proposed algorithms have optimal order are given in Corollary 1. These general theorems are illustrated for univariate numerical integration and function recovery via adaptive algorithms based on linear splines. The key to these adaptive algorithms is performing the analysis for cones of input functions rather than balls. Cones provide a setting where adaption may be beneficial.