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
 53
 54
 55
 56
 57
 58
 59
 60
 61
 62
 63
 64
 65
 66
 67
 68
 69
 70
 71
 72
 73
 74
 75
 76
 77
 78
 79
 80
 81
 82
 83
 84
 85
 86
 87
 88
 89
 90
 91
 92
 93
 94
 95
 96
 97
 98
 99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122

#![feature(
    box_syntax, fn_traits, unboxed_closures, pattern_parentheses, test, specialization,
    optin_builtin_traits
)]
#![allow(improper_ctypes)]

extern crate colored;
extern crate gnuplot;
extern crate itertools_num;
extern crate libc;
extern crate rand;
extern crate test;
#[macro_use]
extern crate log;

mod canvas;
mod minuit;
mod terms;

pub use canvas::Canvas;
pub use terms::{function::Function, graph::Graph, histogram::Histogram};

#[macro_export]
macro_rules! fit {
    ($func:expr, $data:expr) => {
        minuit()
            .fit($func, $data)
            .map(|r| Function::from(($func, r)))
            .map_err(|r| Function::from(($func, r)))
    };
}

pub mod fileio;
pub mod utility;
pub mod prelude {
    pub use self::Dir::*;
    pub use minuit::minuit;

    #[derive(Debug, Clone, Copy, PartialEq)]
    pub enum Dir {
        Right,
        Left,
    }

    #[derive(Debug, Clone, Copy, PartialEq)]
    pub struct ParamInfo {
        pub err: f64,
        pub rng: Option<(f64, f64)>,
    }

    impl ParamInfo {
        pub fn new(err: f64, rng: Option<(f64, f64)>) -> Self {
            assert!(err > 0f64);
            ParamInfo { err, rng }
        }
    }

    #[derive(Debug, Clone, Copy, PartialEq)]
    pub struct Data {
        pub x: f64,
        pub val: f64,
        pub valerr: f64,
    }

    pub trait WithData {
        fn data(&self) -> Vec<Data>;
    }

    #[derive(Debug, Clone, Copy, PartialEq)]
    pub enum MnOut<T> {
        Stdout,
        File(T),
        None,
    }
}

#[cfg(test)]
mod tests {
    use super::utility::*;
    use super::*;
    use rand::{Rng, SeedableRng, StdRng};
    use std::time::{self, UNIX_EPOCH};
    use test::Bencher;

    fn randoms<R: Rng>(n: usize, rng: R) -> f64 {
        RngDist::new(rng, -15.0, 15.0, &|x| {
            (-0.5 * (x - 0.0).powi(2) / (3f64 * 3.0)).exp()
        }).squize(n)
            .take(10000)
            .sum()
    }

    #[bench]
    fn rnd10(b: &mut Bencher) {
        let rng = StdRng::from_seed(&[time::SystemTime::now()
            .duration_since(UNIX_EPOCH)
            .unwrap()
            .as_secs() as usize]);

        b.iter(|| randoms(10, rng))
    }

    #[bench]
    fn rnd100(b: &mut Bencher) {
        let rng = StdRng::from_seed(&[time::SystemTime::now()
            .duration_since(UNIX_EPOCH)
            .unwrap()
            .as_secs() as usize]);

        b.iter(|| randoms(100, rng))
    }

    #[bench]
    fn rnd1000(b: &mut Bencher) {
        let rng = StdRng::from_seed(&[time::SystemTime::now()
            .duration_since(UNIX_EPOCH)
            .unwrap()
            .as_secs() as usize]);

        b.iter(|| randoms(1000, rng))
    }
}