Enable more aggressive optimisation of PRN generation.

Implement Box-Muller method for generating normal distributed
PRN. This seems to be more precise but is currently slightly
slower than the central limit theorem method.