point_process 0.2.3

A crate for simulating random point processes.
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
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

extern crate pointprocesses;
extern crate gnuplot;

#[macro_use]
extern crate ndarray;

use gnuplot::{Figure,Caption,Color,PointSymbol,PointSize};

use pointprocesses::generalized;

fn main() {
    
    square_example();

    circle_example();

}

fn square_example() {
    let lambda = 50.0;

    let close = array![0.0,0.0];
    let far = array![1.0,1.0];
    let domain = generalized::Rectangle::new(close, far);

    let events = generalized::poisson_process(lambda, &domain);

    println!("{:?}", events);

    let mut fg = Figure::new();

    fg.axes2d()
        .points(
            &events.slice(s![..,0]),
            &events.slice(s![..,1]),
            &[
                Caption("Events"),
                Color("red"),
                PointSymbol('O'),
                PointSize(0.8)
            ]
        );
    
    fg.echo_to_file("2d_poisson_square.gnuplot");
}

fn circle_example() {
    let lambda = 500.0;

    let center = array![1.0,1.0];
    let radius = 1.0;
    let domain = generalized::Sphere::new(center, radius);

    let events = generalized::poisson_process(lambda, &domain);

    println!("{:?}", events);

    let mut fg = Figure::new();

    fg.axes2d()
        .points(
            &events.slice(s![..,0]),
            &events.slice(s![..,1]),
            &[
                Caption("Events"),
                Color("red"),
                PointSymbol('O'),
                PointSize(0.8)
            ]
        );
    
    fg.echo_to_file("2d_poisson_circle.gnuplot");
    fg.show();
}