Chris@49
|
1 // Copyright (C) 2008-2013 NICTA (www.nicta.com.au)
|
Chris@49
|
2 // Copyright (C) 2008-2013 Conrad Sanderson
|
Chris@49
|
3 //
|
Chris@49
|
4 // This Source Code Form is subject to the terms of the Mozilla Public
|
Chris@49
|
5 // License, v. 2.0. If a copy of the MPL was not distributed with this
|
Chris@49
|
6 // file, You can obtain one at http://mozilla.org/MPL/2.0/.
|
Chris@49
|
7
|
Chris@49
|
8
|
Chris@49
|
9 //! \addtogroup subview_cube
|
Chris@49
|
10 //! @{
|
Chris@49
|
11
|
Chris@49
|
12
|
Chris@49
|
13 //! Class for storing data required to construct or apply operations to a subcube
|
Chris@49
|
14 //! (i.e. where the subcube starts and ends as well as a reference/pointer to the original cube),
|
Chris@49
|
15 template<typename eT>
|
Chris@49
|
16 class subview_cube : public BaseCube<eT, subview_cube<eT> >
|
Chris@49
|
17 {
|
Chris@49
|
18 public:
|
Chris@49
|
19
|
Chris@49
|
20 typedef eT elem_type;
|
Chris@49
|
21 typedef typename get_pod_type<elem_type>::result pod_type;
|
Chris@49
|
22
|
Chris@49
|
23 arma_aligned const Cube<eT>& m;
|
Chris@49
|
24
|
Chris@49
|
25 const uword aux_row1;
|
Chris@49
|
26 const uword aux_col1;
|
Chris@49
|
27 const uword aux_slice1;
|
Chris@49
|
28
|
Chris@49
|
29 const uword n_rows;
|
Chris@49
|
30 const uword n_cols;
|
Chris@49
|
31 const uword n_elem_slice;
|
Chris@49
|
32 const uword n_slices;
|
Chris@49
|
33 const uword n_elem;
|
Chris@49
|
34
|
Chris@49
|
35
|
Chris@49
|
36 protected:
|
Chris@49
|
37
|
Chris@49
|
38 arma_inline subview_cube(const Cube<eT>& in_m, const uword in_row1, const uword in_col1, const uword in_slice1, const uword in_n_rows, const uword in_n_cols, const uword in_n_slices);
|
Chris@49
|
39
|
Chris@49
|
40
|
Chris@49
|
41 public:
|
Chris@49
|
42
|
Chris@49
|
43 inline ~subview_cube();
|
Chris@49
|
44
|
Chris@49
|
45 inline void operator+= (const eT val);
|
Chris@49
|
46 inline void operator-= (const eT val);
|
Chris@49
|
47 inline void operator*= (const eT val);
|
Chris@49
|
48 inline void operator/= (const eT val);
|
Chris@49
|
49
|
Chris@49
|
50 // deliberately returning void
|
Chris@49
|
51 template<typename T1> inline void operator= (const BaseCube<eT,T1>& x);
|
Chris@49
|
52 template<typename T1> inline void operator+= (const BaseCube<eT,T1>& x);
|
Chris@49
|
53 template<typename T1> inline void operator-= (const BaseCube<eT,T1>& x);
|
Chris@49
|
54 template<typename T1> inline void operator%= (const BaseCube<eT,T1>& x);
|
Chris@49
|
55 template<typename T1> inline void operator/= (const BaseCube<eT,T1>& x);
|
Chris@49
|
56
|
Chris@49
|
57 inline void operator= (const subview_cube& x);
|
Chris@49
|
58 inline void operator+= (const subview_cube& x);
|
Chris@49
|
59 inline void operator-= (const subview_cube& x);
|
Chris@49
|
60 inline void operator%= (const subview_cube& x);
|
Chris@49
|
61 inline void operator/= (const subview_cube& x);
|
Chris@49
|
62
|
Chris@49
|
63 template<typename T1> inline void operator= (const Base<eT,T1>& x);
|
Chris@49
|
64 template<typename T1> inline void operator+= (const Base<eT,T1>& x);
|
Chris@49
|
65 template<typename T1> inline void operator-= (const Base<eT,T1>& x);
|
Chris@49
|
66 template<typename T1> inline void operator%= (const Base<eT,T1>& x);
|
Chris@49
|
67 template<typename T1> inline void operator/= (const Base<eT,T1>& x);
|
Chris@49
|
68
|
Chris@49
|
69 inline static void extract(Cube<eT>& out, const subview_cube& in);
|
Chris@49
|
70 inline static void plus_inplace(Cube<eT>& out, const subview_cube& in);
|
Chris@49
|
71 inline static void minus_inplace(Cube<eT>& out, const subview_cube& in);
|
Chris@49
|
72 inline static void schur_inplace(Cube<eT>& out, const subview_cube& in);
|
Chris@49
|
73 inline static void div_inplace(Cube<eT>& out, const subview_cube& in);
|
Chris@49
|
74
|
Chris@49
|
75 inline static void extract(Mat<eT>& out, const subview_cube& in);
|
Chris@49
|
76 inline static void plus_inplace(Mat<eT>& out, const subview_cube& in);
|
Chris@49
|
77 inline static void minus_inplace(Mat<eT>& out, const subview_cube& in);
|
Chris@49
|
78 inline static void schur_inplace(Mat<eT>& out, const subview_cube& in);
|
Chris@49
|
79 inline static void div_inplace(Mat<eT>& out, const subview_cube& in);
|
Chris@49
|
80
|
Chris@49
|
81 template<typename functor> inline void transform(functor F);
|
Chris@49
|
82 template<typename functor> inline void imbue(functor F);
|
Chris@49
|
83
|
Chris@49
|
84 inline void fill(const eT val);
|
Chris@49
|
85 inline void zeros();
|
Chris@49
|
86 inline void ones();
|
Chris@49
|
87
|
Chris@49
|
88 inline eT at_alt (const uword i) const;
|
Chris@49
|
89
|
Chris@49
|
90 inline eT& operator[](const uword i);
|
Chris@49
|
91 inline eT operator[](const uword i) const;
|
Chris@49
|
92
|
Chris@49
|
93 inline eT& operator()(const uword i);
|
Chris@49
|
94 inline eT operator()(const uword i) const;
|
Chris@49
|
95
|
Chris@49
|
96 arma_inline eT& operator()(const uword in_row, const uword in_col, const uword in_slice);
|
Chris@49
|
97 arma_inline eT operator()(const uword in_row, const uword in_col, const uword in_slice) const;
|
Chris@49
|
98
|
Chris@49
|
99 arma_inline eT& at(const uword in_row, const uword in_col, const uword in_slice);
|
Chris@49
|
100 arma_inline eT at(const uword in_row, const uword in_col, const uword in_slice) const;
|
Chris@49
|
101
|
Chris@49
|
102 arma_inline eT* slice_colptr(const uword in_slice, const uword in_col);
|
Chris@49
|
103 arma_inline const eT* slice_colptr(const uword in_slice, const uword in_col) const;
|
Chris@49
|
104
|
Chris@49
|
105 inline bool check_overlap(const subview_cube& x) const;
|
Chris@49
|
106 inline bool check_overlap(const Mat<eT>& x) const;
|
Chris@49
|
107
|
Chris@49
|
108
|
Chris@49
|
109 private:
|
Chris@49
|
110
|
Chris@49
|
111 friend class Mat<eT>;
|
Chris@49
|
112 friend class Cube<eT>;
|
Chris@49
|
113
|
Chris@49
|
114 subview_cube();
|
Chris@49
|
115 };
|
Chris@49
|
116
|
Chris@49
|
117
|
Chris@49
|
118
|
Chris@49
|
119 //! @}
|