argmin-math 0.5.1

argmin math abstractions
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 crate::ArgminMinMax;
use faer::{mat, unzip, zip, Mat, MatRef};

impl<E: PartialOrd + Copy> ArgminMinMax for Mat<E> {
    #[inline]
    fn max(a: &Self, b: &Self) -> Self {
        faer::zip!(a, b).map(|faer::unzip!(a, b)| {
            let aa = *a;
            let bb = *b;
            //@note(geo-ant) directly cribbed from the nalgebra implementation
            // will have the same problems with NaN values.
            if aa > bb {
                aa
            } else {
                bb
            }
        })
    }

    #[inline]
    fn min(a: &Mat<E>, b: &Mat<E>) -> Mat<E> {
        faer::zip!(a, b).map(|faer::unzip!(a, b)| {
            let aa = *a;
            let bb = *b;
            //@note(geo-ant) directly cribbed from the nalgebra implementation
            // will have the same problems with NaN values.
            if aa < bb {
                aa
            } else {
                bb
            }
        })
    }
}