blas-array2 0.3.0

Parameter-optional BLAS wrapper by ndarray::Array (Ix1 or Ix2).
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

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);
        }
    }
}