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

use crate::algebra::abstr::AbsDiffEq;
use crate::algebra::abstr::Field;
use crate::algebra::abstr::Scalar;
use crate::algebra::linear::matrix::General;

impl<T> AbsDiffEq for General<T>
where
    T: Field + Scalar + AbsDiffEq<Epsilon = T>,
{
    type Epsilon = T;

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

    fn abs_diff_eq(&self, other: &General<T>, epsilon: T) -> bool {
        if self.dim() != other.dim() {
            return false;
        }

        for (a, b) in self.iter().zip(other.iter()) {
            if a.abs_diff_ne(b, epsilon) {
                return false;
            }
        }

        true
    }
}