rsl-interpolation 0.1.18

A re-write of GSL's Interpolation Routines 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
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

use crate::Accelerator;

fn setup_acc() -> Accelerator {
    Accelerator::new()
}

fn setup_xarray() -> [f64; 5] {
    [0.0, 1.0, 2.0, 3.0, 4.0]
}

#[test]
fn test_bsearch_interior_point() {
    let (acc, xarray) = (setup_acc(), setup_xarray());

    let res = acc.bsearch(&xarray, 1.5, 0, 4);
    assert_eq!(res, 1);
}

#[test]
fn test_bsearch_last_value() {
    let (acc, xarray) = (setup_acc(), setup_xarray());

    let res = acc.bsearch(&xarray, 4.0, 0, 4);
    assert_eq!(res, 3);
}

#[test]
fn test_bsearch_first_value() {
    let (acc, xarray) = (setup_acc(), setup_xarray());

    let res = acc.bsearch(&xarray, 0.0, 0, 4);
    assert_eq!(res, 0);
}

#[test]
fn test_bsearch_boundary() {
    let (acc, xarray) = (setup_acc(), setup_xarray());

    let res = acc.bsearch(&xarray, 2.0, 0, 4);
    assert_eq!(res, 2);
}

#[test]
fn test_bsearch_above_bounds() {
    let (acc, xarray) = (setup_acc(), setup_xarray());

    let res = acc.bsearch(&xarray, 10.0, 0, 4);
    assert_eq!(res, 3);
}

#[test]
fn test_bsearch_below_bounds() {
    let (acc, xarray) = (setup_acc(), setup_xarray());

    let res = acc.bsearch(&xarray, -10.0, 0, 4);
    assert_eq!(res, 0);
}

// TODO:
// #[test]
// fn bsearch_vs_lookup_bench() { }

#[test]
fn test_accelerator() {
    let xarray = setup_xarray();
    let mut acc = setup_acc();
    let mut k1 = 0;
    let mut k2 = 0;
    let mut t = false;
    let r = [
        -0.2, 0.0, 0.1, 0.7, 1.0, 1.3, 1.9, 2.0, 2.2, 2.7, 3.0, 3.1, 3.6, 4.0, 4.1, 4.9,
    ];

    // run through all the pairs of points.
    while (k1 < 16) & (k2 < 16) {
        let x = if t { r[k1] } else { r[k2] };
        t = !t;

        if !t {
            k1 = (k1 + 1) % 16;
            if k1 == 0 {
                k2 += 1;
            }
        }

        let i = acc.find(&xarray, x);
        let j = acc.bsearch(&xarray, x, 0, 4);
        assert_eq!(i, j);
    }
}