A Probabilistic Approach to Floating-Point Arithmetic

arXiv1912.00867MaRDI QIDQ6330262FDOQ6330262

Authors: Fredrik Dahlqvist, Rocco Salvia, George A. Constantinides 

Publication date: 2 December 2019

Abstract: Finite-precision floating point arithmetic unavoidably introduces rounding errors which are traditionally bounded using a worst-case analysis. However, worst-case analysis might be overly conservative because worst-case errors can be extremely rare events in practice. Here we develop a probabilistic model of rounding errors with which it becomes possible to estimate the likelihood that the rounding error of an algorithm lies within a given interval. Given an input distribution, we show how to compute the distribution of rounding errors. We do this exactly for low precision arithmetic, for high precision arithmetic we derive a simple approximation. The model is then entirely compositional: given a numerical program written in a simple imperative programming language we can recursively compute the distribution of rounding errors at each step of the computation and propagate it through each program instruction. This is done by applying a formalism originally developed by Kozen to formalize the semantics of probabilistic programs. We then discuss an implementation of the model and use it to perform probabilistic range analyses on some benchmarks.

This page was built for publication: A Probabilistic Approach to Floating-Point Arithmetic

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q6330262)