Universality in Numerical Computations with Random Data. Case Studies
Abstract
The authors present evidence for universality in numerical computations with random data. Given a (possibly stochastic) numerical algorithm with random input data, the time (or number of iterations) to convergence (within a given tolerance) is a random variable, called the halting time. Two-component universality is observed for the fluctuations of the halting time, i.e., the histogram for the halting times, centered by the sample average and scaled by the sample variance, collapses to a universal curve, independent of the input data distribution, as the dimension increases. Thus, up to two components, the sample average and the sample variance, the statistics for the halting time are universally prescribed. The case studies include six standard numerical algorithms, as well as a model of neural computation and decision making. A link to relevant software is provided in for the reader who would like to do computations of his'r own.
Cite
@article{arxiv.1407.3829,
title = {Universality in Numerical Computations with Random Data. Case Studies},
author = {Percy Deift and Govind Menon and Sheehan Olver and Thomas Trogdon},
journal= {arXiv preprint arXiv:1407.3829},
year = {2015}
}