opensrdk-linear-algebra 0.8.5

Standard linear algebra library using blas and lapack for OpenSRDK toolchain.
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
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51

use crate::{number::Number, Matrix};
use rayon::prelude::*;

impl<T> Matrix<T>
where
    T: Number,
{
    /// # Linear product
    pub fn linear_prod(&self, rhs: &Matrix<T>) -> T {
        if !self.is_same_size(rhs) {
            panic!("Dimension mismatch.")
        }

        self.elems
            .par_iter()
            .zip(rhs.elems.par_iter())
            .map(|(&s, &r)| s * r)
            .sum()
    }

    /// # Hamadard product
    pub fn hadamard_prod(self, rhs: &Matrix<T>) -> Matrix<T> {
        if !self.is_same_size(rhs) {
            panic!("Dimension mismatch.")
        }
        let mut slf = self;

        slf.elems
            .par_iter_mut()
            .zip(rhs.elems.par_iter())
            .map(|(s, &r)| {
                *s *= r;
            })
            .collect::<Vec<_>>();

        slf
    }
}

#[cfg(test)]
mod tests {
    use crate::*;
    #[test]
    fn it_works() {
        let a = DiagonalMatrix::<f64>::identity(2).mat();
        let b = DiagonalMatrix::<f64>::identity(2).mat();
        let c = a.hadamard_prod(&b);

        assert_eq!(b[0][0], c[0][0])
    }
}