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

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

use super::location_enum::Point2ArcSegment2Location;

pub fn is_point_2_on_arc_segment_2<T: NumberType>(
    point: &Point2<T>,
    arc_segment: &impl Segment2<T>,
) -> bool {
    let location = locate_point_2_arc_segment_2(point, arc_segment);
    match location {
        Point2ArcSegment2Location::On => true,
        _ => false,
    }
}

pub fn locate_point_2_arc_segment_2<T: NumberType>(
    point: &Point2<T>,
    arc_segment: &impl Segment2<T>,
) -> Point2ArcSegment2Location {
    let center = arc_segment.center();
    let radius = arc_segment.radius();
    let distance = center.distance(point);
    if distance.equals(radius) {
        let vector = *point - center;
        let source = arc_segment.source() - center;
        let target = arc_segment.target() - center;
        let mut radian = source.radian_to(&target);
        if source == target {
            radian = T::pi() * T::from_f64(2.0);
        }
        let radian_source = vector.radian_to(&source);
        let radian_target = vector.radian_to(&target);
        let source_radian = source.radian_to(&vector);
        let target_radian = target.radian_to(&vector);
        if radian_source.equals(T::zero())
            || radian_target.equals(T::zero())
            || source_radian.equals(T::zero())
            || target_radian.equals(T::zero())
        {
            Point2ArcSegment2Location::On
        } else if (source_radian + radian_target).equals(radian) {
            Point2ArcSegment2Location::On
        } else {
            Point2ArcSegment2Location::NotOn
        }
    } else {
        Point2ArcSegment2Location::NotOn
    }
}