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
68
69
70
71
72
73

use sciforge::hub::prelude::*;

fn run_phys(name: &str, params: Vec<(&str, ParameterValue)>) -> RunOutput {
    let mut exp = Experiment::new(DomainType::Physics, 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 snell() {
    let v = scalar(run_phys(
        "optics::refraction::snell",
        vec![
            ("n1", ParameterValue::Scalar(1.0)),
            ("n2", ParameterValue::Scalar(1.5)),
            ("theta1", ParameterValue::Scalar(0.5236)),
        ],
    ));
    assert!(
        v > 0.0 && v < 0.5236,
        "refracted angle smaller in denser medium, got {v}"
    );
}

#[test]
fn thin_lens_equation() {
    let v = scalar(run_phys(
        "thin_lens_equation",
        vec![
            ("focal", ParameterValue::Scalar(0.1)),
            ("object_dist", ParameterValue::Scalar(0.2)),
        ],
    ));
    assert!((v - 0.2).abs() < 1e-8, "2f object => 2f image, got {v}");
}

#[test]
fn brewster_angle() {
    let v = scalar(run_phys(
        "optics::refraction::brewster_angle",
        vec![
            ("n1", ParameterValue::Scalar(1.0)),
            ("n2", ParameterValue::Scalar(1.5)),
        ],
    ));
    assert!(
        v > 0.5 && v < 1.2,
        "Brewster angle for glass ~0.98 rad, got {v}"
    );
}

#[test]
fn rayleigh_criterion() {
    let v = scalar(run_phys(
        "rayleigh_criterion",
        vec![
            ("wavelength", ParameterValue::Scalar(500e-9)),
            ("aperture", ParameterValue::Scalar(0.01)),
        ],
    ));
    assert!(v > 0.0, "angular resolution should be positive, got {v}");
}