Simulations of Levy flights

Levy flights, also known as alpha-stable Levy processes or heavy-tailed statistics, are becoming a commonly used tool in optics. Nonetheless, the different parametrizations and the absence of any analytic expression for the distribution functions (apart from some exceptions) makes it difficult to efficiently simulate such processes. We review and compare three algorithms for the generation of sequences of symmetric stable Levy random variables.