use crate::util::*;
use approx::*;
use blas_array2::blas1::nrm2::NRM2;
use blas_array2::util::*;

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

    #[test]
    fn test_example() {
        for incx in [1, 2] {
            let n = 100;
            let x = random_array::<f64>(1000);
            let x_slc = slice_1d(n, incx);
            let out = NRM2::default().x(x.slice(x_slc)).run().unwrap();
            let expected = f64::sqrt(x.slice(x_slc).view().mapv(|x| x * x).sum());
            assert_relative_eq!(out, expected, epsilon = 1.0e-6);
        }

        for incx in [1, 2] {
            let n = 100;
            let x = random_array::<c64>(1000);
            let x_slc = slice_1d(n, incx);
            let out = NRM2::default().x(x.slice(x_slc)).run().unwrap();
            let expected = f64::sqrt(x.slice(x_slc).view().mapv(|x| (x * x.conj()).re).sum());
            assert_relative_eq!(out, expected, epsilon = 1.0e-6);
        }
    }
}