Function norm_zero_case [src]

Prototype

pub fn norm_zero_case(random: Random, u: f64) f64

Parameters

random: Randomu: f64

Source

  pub fn norm_zero_case(random: Random, u: f64) f64 {
    var x: f64 = 1;
    var y: f64 = 0;

    while (-2.0 * y < x * x) {
        x = @log(random.float(f64)) / norm_r;
        y = @log(random.float(f64));
    }

    if (u < 0) {
        return x - norm_r;
    } else {
        return norm_r - x;
    }
}  