sciforge 0.0.3

A comprehensive scientific computing library in pure Rust with zero dependencies
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 sciforge::hub::prelude::*;

fn run_maths(name: &str, params: Vec<(&str, ParameterValue)>) -> RunOutput {
    let mut exp = Experiment::new(DomainType::Maths, name);
    for (k, v) in params {
        exp = exp.param(k, v);
    }
    ExperimentRunner::new()
        .run(&exp)
        .unwrap_or_else(|_| panic!("dispatch '{name}' failed"))
}

fn scalar(o: RunOutput) -> f64 {
    match o {
        RunOutput::Scalar(v) => v,
        _ => panic!("expected Scalar, got {o:?}"),
    }
}

#[test]
fn quadratic_objective() {
    let v = scalar(run_maths(
        "quadratic_objective",
        vec![
            (
                "h",
                ParameterValue::Matrix(vec![vec![2.0, 0.0], vec![0.0, 2.0]]),
            ),
            ("c", ParameterValue::Vector(vec![0.0, 0.0])),
            ("x", ParameterValue::Vector(vec![0.0, 0.0])),
        ],
    ));
    assert!(v.abs() < 1e-8, "f([0,0])=0 at minimum, got {v}");
}

#[test]
fn project_box() {
    let o = run_maths(
        "project_box",
        vec![
            ("x", ParameterValue::Vector(vec![5.0])),
            ("lower", ParameterValue::Vector(vec![0.0])),
            ("upper", ParameterValue::Vector(vec![3.0])),
        ],
    );
    match o {
        RunOutput::Vector(v) => assert!(
            (v[0] - 3.0).abs() < 1e-8,
            "project(5, [0,3])=3, got {}",
            v[0]
        ),
        other => panic!("expected Vector, got {other:?}"),
    }
}

#[test]
fn lagrangian() {
    let v = scalar(run_maths(
        "lagrangian",
        vec![
            ("f", ParameterValue::Scalar(10.0)),
            ("constraints", ParameterValue::Vector(vec![2.0])),
            ("lambdas", ParameterValue::Vector(vec![3.0])),
        ],
    ));
    assert!((v - 16.0).abs() < 1e-8, "L=f+lambda*g=10+3*2=16, got {v}");
}