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

use alga::general::{Field, JoinSemilattice, Lattice, MeetSemilattice, Real};
use alga::linear::{AffineSpace, EuclideanSpace};

use base::allocator::Allocator;
use base::dimension::DimName;
use base::{DefaultAllocator, Scalar, VectorN};

use geometry::Point;

impl<N: Scalar + Field, D: DimName> AffineSpace for Point<N, D>
where
    N: Scalar + Field,
    DefaultAllocator: Allocator<N, D>,
{
    type Translation = VectorN<N, D>;
}

impl<N: Real, D: DimName> EuclideanSpace for Point<N, D>
where DefaultAllocator: Allocator<N, D>
{
    type Coordinates = VectorN<N, D>;
    type Real = N;

    #[inline]
    fn origin() -> Self {
        Self::origin()
    }

    #[inline]
    fn coordinates(&self) -> Self::Coordinates {
        self.coords.clone()
    }

    #[inline]
    fn from_coordinates(coords: Self::Coordinates) -> Self {
        Self::from(coords)
    }

    #[inline]
    fn scale_by(&self, n: N) -> Self {
        self * n
    }
}

/*
 *
 * Ordering
 *
 */
impl<N, D: DimName> MeetSemilattice for Point<N, D>
where
    N: Scalar + MeetSemilattice,
    DefaultAllocator: Allocator<N, D>,
{
    #[inline]
    fn meet(&self, other: &Self) -> Self {
        Point::from(self.coords.meet(&other.coords))
    }
}

impl<N, D: DimName> JoinSemilattice for Point<N, D>
where
    N: Scalar + JoinSemilattice,
    DefaultAllocator: Allocator<N, D>,
{
    #[inline]
    fn join(&self, other: &Self) -> Self {
        Point::from(self.coords.join(&other.coords))
    }
}

impl<N, D: DimName> Lattice for Point<N, D>
where
    N: Scalar + Lattice,
    DefaultAllocator: Allocator<N, D>,
{
    #[inline]
    fn meet_join(&self, other: &Self) -> (Self, Self) {
        let (meet, join) = self.coords.meet_join(&other.coords);

        (Point::from(meet), Point::from(join))
    }
}