2017-05-30 10:54:14 +02:00
|
|
|
#![cfg_attr(feature = "cargo-clippy", allow(float_cmp, map_clone))]
|
|
|
|
|
2017-05-19 15:54:13 +02:00
|
|
|
use core::iter::Iterator;
|
|
|
|
|
2019-07-08 16:17:04 +02:00
|
|
|
use average::{MeanWithError, Estimate, Merge, assert_almost_eq};
|
2017-05-19 15:54:13 +02:00
|
|
|
|
|
|
|
#[test]
|
|
|
|
fn trivial() {
|
2017-05-28 21:13:47 +02:00
|
|
|
let mut a = MeanWithError::new();
|
2017-05-19 15:54:13 +02:00
|
|
|
assert_eq!(a.len(), 0);
|
|
|
|
a.add(1.0);
|
|
|
|
assert_eq!(a.mean(), 1.0);
|
|
|
|
assert_eq!(a.len(), 1);
|
|
|
|
assert_eq!(a.sample_variance(), 0.0);
|
|
|
|
assert_eq!(a.population_variance(), 0.0);
|
|
|
|
assert_eq!(a.error(), 0.0);
|
|
|
|
a.add(1.0);
|
|
|
|
assert_eq!(a.mean(), 1.0);
|
|
|
|
assert_eq!(a.len(), 2);
|
|
|
|
assert_eq!(a.sample_variance(), 0.0);
|
|
|
|
assert_eq!(a.population_variance(), 0.0);
|
|
|
|
assert_eq!(a.error(), 0.0);
|
2017-05-05 16:40:23 +02:00
|
|
|
}
|
|
|
|
|
|
|
|
#[test]
|
2017-05-19 15:54:13 +02:00
|
|
|
fn simple() {
|
2017-05-28 21:13:47 +02:00
|
|
|
let a: MeanWithError = (1..6).map(f64::from).collect();
|
2017-05-19 15:54:13 +02:00
|
|
|
assert_eq!(a.mean(), 3.0);
|
|
|
|
assert_eq!(a.len(), 5);
|
|
|
|
assert_eq!(a.sample_variance(), 2.5);
|
|
|
|
assert_almost_eq!(a.error(), f64::sqrt(0.5), 1e-16);
|
2017-05-05 16:40:23 +02:00
|
|
|
}
|
|
|
|
|
2019-01-22 14:29:29 +01:00
|
|
|
#[cfg(feature = "serde1")]
|
2017-12-20 22:46:50 +01:00
|
|
|
#[test]
|
|
|
|
fn simple_serde() {
|
|
|
|
let a: MeanWithError = (1..6).map(f64::from).collect();
|
|
|
|
let b = serde_json::to_string(&a).unwrap();
|
|
|
|
assert_eq!(&b, "{\"avg\":{\"avg\":3.0,\"n\":5},\"sum_2\":10.0}");
|
|
|
|
let c: MeanWithError = serde_json::from_str(&b).unwrap();
|
|
|
|
assert_eq!(c.mean(), 3.0);
|
|
|
|
assert_eq!(c.len(), 5);
|
|
|
|
assert_eq!(c.sample_variance(), 2.5);
|
|
|
|
assert_almost_eq!(c.error(), f64::sqrt(0.5), 1e-16);
|
|
|
|
}
|
|
|
|
|
2017-05-05 16:40:23 +02:00
|
|
|
#[test]
|
2017-05-19 15:54:13 +02:00
|
|
|
fn numerically_unstable() {
|
|
|
|
// The naive algorithm fails for this example due to cancelation.
|
|
|
|
let big = 1e9;
|
|
|
|
let sample = &[big + 4., big + 7., big + 13., big + 16.];
|
2018-02-27 01:38:45 +01:00
|
|
|
let a: MeanWithError = sample.iter().collect();
|
2017-05-19 15:54:13 +02:00
|
|
|
assert_eq!(a.sample_variance(), 30.);
|
|
|
|
}
|
|
|
|
|
2017-05-24 10:48:27 +02:00
|
|
|
#[test]
|
|
|
|
fn merge() {
|
|
|
|
let sequence: &[f64] = &[1., 2., 3., 4., 5., 6., 7., 8., 9.];
|
|
|
|
for mid in 0..sequence.len() {
|
|
|
|
let (left, right) = sequence.split_at(mid);
|
2018-02-27 01:38:45 +01:00
|
|
|
let avg_total: MeanWithError = sequence.iter().collect();
|
|
|
|
let mut avg_left: MeanWithError = left.iter().collect();
|
|
|
|
let avg_right: MeanWithError = right.iter().collect();
|
2017-05-24 10:48:27 +02:00
|
|
|
avg_left.merge(&avg_right);
|
|
|
|
assert_eq!(avg_total.len(), avg_left.len());
|
|
|
|
assert_eq!(avg_total.mean(), avg_left.mean());
|
|
|
|
assert_eq!(avg_total.sample_variance(), avg_left.sample_variance());
|
|
|
|
}
|
|
|
|
}
|