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

//! Dispatch handler for environmental chemistry 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 {
        "ozone_formation_rate" => Ok(RunOutput::Scalar(
            chem::environmental::atmosphere::ozone_formation_rate(
                get_f(p, "k")?,
                get_f(p, "no2")?,
                get_f(p, "hv")?,
            ),
        )),
        "ozone_destruction_rate" => Ok(RunOutput::Scalar(
            chem::environmental::atmosphere::ozone_destruction_rate(
                get_f(p, "k")?,
                get_f(p, "o3")?,
                get_f(p, "no")?,
            ),
        )),
        "photochemical_smog_potential" => Ok(RunOutput::Scalar(
            chem::environmental::atmosphere::photochemical_smog_potential(
                get_f(p, "voc")?,
                get_f(p, "nox")?,
            ),
        )),
        "atmospheric_lifetime" => Ok(RunOutput::Scalar(
            chem::environmental::atmosphere::atmospheric_lifetime(
                get_f(p, "concentration")?,
                get_f(p, "removal_rate")?,
            ),
        )),
        "global_warming_potential" => Ok(RunOutput::Scalar(
            chem::environmental::atmosphere::global_warming_potential(
                get_f(p, "rf_gas")?,
                get_f(p, "lifetime_gas")?,
                get_f(p, "rf_co2")?,
                get_f(p, "lifetime_co2")?,
            ),
        )),
        "henry_law_solubility" => Ok(RunOutput::Scalar(
            chem::environmental::atmosphere::henry_law_solubility(
                get_f(p, "kh")?,
                get_f(p, "partial_pressure")?,
            ),
        )),
        "ozone_depletion_potential" => Ok(RunOutput::Scalar(
            chem::environmental::atmosphere::ozone_depletion_potential(
                get_f(p, "cl_atoms")?,
                get_f(p, "lifetime")?,
                get_f(p, "mw")?,
            ),
        )),
        "photolysis_rate_constant" => Ok(RunOutput::Scalar(
            chem::environmental::atmosphere::photolysis_rate_constant(
                get_f(p, "quantum_yield")?,
                get_f(p, "absorption")?,
                get_f(p, "actinic_flux")?,
            ),
        )),

        "biochemical_oxygen_demand" => Ok(RunOutput::Scalar(
            chem::environmental::water::biochemical_oxygen_demand(
                get_f(p, "bod_ultimate")?,
                get_f(p, "k")?,
                get_f(p, "t")?,
            ),
        )),
        "chemical_oxygen_demand" => Ok(RunOutput::Scalar(
            chem::environmental::water::chemical_oxygen_demand(
                get_f(p, "sample_oxygen")?,
                get_f(p, "blank_oxygen")?,
                get_f(p, "volume")?,
            ),
        )),
        "dissolved_oxygen_saturation" => Ok(RunOutput::Scalar(
            chem::environmental::water::dissolved_oxygen_saturation(get_f(p, "t")?),
        )),
        "streeter_phelps" => Ok(RunOutput::Scalar(
            chem::environmental::water::streeter_phelps(
                get_f(p, "d0")?,
                get_f(p, "l0")?,
                get_f(p, "kd")?,
                get_f(p, "kr")?,
                get_f(p, "t")?,
            ),
        )),
        "critical_point_time" => Ok(RunOutput::Scalar(
            chem::environmental::water::critical_point_time(
                get_f(p, "kd")?,
                get_f(p, "kr")?,
                get_f(p, "d0")?,
                get_f(p, "l0")?,
            ),
        )),
        "chlorine_decay" => Ok(RunOutput::Scalar(
            chem::environmental::water::chlorine_decay(
                get_f(p, "c0")?,
                get_f(p, "k")?,
                get_f(p, "t")?,
            ),
        )),
        "ct_disinfection" => Ok(RunOutput::Scalar(
            chem::environmental::water::ct_disinfection(get_f(p, "c")?, get_f(p, "t")?),
        )),
        "hardness_total" => Ok(RunOutput::Scalar(
            chem::environmental::water::hardness_total(get_f(p, "ca_mg_l")?, get_f(p, "mg_mg_l")?),
        )),
        "langelier_saturation_index" => Ok(RunOutput::Scalar(
            chem::environmental::water::langelier_saturation_index(
                get_f(p, "ph")?,
                get_f(p, "ph_s")?,
            ),
        )),
        "total_dissolved_solids_from_conductivity" => Ok(RunOutput::Scalar(
            chem::environmental::water::total_dissolved_solids_from_conductivity(get_f(
                p,
                "conductivity_us",
            )?),
        )),

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