argmin-math 0.5.1

use crate::ArgminRandom;
use faer::{unzipped, Entity, Mat};
use rand::distr::uniform::SampleUniform;

impl<E: Entity + PartialOrd + SampleUniform> ArgminRandom for Mat<E> {
    fn rand_from_range<R: rand::Rng>(min: &Self, max: &Self, rng: &mut R) -> Self {
        assert!(
            min.nrows() != 0 || min.ncols() != 0,
            "internal error: empty matrix unexpected"
        );
        assert_eq!(
            min.shape(),
            max.shape(),
            "internal error: matrices of same shape expected"
        );

        faer::zipped_rw!(min, max).map(|unzipped!(min, max)| {
            let a = min.read();
            let b = max.read();
            //@note(geo-ant) code was copied from the nalgebra implementation
            // to get the exact same behaviour
            // Do not require a < b:
            // We do want to know if a and b are *exactly* the same.
            #[allow(clippy::float_cmp)]
            if a == b {
                a
            } else if a < b {
                rng.random_range(a..b)
            } else {
                rng.random_range(b..a)
            }
        })
    }
}