rinia 0.0.20

Stable numeric foundations for graphics, animation, and simulation.
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
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
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102

use crate::algebra::{ApproxEqAbs, ApproxEqRel};

// Blanket ApproxEqAbs for Arrays
impl<T, U, const N: usize> ApproxEqAbs<[U; N]> for [T; N]
where
    T: ApproxEqAbs<U>,
{
    type Tolerance = T::Tolerance;

    const DEFAULT_TOLERANCE_ABS: Self::Tolerance = T::DEFAULT_TOLERANCE_ABS;

    #[inline]
    fn approx_eq_abs_tol(&self, other: &[U; N], tol: Self::Tolerance) -> bool {
        self.iter().zip(other.iter()).all(|(lhs, rhs)| lhs.approx_eq_abs_tol(rhs, tol))
    }
}

// Blanket ApproxEqRel for Arrays
impl<T, U, const N: usize> ApproxEqRel<[U; N]> for [T; N]
where
    T: ApproxEqRel<U>,
{
    type Tolerance = T::Tolerance;

    const DEFAULT_TOLERANCE_REL: Self::Tolerance = T::DEFAULT_TOLERANCE_REL;

    #[inline]
    fn approx_eq_rel_tol(&self, other: &[U; N], tol: Self::Tolerance) -> bool {
        self.iter().zip(other.iter()).all(|(lhs, rhs)| lhs.approx_eq_rel_tol(rhs, tol))
    }
}

// Blanket ApproxEqAbs for Slices
impl<T, U> ApproxEqAbs<[U]> for [T]
where
    T: ApproxEqAbs<U>,
{
    type Tolerance = T::Tolerance;

    const DEFAULT_TOLERANCE_ABS: Self::Tolerance = T::DEFAULT_TOLERANCE_ABS;

    #[inline]
    fn approx_eq_abs_tol(&self, other: &[U], tol: Self::Tolerance) -> bool {
        self.len() == other.len()
            && self.iter().zip(other.iter()).all(|(lhs, rhs)| lhs.approx_eq_abs_tol(rhs, tol))
    }
}

// Blanket ApproxEqRel for Slices
impl<T, U> ApproxEqRel<[U]> for [T]
where
    T: ApproxEqRel<U>,
{
    type Tolerance = T::Tolerance;

    const DEFAULT_TOLERANCE_REL: Self::Tolerance = T::DEFAULT_TOLERANCE_REL;

    #[inline]
    fn approx_eq_rel_tol(&self, other: &[U], tol: Self::Tolerance) -> bool {
        self.len() == other.len()
            && self.iter().zip(other.iter()).all(|(lhs, rhs)| lhs.approx_eq_rel_tol(rhs, tol))
    }
}

#[cfg(any(feature = "alloc", feature = "std", test))]
mod vec_equality {
    use alloc::vec::Vec;

    use crate::algebra::{ApproxEqAbs, ApproxEqRel};

    // Blanket ApproxEqAbs for Vec
    impl<T, U> ApproxEqAbs<Vec<U>> for Vec<T>
    where
        T: ApproxEqAbs<U>,
    {
        type Tolerance = T::Tolerance;

        const DEFAULT_TOLERANCE_ABS: Self::Tolerance = T::DEFAULT_TOLERANCE_ABS;

        #[inline]
        fn approx_eq_abs_tol(&self, other: &Vec<U>, tol: Self::Tolerance) -> bool {
            self.len() == other.len()
                && self.iter().zip(other.iter()).all(|(lhs, rhs)| lhs.approx_eq_abs_tol(rhs, tol))
        }
    }

    // Blanket ApproxEqRel for Vec
    impl<T, U> ApproxEqRel<Vec<U>> for Vec<T>
    where
        T: ApproxEqRel<U>,
    {
        type Tolerance = T::Tolerance;

        const DEFAULT_TOLERANCE_REL: Self::Tolerance = T::DEFAULT_TOLERANCE_REL;

        #[inline]
        fn approx_eq_rel_tol(&self, other: &Vec<U>, tol: Self::Tolerance) -> bool {
            self.len() == other.len()
                && self.iter().zip(other.iter()).all(|(lhs, rhs)| lhs.approx_eq_rel_tol(rhs, tol))
        }
    }
}