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_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 Base<typename T1::elem_type,T1>&
|
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 X;
|
max@0
|
29 }
|
max@0
|
30
|
max@0
|
31
|
max@0
|
32
|
max@0
|
33 //! Base + scalar
|
max@0
|
34 template<typename T1>
|
max@0
|
35 arma_inline
|
max@0
|
36 const eOp<T1, eop_scalar_plus>
|
max@0
|
37 operator+
|
max@0
|
38 (const Base<typename T1::elem_type,T1>& X, const typename T1::elem_type k)
|
max@0
|
39 {
|
max@0
|
40 arma_extra_debug_sigprint();
|
max@0
|
41
|
max@0
|
42 return eOp<T1, eop_scalar_plus>(X.get_ref(), k);
|
max@0
|
43 }
|
max@0
|
44
|
max@0
|
45
|
max@0
|
46
|
max@0
|
47 //! scalar + Base
|
max@0
|
48 template<typename T1>
|
max@0
|
49 arma_inline
|
max@0
|
50 const eOp<T1, eop_scalar_plus>
|
max@0
|
51 operator+
|
max@0
|
52 (const typename T1::elem_type k, const Base<typename T1::elem_type,T1>& X)
|
max@0
|
53 {
|
max@0
|
54 arma_extra_debug_sigprint();
|
max@0
|
55
|
max@0
|
56 return eOp<T1, eop_scalar_plus>(X.get_ref(), k); // NOTE: order is swapped
|
max@0
|
57 }
|
max@0
|
58
|
max@0
|
59
|
max@0
|
60
|
max@0
|
61 //! non-complex Base + complex scalar (experimental)
|
max@0
|
62 template<typename T1>
|
max@0
|
63 arma_inline
|
max@0
|
64 const mtOp<typename std::complex<typename T1::pod_type>, T1, op_cx_scalar_plus>
|
max@0
|
65 operator+
|
max@0
|
66 (
|
max@0
|
67 const Base<typename T1::pod_type, T1>& X,
|
max@0
|
68 const std::complex<typename T1::pod_type>& k
|
max@0
|
69 )
|
max@0
|
70 {
|
max@0
|
71 arma_extra_debug_sigprint();
|
max@0
|
72
|
max@0
|
73 return mtOp<typename std::complex<typename T1::pod_type>, T1, op_cx_scalar_plus>('j', X.get_ref(), k);
|
max@0
|
74 }
|
max@0
|
75
|
max@0
|
76
|
max@0
|
77
|
max@0
|
78 //! complex scalar + non-complex Base (experimental)
|
max@0
|
79 template<typename T1>
|
max@0
|
80 arma_inline
|
max@0
|
81 const mtOp<typename std::complex<typename T1::pod_type>, T1, op_cx_scalar_plus>
|
max@0
|
82 operator+
|
max@0
|
83 (
|
max@0
|
84 const std::complex<typename T1::pod_type>& k,
|
max@0
|
85 const Base<typename T1::pod_type, T1>& X
|
max@0
|
86 )
|
max@0
|
87 {
|
max@0
|
88 arma_extra_debug_sigprint();
|
max@0
|
89
|
max@0
|
90 return mtOp<typename std::complex<typename T1::pod_type>, T1, op_cx_scalar_plus>('j', X.get_ref(), k); // NOTE: order is swapped
|
max@0
|
91 }
|
max@0
|
92
|
max@0
|
93
|
max@0
|
94
|
max@0
|
95 //! addition of Base objects with same element type
|
max@0
|
96 template<typename T1, typename T2>
|
max@0
|
97 arma_inline
|
max@0
|
98 const eGlue<T1, T2, eglue_plus>
|
max@0
|
99 operator+
|
max@0
|
100 (
|
max@0
|
101 const Base<typename T1::elem_type,T1>& X,
|
max@0
|
102 const Base<typename T1::elem_type,T2>& Y
|
max@0
|
103 )
|
max@0
|
104 {
|
max@0
|
105 arma_extra_debug_sigprint();
|
max@0
|
106
|
max@0
|
107 return eGlue<T1, T2, eglue_plus>(X.get_ref(), Y.get_ref());
|
max@0
|
108 }
|
max@0
|
109
|
max@0
|
110
|
max@0
|
111
|
max@0
|
112 //! addition of Base objects with different element types
|
max@0
|
113 template<typename T1, typename T2>
|
max@0
|
114 inline
|
max@0
|
115 const mtGlue<typename promote_type<typename T1::elem_type, typename T2::elem_type>::result, T1, T2, glue_mixed_plus>
|
max@0
|
116 operator+
|
max@0
|
117 (
|
max@0
|
118 const Base< typename force_different_type<typename T1::elem_type, typename T2::elem_type>::T1_result, T1>& X,
|
max@0
|
119 const Base< typename force_different_type<typename T1::elem_type, typename T2::elem_type>::T2_result, T2>& Y
|
max@0
|
120 )
|
max@0
|
121 {
|
max@0
|
122 arma_extra_debug_sigprint();
|
max@0
|
123
|
max@0
|
124 typedef typename T1::elem_type eT1;
|
max@0
|
125 typedef typename T2::elem_type eT2;
|
max@0
|
126
|
max@0
|
127 typedef typename promote_type<eT1,eT2>::result out_eT;
|
max@0
|
128
|
max@0
|
129 promote_type<eT1,eT2>::check();
|
max@0
|
130
|
max@0
|
131 return mtGlue<out_eT, T1, T2, glue_mixed_plus>( X.get_ref(), Y.get_ref() );
|
max@0
|
132 }
|
max@0
|
133
|
max@0
|
134
|
max@0
|
135
|
max@0
|
136 //! @}
|