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
use std::fmt; use traits::{MatrixRawGet, MatrixRawSet, MatrixShape, MatrixColumnAccess, SameShape}; use matrix::write_mat; impl<T: MatrixShape> MatrixColumnAccess for T { unsafe fn unsafe_col(self, c: usize) -> ColumnAccessor<T> { ColumnAccessor::unsafe_new(self, c) } fn col(self, c: usize) -> ColumnAccessor<T> { ColumnAccessor::new(self, c) } } #[derive(Copy)] pub struct ColumnAccessor<T> { base: T, col: usize } impl<T: MatrixShape> ColumnAccessor<T> { pub unsafe fn unsafe_new(base: T, col: usize) -> ColumnAccessor<T> { ColumnAccessor{ base: base, col: col } } pub fn new(base: T, col: usize) -> ColumnAccessor<T> { assert!(col < base.ncol()); ColumnAccessor{ base: base, col: col } } } impl<T: MatrixShape> MatrixShape for ColumnAccessor<T> { fn nrow(&self) -> usize { self.base.nrow() } fn ncol(&self) -> usize { 1 } } impl<T: MatrixShape> SameShape for ColumnAccessor<T> { fn same_shape(&self, nrow: usize, ncol: usize) -> bool { self.nrow() == nrow && self.ncol() == ncol } } impl<T: MatrixRawGet + MatrixShape> MatrixRawGet for ColumnAccessor<T> { unsafe fn raw_get(&self, r: usize, _: usize) -> f64 { self.base.raw_get(r, self.col) } } impl<T: MatrixRawSet + MatrixShape> MatrixRawSet for ColumnAccessor<T> { unsafe fn raw_set(&self, r: usize, _: usize, v: f64) { self.base.raw_set(r, self.col, v) } } impl<T: MatrixShape + MatrixRawGet> fmt::Display for ColumnAccessor<T> { fn fmt(&self, buf: &mut fmt::Formatter) -> fmt::Result { write_mat(buf, self) } } impl<T: Clone> Clone for ColumnAccessor<T> { fn clone(&self) -> ColumnAccessor<T> { ColumnAccessor{ base: self.base.clone(), col: self.col } } }