bacon-sci 0.16.2

Scientific computing 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

use crate::differentiate::*;

fn exp(f: f64) -> f64 {
    f.exp()
}

fn quadratic(f: f64) -> f64 {
    f.powi(2)
}

fn sin(f: f64) -> f64 {
    f.sin()
}

#[test]
fn deriv() {
    for i in -1000..1000 {
        let i = i as f64 * 0.01;
        let d = derivative(exp, i, 0.1);
        assert!(((d - i.exp()) / i.exp()).abs() < 0.01);
    }

    for i in -1000..1000 {
        let i = i as f64 * 0.01;
        let d = derivative(quadratic, i, 0.1);
        if 2.0 * i.abs() > 0.05 {
            assert!(((d - 2.0 * i) / (2.0 * i)).abs() < 0.05);
        } else {
            assert!((d - 2.0 * i).abs() < 0.05);
        }
    }

    for i in -1000..1000 {
        let i = i as f64 * 0.01;
        let d = derivative(sin, i, 0.1);
        if i.cos().abs() > 0.05 {
            assert!(((d - i.cos()) / i.cos()).abs() < 0.01);
        } else {
            assert!((d - i.cos()).abs() < 0.05);
        }
    }
}

#[test]
fn second_deriv() {
    for i in -1000..1000 {
        let i = i as f64 * 0.01;
        let d = second_derivative(exp, i, 0.1);
        assert!(((d - i.exp()) / i.exp()).abs() < 0.01);
    }

    for i in -1000..1000 {
        let i = i as f64 * 0.01;
        let d = second_derivative(quadratic, i, 0.1);
        assert!(((d - 2.0) / 2.0).abs() < 0.05);
    }

    for i in -1000..1000 {
        let i = i as f64 * 0.01;
        let d = second_derivative(sin, i, 0.1);
        if i.sin().abs() > 0.05 {
            assert!(((d + i.sin()) / i.sin()).abs() < 0.01);
        } else {
            assert!((d + i.sin()).abs() < 0.05);
        }
    }
}