tessera 0.1.0

A library and command line interface for calculating geometric error for 3D Tiles tilesets
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

use crate::{error::TesseraError, geometry::PointPrimitive};

pub fn get_shortest_distance_between_points(
    a: &PointPrimitive,
    b: &PointPrimitive,
) -> Result<f64, TesseraError> {
    let mut shortest_distance = f64::INFINITY;

    for a_point in a.iter_vertices() {
        for b_point in b.iter_vertices() {
            let distance = point_distance_squared(&a_point, &b_point);

            if distance < shortest_distance {
                shortest_distance = distance;
            }
        }
    }

    return Ok(shortest_distance.sqrt());
}

pub fn point_distance_squared(a: &[f32; 3], b: &[f32; 3]) -> f64 {
    let dx: f64 = a[0] as f64 - b[0] as f64;
    let dy: f64 = a[1] as f64 - b[1] as f64;
    let dz: f64 = a[2] as f64 - b[2] as f64;

    return dx * dx + dy * dy + dz * dz;
}

#[cfg(test)]
mod tests {
    use crate::geometry::{PointPrimitive, Vertices};

    use super::*;

    #[test]
    fn test_get_shortest_distance_between_points() {
        let mut a = PointPrimitive::new();
        a.set_vertices(vec![[0.0, 0.0, 0.0]]);

        let mut b = PointPrimitive::new();
        b.set_vertices(vec![[1.0, 1.0, 1.0]]);

        let distance = get_shortest_distance_between_points(&a, &b);
        assert!(distance.is_ok());
        assert_eq!(distance.unwrap(), 3.0_f64.sqrt());
    }

    #[test]
    fn test_get_shortest_distance_between_same_points() {
        let mut a = PointPrimitive::new();
        a.set_vertices(vec![[1.0, 1.0, 1.0]]);

        let mut b = PointPrimitive::new();
        b.set_vertices(vec![[1.0, 1.0, 1.0]]);

        let distance = get_shortest_distance_between_points(&a, &b);
        assert!(distance.is_ok());
        assert_eq!(distance.unwrap(), 0.0_f64);
    }

    #[test]
    fn test_get_shortest_distance_between_points_with_multiple_points() {
        let mut a = PointPrimitive::new();
        a.set_vertices(vec![[0.0, 0.0, 0.0], [1.0, 1.0, 1.0]]);

        let mut b = PointPrimitive::new();
        b.set_vertices(vec![[2.0, 2.0, 2.0], [3.0, 3.0, 3.0]]);

        let distance = get_shortest_distance_between_points(&a, &b);
        assert!(distance.is_ok());
        assert_eq!(distance.unwrap(), 3.0_f64.sqrt());
    }

    #[test]
    fn test_get_shortest_distance_between_points_with_multiple_points_that_share_a_point() {
        let mut a = PointPrimitive::new();
        a.set_vertices(vec![[0.0, 0.0, 0.0], [1.0, 1.0, 1.0]]);

        let mut b = PointPrimitive::new();
        b.set_vertices(vec![[1.0, 1.0, 1.0], [2.0, 2.0, 2.0]]);

        let distance = get_shortest_distance_between_points(&a, &b);
        assert!(distance.is_ok());
        assert_eq!(distance.unwrap(), 0.0_f64);
    }

    #[test]
    fn test_point_distance_squared() {
        let a = [0.0, 0.0, 0.0];
        let b = [1.0, 1.0, 1.0];
        let distance = point_distance_squared(&a, &b);
        assert_eq!(distance, 3.0);
    }
}