r2rs-nmath 0.1.1

Statistics programming for Rust based on R's nmath package
Documentation 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52

mod lecuyer_cmrg;
mod marsaglia_multicarry;
mod mersenne_twister;
mod sexp;
mod snorm;
mod super_duper;
mod wichmann_hill;

pub use self::{
    lecuyer_cmrg::*, marsaglia_multicarry::*, mersenne_twister::*, super_duper::*, wichmann_hill::*,
};
use self::{sexp::exp_rand, snorm::norm_rand};
use crate::traits::RNG;

pub trait NMathRNG {
    fn norm_rand(&mut self) -> f64;
    fn exp_rand(&mut self) -> f64;
}

impl<R> NMathRNG for R
where
    R: RNG + Sized,
{
    fn norm_rand(&mut self) -> f64 {
        norm_rand(self)
    }

    fn exp_rand(&mut self) -> f64 {
        exp_rand(self)
    }
}

static i2_32m1: f64 = 2.328306437080797e-10;

pub fn fixup(x: f64) -> f64 {
    if x <= 0.0 {
        0.5 * i2_32m1
    } else if (1.0 - x) <= 0.0 {
        1.0 - 0.5 * i2_32m1
    } else {
        x
    }
}

#[cfg(test)]
mod tests;

#[cfg(all(test, feature = "enable_proptest"))]
mod proptests;

#[cfg(all(test, feature = "enable_covtest"))]
mod covtests;