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

use format::{Conventional, Diagonal};
use operation::{MultiplySelf, Transpose};
use {Element, Number};

#[cfg(feature = "acceleration")]
mod acceleration;

impl<T> MultiplySelf<Diagonal<T>> for Conventional<T>
where
    T: Element + Number,
{
    #[inline]
    fn multiply_self(&mut self, right: &Diagonal<T>) {
        let (rows, insides, columns) = (self.rows, self.columns, right.columns);
        debug_assert_eq!(insides, right.rows);
        self.resize((rows, columns));
        for j in 0..insides {
            let factor = right[j];
            for i in 0..rows {
                self[(i, j)] = factor * self[(i, j)];
            }
        }
    }
}

impl<T: Element> Transpose for Conventional<T> {
    fn transpose(&self) -> Self {
        let (rows, columns) = (self.rows, self.columns);
        let mut matrix = Conventional::new((columns, rows));
        for i in 0..rows {
            for j in 0..columns {
                matrix.values[i * columns + j] = self.values[j * rows + i];
            }
        }
        matrix
    }
}

#[cfg(test)]
mod tests {
    use prelude::*;

    #[test]
    fn multiply_self() {
        let mut matrix = Conventional::from_vec(
            (3, 2),
            matrix![
                1.0, 4.0;
                2.0, 5.0;
                3.0, 6.0;
            ],
        );
        let right = Diagonal::from_vec((2, 4), vec![2.0, 3.0]);
        matrix.multiply_self(&right);
        assert_eq!(
            &*matrix,
            &*matrix![
                2.0, 12.0, 0.0, 0.0;
                4.0, 15.0, 0.0, 0.0;
                6.0, 18.0, 0.0, 0.0;
            ]
        );
    }

    #[test]
    fn transpose() {
        let matrix = Conventional::from_vec(
            (3, 2),
            matrix![
                1.0, 4.0;
                2.0, 5.0;
                3.0, 6.0;
            ],
        );
        assert_eq!(
            matrix.transpose(),
            Conventional::from_vec(
                (2, 3),
                matrix![
                    1.0, 2.0, 3.0;
                    4.0, 5.0, 6.0;
                ],
            )
        );
    }
}