rene 0.2.0

Computational geometry.
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

use crate::geometries::Point;
use crate::locatable::{Locatable, Location};
use crate::operations::{is_point_in_segment, Orient};
use crate::traits::Elemental;

use super::types::Contour;

impl<Scalar: PartialOrd> Locatable<&Point<Scalar>> for &Contour<Scalar>
where
    Point<Scalar>: PartialEq,
    for<'a> &'a Point<Scalar>: Elemental<Coordinate = &'a Scalar> + Orient,
{
    fn locate(self, point: &Point<Scalar>) -> Location {
        for index in 0..self.vertices.len() - 1 {
            if is_point_in_segment(
                point,
                &self.vertices[index],
                &self.vertices[index + 1],
            ) {
                return Location::Boundary;
            }
        }
        if is_point_in_segment(
            point,
            &self.vertices[self.vertices.len() - 1],
            &self.vertices[0],
        ) {
            Location::Boundary
        } else {
            Location::Exterior
        }
    }
}