modppl 0.3.1

a experimental library for probabilistic programming in Rust.
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

mod pointed_model;
use nalgebra::dvector;
use pointed_model::types_2d;
use rand::rngs::ThreadRng;

use modppl::{Distribution, Real};

#[test]
fn test_uniform2d() {
    let mut rng = ThreadRng::default();
    let bounds = types_2d::Bounds {
        xmin: 0.,
        xmax: 2.5,
        ymin: -1.,
        ymax: 0.25,
    };

    let samples = (0..50000)
        .map(|_| types_2d::uniform_2d.random(&mut rng, bounds))
        .collect::<Vec<types_2d::Point>>();
    for sample in samples {
        assert!(sample[0] >= bounds.xmin);
        assert!(sample[0] <= bounds.xmax);
        assert!(sample[1] >= bounds.ymin);
        assert!(sample[1] <= bounds.ymax);
        approx::assert_abs_diff_eq!(
            types_2d::uniform_2d.logpdf(&sample, bounds),
            -1.1394342831883648,
            epsilon = Real::EPSILON
        );
    }

    assert_eq!(
        types_2d::uniform_2d.logpdf(&dvector![-1., 0.], bounds),
        Real::NEG_INFINITY
    );
}