double Diffusion::get_new_pos(const Walker* walker, int i, int j) { #ifdef RNG_ZIG return DRanNormalZig32() * std; #endif #ifdef RNG_NUMREC return gaussian_deviate(&random_seed) * std; #endif }
double Diffusion::get_new_pos(const Walker* walker, int i, int j) { return gaussian_deviate(&random_seed) * std; }