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