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

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

use crate::base::{SVector, Scalar};

use crate::geometry::Point;

impl<T: Scalar + Field, const D: usize> AffineSpace for Point<T, D>
where
    T: Scalar + Field,
{
    type Translation = SVector<T, D>;
}

impl<T: RealField + simba::scalar::RealField, const D: usize> EuclideanSpace for Point<T, D> {
    type Coordinates = SVector<T, D>;
    type RealField = T;

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

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

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

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

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

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

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

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