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
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155

//! Dispatch handler for gas law functions.

use super::super::params::*;
use crate::hub::domain::chemistry as chem;
use crate::hub::domain::common::errors::{HubError, HubResult};
use crate::hub::engine::experience::runner::RunOutput;

pub(super) fn dispatch(func: &str, p: &Params) -> HubResult<RunOutput> {
    match func {
        "ideal_gas_pressure" => Ok(RunOutput::Scalar(
            chem::gas_laws::ideal::ideal_gas_pressure(
                get_f(p, "n")?,
                get_f(p, "t")?,
                get_f(p, "v")?,
            ),
        )),
        "ideal_gas_volume" => Ok(RunOutput::Scalar(chem::gas_laws::ideal::ideal_gas_volume(
            get_f(p, "n")?,
            get_f(p, "t")?,
            get_f(p, "p")?,
        ))),
        "ideal_gas_temperature" => Ok(RunOutput::Scalar(
            chem::gas_laws::ideal::ideal_gas_temperature(
                get_f(p, "p")?,
                get_f(p, "v")?,
                get_f(p, "n")?,
            ),
        )),
        "boyles_law" => Ok(RunOutput::Scalar(chem::gas_laws::ideal::boyles_law(
            get_f(p, "p1")?,
            get_f(p, "v1")?,
            get_f(p, "v2")?,
        ))),
        "charles_law" => Ok(RunOutput::Scalar(chem::gas_laws::ideal::charles_law(
            get_f(p, "v1")?,
            get_f(p, "t1")?,
            get_f(p, "t2")?,
        ))),
        "gay_lussac_law" => Ok(RunOutput::Scalar(chem::gas_laws::ideal::gay_lussac_law(
            get_f(p, "p1")?,
            get_f(p, "t1")?,
            get_f(p, "t2")?,
        ))),
        "combined_gas_law" => Ok(RunOutput::Scalar(chem::gas_laws::ideal::combined_gas_law(
            get_f(p, "p1")?,
            get_f(p, "v1")?,
            get_f(p, "t1")?,
            get_f(p, "t2")?,
            get_f(p, "p2")?,
        ))),
        "dalton_partial_pressure" => Ok(RunOutput::Scalar(
            chem::gas_laws::ideal::dalton_partial_pressure(
                get_f(p, "mole_fraction")?,
                get_f(p, "total_pressure")?,
            ),
        )),
        "grahams_law_effusion" => Ok(RunOutput::Scalar(
            chem::gas_laws::ideal::grahams_law_effusion(get_f(p, "m1")?, get_f(p, "m2")?),
        )),
        "gas_density" => Ok(RunOutput::Scalar(chem::gas_laws::ideal::gas_density(
            get_f(p, "p")?,
            get_f(p, "mw")?,
            get_f(p, "t")?,
        ))),
        "rms_speed" => Ok(RunOutput::Scalar(chem::gas_laws::ideal::rms_speed(
            get_f(p, "t")?,
            get_f(p, "mw")?,
        ))),
        "mean_speed" => Ok(RunOutput::Scalar(chem::gas_laws::ideal::mean_speed(
            get_f(p, "t")?,
            get_f(p, "mw")?,
        ))),
        "most_probable_speed" => Ok(RunOutput::Scalar(
            chem::gas_laws::ideal::most_probable_speed(get_f(p, "t")?, get_f(p, "mw")?),
        )),
        "mean_free_path" => Ok(RunOutput::Scalar(chem::gas_laws::ideal::mean_free_path(
            get_f(p, "d")?,
            get_f(p, "n_over_v")?,
        ))),

        "van_der_waals_pressure" => Ok(RunOutput::Scalar(
            chem::gas_laws::real::van_der_waals_pressure(
                get_f(p, "n")?,
                get_f(p, "t")?,
                get_f(p, "v")?,
                get_f(p, "a")?,
                get_f(p, "b")?,
            ),
        )),
        "redlich_kwong_pressure" => Ok(RunOutput::Scalar(
            chem::gas_laws::real::redlich_kwong_pressure(
                get_f(p, "n")?,
                get_f(p, "t")?,
                get_f(p, "v")?,
                get_f(p, "a")?,
                get_f(p, "b")?,
            ),
        )),
        "compressibility_factor" => Ok(RunOutput::Scalar(
            chem::gas_laws::real::compressibility_factor(
                get_f(p, "p")?,
                get_f(p, "v")?,
                get_f(p, "n")?,
                get_f(p, "t")?,
            ),
        )),
        "virial_equation_b" => Ok(RunOutput::Scalar(chem::gas_laws::real::virial_equation_b(
            get_f(p, "p")?,
            get_f(p, "t")?,
            get_f(p, "b")?,
        ))),
        "boyle_temperature" => Ok(RunOutput::Scalar(chem::gas_laws::real::boyle_temperature(
            get_f(p, "a")?,
            get_f(p, "b")?,
        ))),
        "critical_temperature" => Ok(RunOutput::Scalar(
            chem::gas_laws::real::critical_temperature(get_f(p, "a")?, get_f(p, "b")?),
        )),
        "critical_pressure" => Ok(RunOutput::Scalar(chem::gas_laws::real::critical_pressure(
            get_f(p, "a")?,
            get_f(p, "b")?,
        ))),
        "critical_volume" => Ok(RunOutput::Scalar(chem::gas_laws::real::critical_volume(
            get_f(p, "b")?,
        ))),
        "reduced_temperature" => Ok(RunOutput::Scalar(
            chem::gas_laws::real::reduced_temperature(get_f(p, "t")?, get_f(p, "tc")?),
        )),
        "reduced_pressure" => Ok(RunOutput::Scalar(chem::gas_laws::real::reduced_pressure(
            get_f(p, "p")?,
            get_f(p, "pc")?,
        ))),
        "peng_robinson_pressure" => Ok(RunOutput::Scalar(
            chem::gas_laws::real::peng_robinson_pressure(
                get_f(p, "t")?,
                get_f(p, "vm")?,
                get_f(p, "a")?,
                get_f(p, "b")?,
            ),
        )),
        "fugacity_coefficient" => Ok(RunOutput::Scalar(
            chem::gas_laws::real::fugacity_coefficient(
                get_f(p, "z")?,
                get_f(p, "b_prime")?,
                get_f(p, "p")?,
            ),
        )),
        "acentric_factor" => Ok(RunOutput::Scalar(chem::gas_laws::real::acentric_factor(
            get_f(p, "p_sat")?,
            get_f(p, "pc")?,
        ))),

        _ => Err(HubError::InvalidInput(format!("unknown function: {func}"))),
    }
}