extern crate peroxide;
use peroxide::fuga::*;
#[test]
fn test_max_pool_1d() {
let w = WeightedUniform::from_max_pool_1d(f, (-2f64, 3f64), 10, 1e-5);
// w.weights().print();
println!("{:?}", w.intervals());
// w.sample(100000).mean().print();
// w.mean().print();
let x = linspace(w.intervals()[0].0, w.intervals()[8].1, 11);
let y1 = x.fmap(f);
let y2 = x.fmap(|t| w.weight_at(t));
for (a, b) in y1.iter().zip(y2.iter()) {
println!("a: {}, b: {}", a, b);
assert!(a <= b);
}
}
fn f(x: f64) -> f64 {
if x.abs() < 1f64 {
(1f64 - x.abs())
} else {
0f64
}
}