opensrdk-linear-algebra 0.9.1

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
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87

use std::ops::{Index, IndexMut};

use crate::{sparse::SparseTensor, Number, TensorError};

impl<T> Index<&[usize]> for SparseTensor<T>
where
    T: Number,
{
    type Output = T;

    fn index(&self, index: &[usize]) -> &Self::Output {
        if index.len() != self.sizes.len() {
            panic!("{}", TensorError::RankMismatch);
        }
        for (rank, &d) in index.iter().enumerate() {
            if self.sizes[rank] <= d {
                panic!("{}", TensorError::OutOfRange);
            }
        }

        self.elems.get(index).unwrap_or(&self.default)
    }
}

impl<T> IndexMut<&[usize]> for SparseTensor<T>
where
    T: Number,
{
    fn index_mut(&mut self, index: &[usize]) -> &mut Self::Output {
        if index.len() != self.sizes.len() {
            panic!("{}", TensorError::RankMismatch);
        }
        for (rank, &d) in index.iter().enumerate() {
            if self.sizes[rank] <= d {
                panic!("{}", TensorError::OutOfRange);
            }
        }

        self.elems.entry(index.to_vec()).or_default()
    }
}

#[cfg(test)]
mod tests {
    use crate::{sparse::SparseTensor, *};

    #[test]
    fn index() {
        let mut tensor = SparseTensor::new(vec![2, 3]);
        tensor[&[0, 0]] = 1.0;
        tensor[&[1, 1]] = 2.0;
        tensor[&[1, 2]] = 3.0;

        assert_eq!(tensor[&[0, 0]], 1.0);
        assert_eq!(tensor[&[0, 1]], 0.0);
        assert_eq!(tensor[&[0, 2]], 0.0);
        assert_eq!(tensor[&[1, 0]], 0.0);
        assert_eq!(tensor[&[1, 1]], 2.0);
        assert_eq!(tensor[&[1, 2]], 3.0);
    }

    #[test]
    fn index_mut() {
        let mut tensor = SparseTensor::new(vec![2, 3]);
        tensor[&[0, 0]] = 1.0;
        tensor[&[1, 1]] = 2.0;
        tensor[&[1, 2]] = 3.0;

        assert_eq!(tensor[&[0, 0]], 1.0);
        assert_eq!(tensor[&[0, 1]], 0.0);
        assert_eq!(tensor[&[0, 2]], 0.0);
        assert_eq!(tensor[&[1, 0]], 0.0);
        assert_eq!(tensor[&[1, 1]], 2.0);
        assert_eq!(tensor[&[1, 2]], 3.0);

        tensor[&[0, 0]] = 0.0;
        tensor[&[1, 1]] = 0.0;
        tensor[&[1, 2]] = 0.0;

        assert_eq!(tensor[&[0, 0]], 0.0);
        assert_eq!(tensor[&[0, 1]], 0.0);
        assert_eq!(tensor[&[0, 2]], 0.0);
        assert_eq!(tensor[&[1, 0]], 0.0);
        assert_eq!(tensor[&[1, 1]], 0.0);
        assert_eq!(tensor[&[1, 2]], 0.0);
    }
}