mtl 0.1.5

Matrix template library. Dense2D matrix, SparseMatrix.
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
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111

use std::fmt;
use std::ops;

pub struct MatrixElementWiseView<'a, T: 'a> {
    pub inner: &'a [T],
    pub dim: (usize, usize)
}


// sub matrix view
pub struct MatrixSubView<'a, T: 'a> {
    pub inner: &'a [T],
    pub orig_dim: (usize, usize),
    pub offset: (usize, usize),
    pub dim: (usize, usize)
}
pub struct MatrixMutSubView<'a, T: 'a> {
    pub inner: &'a mut [T],
    pub orig_dim: (usize, usize),
    pub offset: (usize, usize),
    pub dim: (usize, usize)
}


impl<'a, T> MatrixSubView<'a, T> {
    #[inline]
    fn map_index(&self, row: usize, col: usize) -> Option<usize> {
        if row >= self.dim.0 {
            None
        } else if col >= self.dim.1 {
            None
        } else {
            Some((self.offset.0 + row) * self.orig_dim.1 + self.offset.1 + col)
        }
    }
}

impl<'a, T> ops::Index<(usize, usize)> for MatrixSubView<'a, T> {
    type Output = T;

    #[inline]
    fn index(&self, (row, col): (usize, usize)) -> &T {
        match self.map_index(row, col) {
            Some(i) => &self.inner[i],
            None    => panic!("index over range")
        }
    }
}


impl<'a, T> ops::Index<usize> for MatrixSubView<'a, T> {
    type Output = [T];

    #[inline]
    fn index<'b>(&'b self, row: usize) -> &'b [T] {
        assert!(row < self.dim.0);
        let start = self.map_index(row, 0).unwrap();
        let end  = self.map_index(row, self.dim.1 - 1).unwrap();
        &self.inner[start .. end + 1]
    }
}

// impl<'a, T> ops::IndexMut<usize> for MatrixSubView<'a, T> {
//     #[inline]
//     fn index_mut<'a>(&'a mut self, row: usize) -> &'a mut [T] {
//         assert!(row < self.dim.0);
//         &mut self.data[
//     }
// }



impl<'a, T: fmt::Debug + Copy> fmt::Debug for MatrixSubView<'a, T> {
    fn fmt(&self, f: &mut fmt::Formatter) -> fmt::Result {
        let mut view = Vec::new();
        for j in 0 .. self.dim.0 {
            for i in 0 .. self.dim.1 {
                view.push(self[j][i]);
            }
        }
        try!(write!(f, "<MatrixSubView dim={:?}, {:?}>", self.dim, view));
        Ok(())
    }
}

impl<'a, T: fmt::Display> fmt::Display for MatrixSubView<'a, T> {
    fn fmt(&self, f: &mut fmt::Formatter) -> fmt::Result {
        try!(writeln!(f, ""));
        for j in 0 .. self.dim.0 {
            if j == 0 {
                try!(write!(f, "[["));
            } else {
                try!(write!(f, " ["));
            }
            for i in 0 .. self.dim.1 {
                try!(write!(f, "{:-4.}", self[j][i]));
                if i == self.dim.1 - 1 {
                    try!(write!(f, "]"));
                } else {
                    try!(write!(f, ", "));
                }
            }
            if j == self.dim.0 - 1 {
                try!(write!(f, "]"));
            } else {
                try!(writeln!(f, ""));
            }
        }
        Ok(())
    }
}