ommx 2.5.1

Open Mathematical prograMming eXchange (OMMX)
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

use ::approx::AbsDiffEq;

use super::*;

fn compare<T>(a: &Sampled<T>, b: &Sampled<T>, mut f: impl FnMut(&T, &T) -> bool) -> bool {
    if a.offsets.len() != b.offsets.len() {
        return false;
    }
    for (id, offset) in a.offsets.iter() {
        debug_assert!(*offset < a.data.len());
        let Some(other_offset) = b.offsets.get(id) else {
            return false;
        };
        debug_assert!(*other_offset < b.data.len());
        if !f(&a.data[*offset], &b.data[*other_offset]) {
            return false;
        }
    }
    true
}

impl<T> PartialEq for Sampled<T>
where
    T: PartialEq,
{
    fn eq(&self, other: &Self) -> bool {
        compare(self, other, |a, b| a == b)
    }
}

impl<T> AbsDiffEq for Sampled<T>
where
    T: AbsDiffEq,
    <T as AbsDiffEq>::Epsilon: Clone,
{
    type Epsilon = <T as AbsDiffEq>::Epsilon;

    fn default_epsilon() -> Self::Epsilon {
        T::default_epsilon()
    }

    fn abs_diff_eq(&self, other: &Self, epsilon: Self::Epsilon) -> bool {
        compare(self, other, |a, b| a.abs_diff_eq(b, epsilon.clone()))
    }
}

#[cfg(test)]
mod tests {
    use super::*;

    #[test]
    fn test_partial_eq() {
        // Different as struct, but mathematical equal
        let a = Sampled::new(
            [[SampleID(1), SampleID(2)], [SampleID(5), SampleID(7)]],
            [1, 2],
        )
        .unwrap();
        let b = Sampled::new(
            [[SampleID(7), SampleID(5)], [SampleID(1), SampleID(2)]],
            [2, 1],
        )
        .unwrap();
        assert_eq!(a, b);
    }
}