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
use super::{Point, PointLocation};
use crate::{Error, Orientation, PolygonScalar, PolygonScalarRef};
use claim::debug_assert_ok;
use std::ops::*;
pub struct Triangle<T>([Point<T, 2>; 3]);
impl<T> Triangle<T>
where
T: PolygonScalar,
for<'a> &'a T: PolygonScalarRef<&'a T, T>,
{
pub fn new(pts: [Point<T, 2>; 3]) -> Triangle<T> {
let triangle = Triangle(pts);
debug_assert_ok!(triangle.validate());
triangle
}
pub fn validate(&self) -> Result<(), Error> {
self.view().validate()
}
pub fn locate(&self, pt: &Point<T, 2>) -> PointLocation {
self.view().locate(pt)
}
pub fn view(&'_ self) -> TriangleView<'_, T> {
TriangleView([&self.0[0], &self.0[1], &self.0[2]])
}
}
pub struct TriangleView<'a, T>([&'a Point<T, 2>; 3]);
impl<'a, T> TriangleView<'a, T>
where
T: PolygonScalar,
{
pub fn new(pts: [&'a Point<T, 2>; 3]) -> TriangleView<'a, T> {
let triangle = TriangleView(pts);
debug_assert_ok!(triangle.validate());
triangle
}
pub fn validate(&self) -> Result<(), Error> {
let arr = &self.0;
if arr.index(0).orientation(&arr[1], &arr[2]) != Orientation::CounterClockWise {
Err(Error::ClockWiseViolation)
} else {
Ok(())
}
}
pub fn locate(&self, pt: &Point<T, 2>) -> PointLocation {
use Orientation::*;
debug_assert_ok!(self.validate());
let [a, b, c] = self.0;
let ab = a.orientation(b, pt);
let bc = b.orientation(c, pt);
let ca = c.orientation(a, pt);
if ab == ClockWise || bc == ClockWise || ca == ClockWise {
PointLocation::Outside
} else if ab == CoLinear || bc == CoLinear || ca == CoLinear {
PointLocation::OnBoundary
} else {
PointLocation::Inside
}
}
}