The Ultimate Box-Muller transform - American History Information Guide and Reference

A Box-Muller transform is a method of generating pairs of independent standard normally distributed (zero expectation, unit variance) random numbers, given a source of uniformly distributed random numbers. There are two kinds:

(1) Given r and φ independently uniformly distributed in (0,1], compute:

$z_0 = \cos(2 \pi \varphi) \cdot \sqrt{-2 \ln r}$

and

$z_1 = \sin(2 \pi \varphi) \cdot \sqrt{-2 \ln r}.$

(2) Given x and y independently uniformly distributed in [−1,1], set R = x² + y². If R = 0 or R > 1, throw them away and try another pair (x, y). Then, for these filtered points, compute:

$z_0 = x \cdot \sqrt{\frac{-2 \ln R}{R}}$

and

$z_1 = y \cdot \sqrt{\frac{-2 \ln R}{R}}.$

The second method is typically faster because it uses only one transcendental function instead of at least two, even though it throws away 1 − π/4 ≈ 21.46% of the total input uniformly distributed random number pairs generated, i.e. throws away 4/π − 1 ≈ 0.2732 uniformly distributed random number pairs per Gaussian random number pair generated, requiring 4/π ≈ 1.2732 input random numbers per output random number.

External links

Generating Gaussian Random Numbers

Categories: Transforms | Sampling techniques

Your American History Reference Guide! - Box-Muller transform

Box-Muller transform

External links

Your American History Reference Guide!
- Box-Muller transform