average 0.17.0

Calculate statistics iteratively
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

use bencher::{benchmark_group, benchmark_main, Bencher};

/// Create a random vector by sampling from a normal distribution.
fn initialize_vec() -> Vec<f64> {
    use rand::SeedableRng;
    use rand_distr::{Distribution, Normal};
    let normal = Normal::new(2.0, 3.0).unwrap();
    let n = 1_000_000;
    let mut values = Vec::with_capacity(n);
    let mut rng = rand::rngs::Xoshiro256PlusPlus::from_seed([
        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,
    ]);
    for _ in 0..n {
        values.push(normal.sample(&mut rng));
    }
    values
}

fn bench_average(b: &mut Bencher) {
    let values = initialize_vec();
    b.iter(|| {
        let m: average::MeanWithError = values.iter().copied().collect();
        m
    });
}

fn bench_stats(b: &mut Bencher) {
    let values = initialize_vec();
    b.iter(|| {
        let m: stats::OnlineStats = values.iter().copied().collect();
        m
    });
}

benchmark_group!(benches, bench_average, bench_stats);
benchmark_main!(benches);