flavio 0.5.0

flavio welcomes you
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

use super::{NewtonRaphson, SecondOrder, TensorRank0};

const TOLERANCE: TensorRank0 = 1e-6;

#[test]
fn linear() {
    assert!(
        NewtonRaphson {
            ..Default::default()
        }
        .minimize(
            |x: &TensorRank0| Ok(*x),
            |_: &TensorRank0| Ok(1.0),
            1.0,
            None,
            None
        )
        .unwrap()
        .abs()
            < TOLERANCE
    )
}

#[test]
fn quadratic() {
    assert!(
        NewtonRaphson {
            ..Default::default()
        }
        .minimize(
            |x: &TensorRank0| Ok(x.powi(2) / 2.0),
            |x: &TensorRank0| Ok(*x),
            1.0,
            None,
            None
        )
        .unwrap()
        .abs()
            < TOLERANCE
    )
}

#[test]
fn sin() {
    assert!(
        NewtonRaphson {
            ..Default::default()
        }
        .minimize(
            |x: &TensorRank0| Ok(x.sin()),
            |x: &TensorRank0| Ok(x.cos()),
            1.0,
            None,
            None
        )
        .unwrap()
        .abs()
            < TOLERANCE
    )
}

#[test]
#[should_panic(expected = "The obtained solution is not a minimum.")]
fn sin_max() {
    NewtonRaphson {
        ..Default::default()
    }
    .minimize(
        |x: &TensorRank0| Ok(x.sin()),
        |x: &TensorRank0| Ok(x.cos()),
        3.0,
        None,
        None,
    )
    .unwrap();
}