ommx 3.0.0-alpha.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

use super::*;
use ::approx::AbsDiffEq;

/// Compare two polynomial by maximum coefficient difference.
impl<M: Monomial> AbsDiffEq for PolynomialBase<M> {
    type Epsilon = crate::ATol;
    fn default_epsilon() -> Self::Epsilon {
        crate::ATol::default()
    }
    fn abs_diff_eq(&self, other: &Self, epsilon: Self::Epsilon) -> bool {
        let residual = self - other;
        residual
            .max_coefficient_abs()
            .is_none_or(|max| max <= *epsilon)
    }
}

#[cfg(test)]
mod tests {
    use super::*;
    use crate::linear;
    use ::approx::assert_abs_diff_eq;

    #[test]
    fn test_abs_diff_eq() {
        let zero = PolynomialBase::default();
        let small = PolynomialBase {
            terms: [(linear!(1), 1e-11.try_into().unwrap())]
                .into_iter()
                .collect(),
        };
        assert_abs_diff_eq!(small, zero);
    }
}