Generation of Gaussian distributed random numbers by using a numerical inversion method

Abstract We describe a vectorizable implementation of a numerical inversion method to generate approximately Gaussian distributed random numbers. This method is, on the CRAY-YMP computer, several times faster than the standard Box-Muller-Wiener algorithm. The validity of the approximation is discussed.