Computing free convolutions via contour integrals

Alice Cortinovis,Lexing Ying
DOI: https://doi.org/10.1142/s2010326324500242
2024-10-12
Random Matrices Theory and Application
Abstract:Random Matrices: Theory and Applications, Ahead of Print. This work proposes algorithms for computing additive and multiplicative free convolutions of two given probability measures. We consider probability measures with compact support whose free convolution results in a probability measure with a density function that exhibits a square-root decay at the boundary (for example, the semicircle distribution or the Marchenko–Pastur distribution). A key ingredient of our method is rewriting the intermediate quantities of the free convolution using the Cauchy integral formula and then discretizing these integrals using the trapezoidal quadrature rule, which converges exponentially fast under suitable analyticity properties of the functions to be integrated.
statistics & probability,physics, mathematical
What problem does this paper attempt to address?