rcgal 0.2.2

Rust Computational Geometry Algorithms Library.
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

use crate::kernel::{number_type::NumberType, point_2::Point2, segment_2::Segment2};

use super::location_enum::Point2Segment2Location;

pub fn is_point_2_on_line_segment_2<T: NumberType>(
    point_2: &Point2<T>,
    segment_2: &impl Segment2<T>,
) -> bool {
    let location = locate_point_2_line_segment_2(point_2, segment_2);
    match location {
        Point2Segment2Location::On => true,
        _ => false,
    }
}

pub fn locate_point_2_line_segment_2<T: NumberType>(
    point_2: &Point2<T>,
    segment_2: &impl Segment2<T>,
) -> Point2Segment2Location {
    let source = segment_2.source();
    let target = segment_2.target();
    let vec_ab = target - source;
    let vec_ac = *point_2 - source;
    let area = vec_ab.cross(&vec_ac);
    if area > T::default_eps() {
        Point2Segment2Location::Left
    } else if area < -T::default_eps() {
        Point2Segment2Location::Right
    } else {
        let dot_ab_ac = vec_ab.dot(&vec_ac);
        let dot_ab_ab = vec_ab.dot(&vec_ab);
        if dot_ab_ac <= dot_ab_ab && dot_ab_ac >= -T::default_eps() {
            Point2Segment2Location::On
        } else {
            Point2Segment2Location::Collinear
        }
    }
}

#[cfg(test)]
mod tests {

    use crate::kernel::line_segment_2::LineSegment2;

    use super::*;

    #[test]
    fn test_is_point_2_on_line_segment_2() {
        let p1 = Point2::new(0.0, 0.0);
        let p2 = Point2::new(10.0, 10.0);
        let s = LineSegment2::new(p1, p2);

        assert_eq!(is_point_2_on_line_segment_2(&p1, &s), true);

        let p1 = Point2::new(0.0, 0.0);
        let p2 = Point2::new(10.0, 10.0);
        let s = LineSegment2::new(p1, p2);

        let p3 = Point2::new(5.0, 5.0);
        assert_eq!(is_point_2_on_line_segment_2(&p3, &s), true);

        let p3 = Point2::new(5.0, 5.1);
        assert_eq!(is_point_2_on_line_segment_2(&p3, &s), false);

        let p3 = Point2::new(f64::sqrt(3.0), f64::sqrt(3.0));
        assert_eq!(is_point_2_on_line_segment_2(&p3, &s), true);
    }

    #[test]
    fn test_locate_point_2_line_segment_2() {
        let p1 = Point2::new(0.0, 0.0);
        let p2 = Point2::new(10.0, 10.0);
        let s = LineSegment2::new(p1, p2);

        let p3 = Point2::new(5.0, 5.0);
        assert_eq!(
            locate_point_2_line_segment_2(&p3, &s),
            Point2Segment2Location::On
        );

        let p3 = Point2::new(5.0, 5.1);
        assert_eq!(
            locate_point_2_line_segment_2(&p3, &s),
            Point2Segment2Location::Left
        );

        let p3 = Point2::new(5.0, 4.9);
        assert_eq!(
            locate_point_2_line_segment_2(&p3, &s),
            Point2Segment2Location::Right
        );

        let p3 = Point2::new(0.0, 0.0);
        assert_eq!(
            locate_point_2_line_segment_2(&p3, &s),
            Point2Segment2Location::On
        );

        let p3 = Point2::new(10.0, 10.0);
        assert_eq!(
            locate_point_2_line_segment_2(&p3, &s),
            Point2Segment2Location::On
        );

        let p3 = Point2::new(10.0, 10.1);
        assert_eq!(
            locate_point_2_line_segment_2(&p3, &s),
            Point2Segment2Location::Left
        );

        let p3 = Point2::new(10.0, 9.9);
        assert_eq!(
            locate_point_2_line_segment_2(&p3, &s),
            Point2Segment2Location::Right
        );

        let p3 = Point2::new(-10.0, -10.0);
        assert_eq!(
            locate_point_2_line_segment_2(&p3, &s),
            Point2Segment2Location::Collinear
        );
    }
}