Chris@49
|
1 // Copyright (C) 2008-2012 NICTA (www.nicta.com.au)
|
Chris@49
|
2 // Copyright (C) 2008-2012 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 Op
|
Chris@49
|
10 //! @{
|
Chris@49
|
11
|
Chris@49
|
12
|
Chris@49
|
13
|
Chris@49
|
14 //! Class for storing data required for delayed unary operations,
|
Chris@49
|
15 //! such as the operand (e.g. the matrix to which the operation is to be applied) and the unary operator (e.g. inverse).
|
Chris@49
|
16 //! The operand is stored as a reference (which can be optimised away),
|
Chris@49
|
17 //! while the operator is "stored" through the template definition (op_type).
|
Chris@49
|
18 //! The operands can be 'Mat', 'Row', 'Col', 'Op', and 'Glue'.
|
Chris@49
|
19 //! Note that as 'Glue' can be one of the operands, more than one matrix can be stored.
|
Chris@49
|
20 //!
|
Chris@49
|
21 //! For example, we could have:
|
Chris@49
|
22 //! Op< Glue< Mat, Mat, glue_times >, op_htrans >
|
Chris@49
|
23
|
Chris@49
|
24 template<typename T1, typename op_type>
|
Chris@49
|
25 class Op : public Base<typename T1::elem_type, Op<T1, op_type> >
|
Chris@49
|
26 {
|
Chris@49
|
27 public:
|
Chris@49
|
28
|
Chris@49
|
29 typedef typename T1::elem_type elem_type;
|
Chris@49
|
30 typedef typename get_pod_type<elem_type>::result pod_type;
|
Chris@49
|
31
|
Chris@49
|
32 static const bool is_row = ( T1::is_col && (is_same_type<op_type, op_strans>::value || is_same_type<op_type, op_htrans>::value || is_same_type<op_type, op_htrans2>::value) );
|
Chris@49
|
33 static const bool is_col = ( T1::is_row && (is_same_type<op_type, op_strans>::value || is_same_type<op_type, op_htrans>::value || is_same_type<op_type, op_htrans2>::value) ) || (is_same_type<op_type, op_diagvec>::value);
|
Chris@49
|
34
|
Chris@49
|
35 inline explicit Op(const T1& in_m);
|
Chris@49
|
36 inline Op(const T1& in_m, const elem_type in_aux);
|
Chris@49
|
37 inline Op(const T1& in_m, const elem_type in_aux, const uword in_aux_uword_a, const uword in_aux_uword_b);
|
Chris@49
|
38 inline Op(const T1& in_m, const uword in_aux_uword_a, const uword in_aux_uword_b);
|
Chris@49
|
39 inline Op(const T1& in_m, const uword in_aux_uword_a, const uword in_aux_uword_b, const uword in_aux_uword_c, const char junk);
|
Chris@49
|
40 inline ~Op();
|
Chris@49
|
41
|
Chris@49
|
42
|
Chris@49
|
43 arma_aligned const T1& m; //!< storage of reference to the operand (eg. a matrix)
|
Chris@49
|
44 arma_aligned elem_type aux; //!< storage of auxiliary data, user defined format
|
Chris@49
|
45 arma_aligned uword aux_uword_a; //!< storage of auxiliary data, uword format
|
Chris@49
|
46 arma_aligned uword aux_uword_b; //!< storage of auxiliary data, uword format
|
Chris@49
|
47 arma_aligned uword aux_uword_c; //!< storage of auxiliary data, uword format
|
Chris@49
|
48
|
Chris@49
|
49 };
|
Chris@49
|
50
|
Chris@49
|
51
|
Chris@49
|
52
|
Chris@49
|
53 //! @}
|