petekio 0.2.8

Subsurface data ingestion + structure layer: surfaces, wells, points, polygons with loading, interpolation, and statistics.
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

//! Length units. petekIO never guesses a unit — it is carried on the project
//! and conversions live here.

/// A length unit. Carried on a `GeoData` project; surfaces/wells inherit it.
#[derive(Debug, Clone, Copy, PartialEq, Eq, serde::Serialize, serde::Deserialize)]
pub enum Unit {
    Feet,
    Metres,
}

impl Unit {
    /// Metres per one of `self`.
    pub fn metres_per_unit(self) -> f64 {
        match self {
            Unit::Feet => 0.3048,
            Unit::Metres => 1.0,
        }
    }

    /// Convert `value` (in `self` units) to `to` units.
    pub fn convert(self, value: f64, to: Unit) -> f64 {
        value * self.metres_per_unit() / to.metres_per_unit()
    }

    /// Convert a planar area expressed in `self`² (e.g. m² or ft²) to acres.
    /// One acre = 4046.8564224 m². Backs `ModelInputs::reservoir_area_acres`.
    pub fn area_to_acres(self, area_in_unit_sq: f64) -> f64 {
        const SQ_METRES_PER_ACRE: f64 = 4_046.856_422_4;
        let m = self.metres_per_unit();
        area_in_unit_sq * m * m / SQ_METRES_PER_ACRE
    }
}

#[cfg(test)]
mod tests {
    use super::*;
    use approx::assert_relative_eq;

    #[test]
    fn feet_to_metres_roundtrip() {
        assert_relative_eq!(Unit::Feet.convert(100.0, Unit::Metres), 30.48);
        assert_relative_eq!(Unit::Metres.convert(30.48, Unit::Feet), 100.0);
        assert_relative_eq!(Unit::Feet.convert(42.0, Unit::Feet), 42.0);
    }

    #[test]
    fn area_to_acres_hand_calc() {
        // 43560 ft² = 1 acre exactly.
        assert_relative_eq!(Unit::Feet.area_to_acres(43_560.0), 1.0, epsilon = 1e-9);
        // 4046.8564224 m² = 1 acre.
        assert_relative_eq!(
            Unit::Metres.area_to_acres(4_046.856_422_4),
            1.0,
            epsilon = 1e-9
        );
    }
}