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_minus
|
max@0
|
15 //! @{
|
max@0
|
16
|
max@0
|
17
|
max@0
|
18
|
max@0
|
19 //! unary -
|
max@0
|
20 template<typename T1>
|
max@0
|
21 arma_inline
|
max@0
|
22 const eOpCube<T1, eop_neg>
|
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 eOpCube<T1, eop_neg>(X.get_ref());
|
max@0
|
31 }
|
max@0
|
32
|
max@0
|
33
|
max@0
|
34
|
max@0
|
35 //! cancellation of two consecutive negations: -(-T1)
|
max@0
|
36 template<typename T1>
|
max@0
|
37 arma_inline
|
max@0
|
38 const T1&
|
max@0
|
39 operator-
|
max@0
|
40 (
|
max@0
|
41 const eOpCube<T1, eop_neg>& X
|
max@0
|
42 )
|
max@0
|
43 {
|
max@0
|
44 arma_extra_debug_sigprint();
|
max@0
|
45
|
max@0
|
46 return X.m;
|
max@0
|
47 }
|
max@0
|
48
|
max@0
|
49
|
max@0
|
50
|
max@0
|
51 //! BaseCube - scalar
|
max@0
|
52 template<typename T1>
|
max@0
|
53 arma_inline
|
max@0
|
54 const eOpCube<T1, eop_scalar_minus_post>
|
max@0
|
55 operator-
|
max@0
|
56 (
|
max@0
|
57 const BaseCube<typename T1::elem_type,T1>& X,
|
max@0
|
58 const typename T1::elem_type k
|
max@0
|
59 )
|
max@0
|
60 {
|
max@0
|
61 arma_extra_debug_sigprint();
|
max@0
|
62
|
max@0
|
63 return eOpCube<T1, eop_scalar_minus_post>(X.get_ref(), k);
|
max@0
|
64 }
|
max@0
|
65
|
max@0
|
66
|
max@0
|
67
|
max@0
|
68 //! scalar - BaseCube
|
max@0
|
69 template<typename T1>
|
max@0
|
70 arma_inline
|
max@0
|
71 const eOpCube<T1, eop_scalar_minus_pre>
|
max@0
|
72 operator-
|
max@0
|
73 (
|
max@0
|
74 const typename T1::elem_type k,
|
max@0
|
75 const BaseCube<typename T1::elem_type,T1>& X
|
max@0
|
76 )
|
max@0
|
77 {
|
max@0
|
78 arma_extra_debug_sigprint();
|
max@0
|
79
|
max@0
|
80 return eOpCube<T1, eop_scalar_minus_pre>(X.get_ref(), k);
|
max@0
|
81 }
|
max@0
|
82
|
max@0
|
83
|
max@0
|
84
|
max@0
|
85 //! complex scalar - non-complex BaseCube (experimental)
|
max@0
|
86 template<typename T1>
|
max@0
|
87 arma_inline
|
max@0
|
88 const mtOpCube<typename std::complex<typename T1::pod_type>, T1, op_cx_scalar_minus_pre>
|
max@0
|
89 operator-
|
max@0
|
90 (
|
max@0
|
91 const std::complex<typename T1::pod_type>& k,
|
max@0
|
92 const BaseCube<typename T1::pod_type, T1>& X
|
max@0
|
93 )
|
max@0
|
94 {
|
max@0
|
95 arma_extra_debug_sigprint();
|
max@0
|
96
|
max@0
|
97 return mtOpCube<typename std::complex<typename T1::pod_type>, T1, op_cx_scalar_minus_pre>('j', X.get_ref(), k);
|
max@0
|
98 }
|
max@0
|
99
|
max@0
|
100
|
max@0
|
101
|
max@0
|
102 //! non-complex BaseCube - complex scalar (experimental)
|
max@0
|
103 template<typename T1>
|
max@0
|
104 arma_inline
|
max@0
|
105 const mtOpCube<typename std::complex<typename T1::pod_type>, T1, op_cx_scalar_minus_post>
|
max@0
|
106 operator-
|
max@0
|
107 (
|
max@0
|
108 const BaseCube<typename T1::pod_type, T1>& X,
|
max@0
|
109 const std::complex<typename T1::pod_type>& k
|
max@0
|
110 )
|
max@0
|
111 {
|
max@0
|
112 arma_extra_debug_sigprint();
|
max@0
|
113
|
max@0
|
114 return mtOpCube<typename std::complex<typename T1::pod_type>, T1, op_cx_scalar_minus_post>('j', X.get_ref(), k);
|
max@0
|
115 }
|
max@0
|
116
|
max@0
|
117
|
max@0
|
118
|
max@0
|
119 //! subtraction of BaseCube objects with same element type
|
max@0
|
120 template<typename T1, typename T2>
|
max@0
|
121 arma_inline
|
max@0
|
122 const eGlueCube<T1, T2, eglue_minus>
|
max@0
|
123 operator-
|
max@0
|
124 (
|
max@0
|
125 const BaseCube<typename T1::elem_type,T1>& X,
|
max@0
|
126 const BaseCube<typename T1::elem_type,T2>& Y
|
max@0
|
127 )
|
max@0
|
128 {
|
max@0
|
129 arma_extra_debug_sigprint();
|
max@0
|
130
|
max@0
|
131 return eGlueCube<T1, T2, eglue_minus>(X.get_ref(), Y.get_ref());
|
max@0
|
132 }
|
max@0
|
133
|
max@0
|
134
|
max@0
|
135
|
max@0
|
136 //! subtraction of BaseCube objects with different element types
|
max@0
|
137 template<typename T1, typename T2>
|
max@0
|
138 inline
|
max@0
|
139 const mtGlueCube<typename promote_type<typename T1::elem_type, typename T2::elem_type>::result, T1, T2, glue_mixed_minus>
|
max@0
|
140 operator-
|
max@0
|
141 (
|
max@0
|
142 const BaseCube< typename force_different_type<typename T1::elem_type, typename T2::elem_type>::T1_result, T1>& X,
|
max@0
|
143 const BaseCube< typename force_different_type<typename T1::elem_type, typename T2::elem_type>::T2_result, T2>& Y
|
max@0
|
144 )
|
max@0
|
145 {
|
max@0
|
146 arma_extra_debug_sigprint();
|
max@0
|
147
|
max@0
|
148 typedef typename T1::elem_type eT1;
|
max@0
|
149 typedef typename T2::elem_type eT2;
|
max@0
|
150
|
max@0
|
151 typedef typename promote_type<eT1,eT2>::result out_eT;
|
max@0
|
152
|
max@0
|
153 promote_type<eT1,eT2>::check();
|
max@0
|
154
|
max@0
|
155 return mtGlueCube<out_eT, T1, T2, glue_mixed_minus>( X.get_ref(), Y.get_ref() );
|
max@0
|
156 }
|
max@0
|
157
|
max@0
|
158
|
max@0
|
159
|
max@0
|
160 //! @}
|