max@0
|
1 // Copyright (C) 2008-2010 NICTA (www.nicta.com.au)
|
max@0
|
2 // Copyright (C) 2008-2010 Conrad Sanderson
|
max@0
|
3 //
|
max@0
|
4 // This file is part of the Armadillo C++ library.
|
max@0
|
5 // It is provided without any warranty of fitness
|
max@0
|
6 // for any purpose. You can redistribute this file
|
max@0
|
7 // and/or modify it under the terms of the GNU
|
max@0
|
8 // Lesser General Public License (LGPL) as published
|
max@0
|
9 // by the Free Software Foundation, either version 3
|
max@0
|
10 // of the License or (at your option) any later version.
|
max@0
|
11 // (see http://www.opensource.org/licenses for more info)
|
max@0
|
12
|
max@0
|
13
|
max@0
|
14 //! \addtogroup operator_cube_plus
|
max@0
|
15 //! @{
|
max@0
|
16
|
max@0
|
17
|
max@0
|
18
|
max@0
|
19 //! unary plus operation (does nothing, but is required for completeness)
|
max@0
|
20 template<typename T1>
|
max@0
|
21 arma_inline
|
max@0
|
22 const BaseCube<typename T1::elem_type,T1>&
|
max@0
|
23 operator+
|
max@0
|
24 (
|
max@0
|
25 const BaseCube<typename T1::elem_type,T1>& X
|
max@0
|
26 )
|
max@0
|
27 {
|
max@0
|
28 arma_extra_debug_sigprint();
|
max@0
|
29
|
max@0
|
30 return X;
|
max@0
|
31 }
|
max@0
|
32
|
max@0
|
33
|
max@0
|
34
|
max@0
|
35 //! BaseCube + scalar
|
max@0
|
36 template<typename T1>
|
max@0
|
37 arma_inline
|
max@0
|
38 const eOpCube<T1, eop_scalar_plus>
|
max@0
|
39 operator+
|
max@0
|
40 (
|
max@0
|
41 const BaseCube<typename T1::elem_type,T1>& X,
|
max@0
|
42 const typename T1::elem_type k
|
max@0
|
43 )
|
max@0
|
44 {
|
max@0
|
45 arma_extra_debug_sigprint();
|
max@0
|
46
|
max@0
|
47 return eOpCube<T1, eop_scalar_plus>(X.get_ref(), k);
|
max@0
|
48 }
|
max@0
|
49
|
max@0
|
50
|
max@0
|
51
|
max@0
|
52 //! scalar + BaseCube
|
max@0
|
53 template<typename T1>
|
max@0
|
54 arma_inline
|
max@0
|
55 const eOpCube<T1, eop_scalar_plus>
|
max@0
|
56 operator+
|
max@0
|
57 (
|
max@0
|
58 const typename T1::elem_type k,
|
max@0
|
59 const BaseCube<typename T1::elem_type,T1>& X
|
max@0
|
60 )
|
max@0
|
61 {
|
max@0
|
62 arma_extra_debug_sigprint();
|
max@0
|
63
|
max@0
|
64 return eOpCube<T1, eop_scalar_plus>(X.get_ref(), k);
|
max@0
|
65 }
|
max@0
|
66
|
max@0
|
67
|
max@0
|
68
|
max@0
|
69 //! non-complex BaseCube + complex scalar (experimental)
|
max@0
|
70 template<typename T1>
|
max@0
|
71 arma_inline
|
max@0
|
72 const mtOpCube<typename std::complex<typename T1::pod_type>, T1, op_cx_scalar_plus>
|
max@0
|
73 operator+
|
max@0
|
74 (
|
max@0
|
75 const BaseCube<typename T1::pod_type, T1>& X,
|
max@0
|
76 const std::complex<typename T1::pod_type>& k
|
max@0
|
77 )
|
max@0
|
78 {
|
max@0
|
79 arma_extra_debug_sigprint();
|
max@0
|
80
|
max@0
|
81 return mtOpCube<typename std::complex<typename T1::pod_type>, T1, op_cx_scalar_plus>('j', X.get_ref(), k);
|
max@0
|
82 }
|
max@0
|
83
|
max@0
|
84
|
max@0
|
85
|
max@0
|
86 //! complex scalar + non-complex BaseCube (experimental)
|
max@0
|
87 template<typename T1>
|
max@0
|
88 arma_inline
|
max@0
|
89 const mtOpCube<typename std::complex<typename T1::pod_type>, T1, op_cx_scalar_plus>
|
max@0
|
90 operator+
|
max@0
|
91 (
|
max@0
|
92 const std::complex<typename T1::pod_type>& k,
|
max@0
|
93 const BaseCube<typename T1::pod_type, T1>& X
|
max@0
|
94 )
|
max@0
|
95 {
|
max@0
|
96 arma_extra_debug_sigprint();
|
max@0
|
97
|
max@0
|
98 return mtOpCube<typename std::complex<typename T1::pod_type>, T1, op_cx_scalar_plus>('j', X.get_ref(), k); // NOTE: order is swapped
|
max@0
|
99 }
|
max@0
|
100
|
max@0
|
101
|
max@0
|
102
|
max@0
|
103 //! addition of BaseCube objects with same element type
|
max@0
|
104 template<typename T1, typename T2>
|
max@0
|
105 arma_inline
|
max@0
|
106 const eGlueCube<T1, T2, eglue_plus>
|
max@0
|
107 operator+
|
max@0
|
108 (
|
max@0
|
109 const BaseCube<typename T1::elem_type,T1>& X,
|
max@0
|
110 const BaseCube<typename T1::elem_type,T2>& Y
|
max@0
|
111 )
|
max@0
|
112 {
|
max@0
|
113 arma_extra_debug_sigprint();
|
max@0
|
114
|
max@0
|
115 return eGlueCube<T1, T2, eglue_plus>(X.get_ref(), Y.get_ref());
|
max@0
|
116 }
|
max@0
|
117
|
max@0
|
118
|
max@0
|
119
|
max@0
|
120 //! addition of BaseCube objects with different element types
|
max@0
|
121 template<typename T1, typename T2>
|
max@0
|
122 inline
|
max@0
|
123 const mtGlueCube<typename promote_type<typename T1::elem_type, typename T2::elem_type>::result, T1, T2, glue_mixed_plus>
|
max@0
|
124 operator+
|
max@0
|
125 (
|
max@0
|
126 const BaseCube< typename force_different_type<typename T1::elem_type, typename T2::elem_type>::T1_result, T1>& X,
|
max@0
|
127 const BaseCube< typename force_different_type<typename T1::elem_type, typename T2::elem_type>::T2_result, T2>& Y
|
max@0
|
128 )
|
max@0
|
129 {
|
max@0
|
130 arma_extra_debug_sigprint();
|
max@0
|
131
|
max@0
|
132 typedef typename T1::elem_type eT1;
|
max@0
|
133 typedef typename T2::elem_type eT2;
|
max@0
|
134
|
max@0
|
135 typedef typename promote_type<eT1,eT2>::result out_eT;
|
max@0
|
136
|
max@0
|
137 promote_type<eT1,eT2>::check();
|
max@0
|
138
|
max@0
|
139 return mtGlueCube<out_eT, T1, T2, glue_mixed_plus>( X.get_ref(), Y.get_ref() );
|
max@0
|
140 }
|
max@0
|
141
|
max@0
|
142
|
max@0
|
143
|
max@0
|
144 //! @}
|