oxiproj-core 0.1.0

Foundation types for OxiProj: coordinates, errors, ellipsoids, datums, and units.
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

//! `pj_tsfn` ported from PROJ `src/tsfn.cpp`.

/// Ported from src/tsfn.cpp (`pj_tsfn`).
///
/// Computes the small `t` function used by conformal projections,
/// `exp(e*atanh(e*sin(phi))) * (cos(phi)/(1+sin(phi)))`, with a numerically
/// stable branch selected on the sign of `sin(phi)`.
pub fn pj_tsfn(phi: f64, sinphi: f64, e: f64) -> f64 {
    let cosphi = phi.cos();
    (e * (e * sinphi).atanh()).exp()
        * if sinphi > 0.0 {
            cosphi / (1.0 + sinphi)
        } else {
            (1.0 - sinphi) / cosphi
        }
}

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

    #[test]
    fn tsfn_sphere_matches_closed_form() {
        let phi = 0.5_f64;
        let t = pj_tsfn(phi, phi.sin(), 0.0);
        // exp(0) == 1, so only the trig branch remains.
        let expected = phi.cos() / (1.0 + phi.sin());
        assert!((t - expected).abs() < 1e-12);
    }

    #[test]
    fn tsfn_wgs84_is_positive() {
        let phi = 0.6_f64;
        let e = 0.0818191908426;
        let t = pj_tsfn(phi, phi.sin(), e);
        assert!(t > 0.0);
    }
}