🤖 AI Summary
This paper addresses the stability of discrete approximations for one-dimensional continuous probability distributions under arithmetic operations (addition, subtraction, multiplication, division). We propose a divide-and-conquer quantization algorithm that recursively partitions the support interval and allocates probability mass, enabling efficient discretization of arbitrary continuous distributions with finite mean. Theoretically, we derive the first tight Wasserstein-1 error upper bound applicable to this broad class of distributions and prove that our method achieves optimal convergence rates across multiple settings. Empirically, our approach significantly improves numerical stability and error robustness in distributional arithmetic compared to existing discretization methods. By providing a theoretically grounded and empirically reliable discrete representation, the method establishes a sound foundation for uncertainty propagation, stochastic optimization, and other tasks relying on accurate distributional computations.
📝 Abstract
This article studies a general divide-and-conquer algorithm for approximating continuous one-dimensional probability distributions with finite mean. The article presents a numerical study that compares pre-existing approximation schemes with a special focus on the stability of the discrete approximations when they undergo arithmetic operations. The main results are a simple upper bound of the approximation error in terms of the Wasserstein-1 distance that is valid for all continuous distributions with finite mean. In many use-cases, the studied method achieve optimal rate of convergence, and numerical experiments show that the algorithm is more stable than pre-existing approximation schemes in the context of arithmetic operations.