Eurobot2012_Secondary

Fork of Eurobot_2012_Secondary by Shuto Naruse

Committer:
narshu
Date:
Wed Oct 17 22:25:31 2012 +0000
Revision:
1:cc2a9eb0bd55
Commit before publishing

Who changed what in which revision?

UserRevisionLine numberNew contents of line
narshu 1:cc2a9eb0bd55 1 /*
narshu 1:cc2a9eb0bd55 2 * Tiny Vector Matrix Library
narshu 1:cc2a9eb0bd55 3 * Dense Vector Matrix Libary of Tiny size using Expression Templates
narshu 1:cc2a9eb0bd55 4 *
narshu 1:cc2a9eb0bd55 5 * Copyright (C) 2001 - 2007 Olaf Petzold <opetzold@users.sourceforge.net>
narshu 1:cc2a9eb0bd55 6 *
narshu 1:cc2a9eb0bd55 7 * This library is free software; you can redistribute it and/or
narshu 1:cc2a9eb0bd55 8 * modify it under the terms of the GNU lesser General Public
narshu 1:cc2a9eb0bd55 9 * License as published by the Free Software Foundation; either
narshu 1:cc2a9eb0bd55 10 * version 2.1 of the License, or (at your option) any later version.
narshu 1:cc2a9eb0bd55 11 *
narshu 1:cc2a9eb0bd55 12 * This library is distributed in the hope that it will be useful,
narshu 1:cc2a9eb0bd55 13 * but WITHOUT ANY WARRANTY; without even the implied warranty of
narshu 1:cc2a9eb0bd55 14 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
narshu 1:cc2a9eb0bd55 15 * lesser General Public License for more details.
narshu 1:cc2a9eb0bd55 16 *
narshu 1:cc2a9eb0bd55 17 * You should have received a copy of the GNU lesser General Public
narshu 1:cc2a9eb0bd55 18 * License along with this library; if not, write to the Free Software
narshu 1:cc2a9eb0bd55 19 * Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
narshu 1:cc2a9eb0bd55 20 *
narshu 1:cc2a9eb0bd55 21 * $Id: MatrixBinaryFunctions.h,v 1.16 2007-06-23 15:58:58 opetzold Exp $
narshu 1:cc2a9eb0bd55 22 */
narshu 1:cc2a9eb0bd55 23
narshu 1:cc2a9eb0bd55 24 #ifndef TVMET_MATRIX_BINARY_FUNCTIONS_H
narshu 1:cc2a9eb0bd55 25 #define TVMET_MATRIX_BINARY_FUNCTIONS_H
narshu 1:cc2a9eb0bd55 26
narshu 1:cc2a9eb0bd55 27 namespace tvmet {
narshu 1:cc2a9eb0bd55 28
narshu 1:cc2a9eb0bd55 29 /*********************************************************
narshu 1:cc2a9eb0bd55 30 * PART I: DECLARATION
narshu 1:cc2a9eb0bd55 31 *********************************************************/
narshu 1:cc2a9eb0bd55 32
narshu 1:cc2a9eb0bd55 33 /*
narshu 1:cc2a9eb0bd55 34 * binary_function(Matrix<T1, Rows, Cols>, Matrix<T2, Rows, Cols>)
narshu 1:cc2a9eb0bd55 35 * binary_function(Matrix<T1, Rows, Cols>, XprMatrix<E, Rows, Cols>)
narshu 1:cc2a9eb0bd55 36 * binary_function(XprMatrix<E, Rows, Cols>, Matrix<T, Rows, Cols>)
narshu 1:cc2a9eb0bd55 37 */
narshu 1:cc2a9eb0bd55 38 #define TVMET_DECLARE_MACRO(NAME) \
narshu 1:cc2a9eb0bd55 39 template<class T1, class T2, std::size_t Rows, std::size_t Cols> \
narshu 1:cc2a9eb0bd55 40 XprMatrix< \
narshu 1:cc2a9eb0bd55 41 XprBinOp< \
narshu 1:cc2a9eb0bd55 42 Fcnl_##NAME<T1, T2>, \
narshu 1:cc2a9eb0bd55 43 MatrixConstReference<T1, Rows, Cols>, \
narshu 1:cc2a9eb0bd55 44 MatrixConstReference<T2, Rows, Cols> \
narshu 1:cc2a9eb0bd55 45 >, \
narshu 1:cc2a9eb0bd55 46 Rows, Cols \
narshu 1:cc2a9eb0bd55 47 > \
narshu 1:cc2a9eb0bd55 48 NAME(const Matrix<T1, Rows, Cols>& lhs, \
narshu 1:cc2a9eb0bd55 49 const Matrix<T2, Cols, Cols>& rhs) TVMET_CXX_ALWAYS_INLINE; \
narshu 1:cc2a9eb0bd55 50 \
narshu 1:cc2a9eb0bd55 51 template<class E, class T, std::size_t Rows, std::size_t Cols> \
narshu 1:cc2a9eb0bd55 52 XprMatrix< \
narshu 1:cc2a9eb0bd55 53 XprBinOp< \
narshu 1:cc2a9eb0bd55 54 Fcnl_##NAME<typename E::value_type, T>, \
narshu 1:cc2a9eb0bd55 55 MatrixConstReference<T, Rows, Cols>, \
narshu 1:cc2a9eb0bd55 56 XprMatrix<E, Rows, Cols> \
narshu 1:cc2a9eb0bd55 57 >, \
narshu 1:cc2a9eb0bd55 58 Rows, Cols \
narshu 1:cc2a9eb0bd55 59 > \
narshu 1:cc2a9eb0bd55 60 NAME(const XprMatrix<E, Rows, Cols>& lhs, \
narshu 1:cc2a9eb0bd55 61 const Matrix<T, Rows, Cols>& rhs) TVMET_CXX_ALWAYS_INLINE; \
narshu 1:cc2a9eb0bd55 62 \
narshu 1:cc2a9eb0bd55 63 template<class E, class T, std::size_t Rows, std::size_t Cols> \
narshu 1:cc2a9eb0bd55 64 XprMatrix< \
narshu 1:cc2a9eb0bd55 65 XprBinOp< \
narshu 1:cc2a9eb0bd55 66 Fcnl_##NAME<T, typename E::value_type>, \
narshu 1:cc2a9eb0bd55 67 MatrixConstReference<T, Rows, Cols>, \
narshu 1:cc2a9eb0bd55 68 XprMatrix<E, Rows, Cols> \
narshu 1:cc2a9eb0bd55 69 >, \
narshu 1:cc2a9eb0bd55 70 Rows, Cols \
narshu 1:cc2a9eb0bd55 71 > \
narshu 1:cc2a9eb0bd55 72 NAME(const Matrix<T, Rows, Cols>& lhs, \
narshu 1:cc2a9eb0bd55 73 const XprMatrix<E, Rows, Cols>& rhs) TVMET_CXX_ALWAYS_INLINE;
narshu 1:cc2a9eb0bd55 74
narshu 1:cc2a9eb0bd55 75 TVMET_DECLARE_MACRO(atan2)
narshu 1:cc2a9eb0bd55 76 TVMET_DECLARE_MACRO(drem)
narshu 1:cc2a9eb0bd55 77 TVMET_DECLARE_MACRO(fmod)
narshu 1:cc2a9eb0bd55 78 TVMET_DECLARE_MACRO(hypot)
narshu 1:cc2a9eb0bd55 79 TVMET_DECLARE_MACRO(jn)
narshu 1:cc2a9eb0bd55 80 TVMET_DECLARE_MACRO(yn)
narshu 1:cc2a9eb0bd55 81 TVMET_DECLARE_MACRO(pow)
narshu 1:cc2a9eb0bd55 82 #if defined(TVMET_HAVE_COMPLEX)
narshu 1:cc2a9eb0bd55 83 TVMET_DECLARE_MACRO(polar)
narshu 1:cc2a9eb0bd55 84 #endif
narshu 1:cc2a9eb0bd55 85
narshu 1:cc2a9eb0bd55 86 #undef TVMET_DECLARE_MACRO
narshu 1:cc2a9eb0bd55 87
narshu 1:cc2a9eb0bd55 88
narshu 1:cc2a9eb0bd55 89 /*
narshu 1:cc2a9eb0bd55 90 * binary_function(Matrix<T, Rows, Cols>, POD)
narshu 1:cc2a9eb0bd55 91 */
narshu 1:cc2a9eb0bd55 92 #define TVMET_DECLARE_MACRO(NAME, TP) \
narshu 1:cc2a9eb0bd55 93 template<class T, std::size_t Rows, std::size_t Cols> \
narshu 1:cc2a9eb0bd55 94 XprMatrix< \
narshu 1:cc2a9eb0bd55 95 XprBinOp< \
narshu 1:cc2a9eb0bd55 96 Fcnl_##NAME<T, TP >, \
narshu 1:cc2a9eb0bd55 97 MatrixConstReference<T, Rows, Cols>, \
narshu 1:cc2a9eb0bd55 98 XprLiteral< TP > \
narshu 1:cc2a9eb0bd55 99 >, \
narshu 1:cc2a9eb0bd55 100 Rows, Cols \
narshu 1:cc2a9eb0bd55 101 > \
narshu 1:cc2a9eb0bd55 102 NAME(const Matrix<T, Rows, Cols>& lhs, TP rhs) TVMET_CXX_ALWAYS_INLINE;
narshu 1:cc2a9eb0bd55 103
narshu 1:cc2a9eb0bd55 104 TVMET_DECLARE_MACRO(atan2, int)
narshu 1:cc2a9eb0bd55 105 TVMET_DECLARE_MACRO(drem, int)
narshu 1:cc2a9eb0bd55 106 TVMET_DECLARE_MACRO(fmod, int)
narshu 1:cc2a9eb0bd55 107 TVMET_DECLARE_MACRO(hypot, int)
narshu 1:cc2a9eb0bd55 108 TVMET_DECLARE_MACRO(jn, int)
narshu 1:cc2a9eb0bd55 109 TVMET_DECLARE_MACRO(yn, int)
narshu 1:cc2a9eb0bd55 110 TVMET_DECLARE_MACRO(pow, int)
narshu 1:cc2a9eb0bd55 111
narshu 1:cc2a9eb0bd55 112 #if defined(TVMET_HAVE_LONG_LONG)
narshu 1:cc2a9eb0bd55 113 TVMET_DECLARE_MACRO(atan2, long long int)
narshu 1:cc2a9eb0bd55 114 TVMET_DECLARE_MACRO(drem, long long int)
narshu 1:cc2a9eb0bd55 115 TVMET_DECLARE_MACRO(fmod, long long int)
narshu 1:cc2a9eb0bd55 116 TVMET_DECLARE_MACRO(hypot, long long int)
narshu 1:cc2a9eb0bd55 117 TVMET_DECLARE_MACRO(jn, long long int)
narshu 1:cc2a9eb0bd55 118 TVMET_DECLARE_MACRO(yn, long long int)
narshu 1:cc2a9eb0bd55 119 TVMET_DECLARE_MACRO(pow, long long int)
narshu 1:cc2a9eb0bd55 120 #endif // defined(TVMET_HAVE_LONG_LONG)
narshu 1:cc2a9eb0bd55 121
narshu 1:cc2a9eb0bd55 122 TVMET_DECLARE_MACRO(atan2, float)
narshu 1:cc2a9eb0bd55 123 TVMET_DECLARE_MACRO(drem, float)
narshu 1:cc2a9eb0bd55 124 TVMET_DECLARE_MACRO(fmod, float)
narshu 1:cc2a9eb0bd55 125 TVMET_DECLARE_MACRO(hypot, float)
narshu 1:cc2a9eb0bd55 126 TVMET_DECLARE_MACRO(jn, float)
narshu 1:cc2a9eb0bd55 127 TVMET_DECLARE_MACRO(yn, float)
narshu 1:cc2a9eb0bd55 128 TVMET_DECLARE_MACRO(pow, float)
narshu 1:cc2a9eb0bd55 129
narshu 1:cc2a9eb0bd55 130 TVMET_DECLARE_MACRO(atan2, double)
narshu 1:cc2a9eb0bd55 131 TVMET_DECLARE_MACRO(drem, double)
narshu 1:cc2a9eb0bd55 132 TVMET_DECLARE_MACRO(fmod, double)
narshu 1:cc2a9eb0bd55 133 TVMET_DECLARE_MACRO(hypot, double)
narshu 1:cc2a9eb0bd55 134 TVMET_DECLARE_MACRO(jn, double)
narshu 1:cc2a9eb0bd55 135 TVMET_DECLARE_MACRO(yn, double)
narshu 1:cc2a9eb0bd55 136 TVMET_DECLARE_MACRO(pow, double)
narshu 1:cc2a9eb0bd55 137
narshu 1:cc2a9eb0bd55 138 #if defined(TVMET_HAVE_LONG_DOUBLE)
narshu 1:cc2a9eb0bd55 139 TVMET_DECLARE_MACRO(atan2, long double)
narshu 1:cc2a9eb0bd55 140 TVMET_DECLARE_MACRO(drem, long double)
narshu 1:cc2a9eb0bd55 141 TVMET_DECLARE_MACRO(fmod, long double)
narshu 1:cc2a9eb0bd55 142 TVMET_DECLARE_MACRO(hypot, long double)
narshu 1:cc2a9eb0bd55 143 TVMET_DECLARE_MACRO(jn, long double)
narshu 1:cc2a9eb0bd55 144 TVMET_DECLARE_MACRO(yn, long double)
narshu 1:cc2a9eb0bd55 145 TVMET_DECLARE_MACRO(pow, long double)
narshu 1:cc2a9eb0bd55 146 #endif // defined(TVMET_HAVE_LONG_DOUBLE)
narshu 1:cc2a9eb0bd55 147
narshu 1:cc2a9eb0bd55 148 #undef TVMET_DECLARE_MACRO
narshu 1:cc2a9eb0bd55 149
narshu 1:cc2a9eb0bd55 150
narshu 1:cc2a9eb0bd55 151 /*
narshu 1:cc2a9eb0bd55 152 * complex math
narshu 1:cc2a9eb0bd55 153 */
narshu 1:cc2a9eb0bd55 154
narshu 1:cc2a9eb0bd55 155 #if defined(TVMET_HAVE_COMPLEX) && defined(TVMET_HAVE_COMPLEX_MATH1)
narshu 1:cc2a9eb0bd55 156 template<class T, std::size_t Rows, std::size_t Cols>
narshu 1:cc2a9eb0bd55 157 XprMatrix<
narshu 1:cc2a9eb0bd55 158 XprBinOp<
narshu 1:cc2a9eb0bd55 159 Fcnl_pow<T, std::complex<T> >,
narshu 1:cc2a9eb0bd55 160 MatrixConstReference<T, Rows, Cols>,
narshu 1:cc2a9eb0bd55 161 XprLiteral< std::complex<T> >
narshu 1:cc2a9eb0bd55 162 >,
narshu 1:cc2a9eb0bd55 163 Rows, Cols
narshu 1:cc2a9eb0bd55 164 >
narshu 1:cc2a9eb0bd55 165 pow(const Matrix<T, Rows, Cols>& lhs,
narshu 1:cc2a9eb0bd55 166 const std::complex<T>& rhs) TVMET_CXX_ALWAYS_INLINE;
narshu 1:cc2a9eb0bd55 167
narshu 1:cc2a9eb0bd55 168
narshu 1:cc2a9eb0bd55 169 template<class T, std::size_t Rows, std::size_t Cols>
narshu 1:cc2a9eb0bd55 170 XprMatrix<
narshu 1:cc2a9eb0bd55 171 XprBinOp<
narshu 1:cc2a9eb0bd55 172 Fcnl_pow< std::complex<T>, std::complex<T> >,
narshu 1:cc2a9eb0bd55 173 MatrixConstReference<std::complex<T>, Rows, Cols>,
narshu 1:cc2a9eb0bd55 174 XprLiteral< std::complex<T> >
narshu 1:cc2a9eb0bd55 175 >,
narshu 1:cc2a9eb0bd55 176 Rows, Cols
narshu 1:cc2a9eb0bd55 177 >
narshu 1:cc2a9eb0bd55 178 pow(const Matrix<std::complex<T>, Rows, Cols>& lhs,
narshu 1:cc2a9eb0bd55 179 const std::complex<T>& rhs) TVMET_CXX_ALWAYS_INLINE;
narshu 1:cc2a9eb0bd55 180
narshu 1:cc2a9eb0bd55 181
narshu 1:cc2a9eb0bd55 182 /**
narshu 1:cc2a9eb0bd55 183 * \fn pow(const Matrix<std::complex<T>, Rows, Cols>& lhs, const T& rhs)
narshu 1:cc2a9eb0bd55 184 * \ingroup _binary_function
narshu 1:cc2a9eb0bd55 185 */
narshu 1:cc2a9eb0bd55 186 template<class T, std::size_t Rows, std::size_t Cols>
narshu 1:cc2a9eb0bd55 187 XprMatrix<
narshu 1:cc2a9eb0bd55 188 XprBinOp<
narshu 1:cc2a9eb0bd55 189 Fcnl_pow<std::complex<T>, T>,
narshu 1:cc2a9eb0bd55 190 MatrixConstReference<std::complex<T>, Rows, Cols>,
narshu 1:cc2a9eb0bd55 191 XprLiteral<T>
narshu 1:cc2a9eb0bd55 192 >,
narshu 1:cc2a9eb0bd55 193 Rows, Cols
narshu 1:cc2a9eb0bd55 194 >
narshu 1:cc2a9eb0bd55 195 pow(const Matrix<std::complex<T>, Rows, Cols>& lhs,
narshu 1:cc2a9eb0bd55 196 const T& rhs) TVMET_CXX_ALWAYS_INLINE;
narshu 1:cc2a9eb0bd55 197
narshu 1:cc2a9eb0bd55 198
narshu 1:cc2a9eb0bd55 199 /**
narshu 1:cc2a9eb0bd55 200 * \fn pow(const Matrix<std::complex<T>, Rows, Cols>& lhs, int rhs)
narshu 1:cc2a9eb0bd55 201 * \ingroup _binary_function
narshu 1:cc2a9eb0bd55 202 */
narshu 1:cc2a9eb0bd55 203 template<class T, std::size_t Rows, std::size_t Cols>
narshu 1:cc2a9eb0bd55 204 XprMatrix<
narshu 1:cc2a9eb0bd55 205 XprBinOp<
narshu 1:cc2a9eb0bd55 206 Fcnl_pow<std::complex<T>, int>,
narshu 1:cc2a9eb0bd55 207 MatrixConstReference<std::complex<T>, Rows, Cols>,
narshu 1:cc2a9eb0bd55 208 XprLiteral<int>
narshu 1:cc2a9eb0bd55 209 >,
narshu 1:cc2a9eb0bd55 210 Rows, Cols
narshu 1:cc2a9eb0bd55 211 >
narshu 1:cc2a9eb0bd55 212 pow(const Matrix<std::complex<T>, Rows, Cols>& lhs,
narshu 1:cc2a9eb0bd55 213 int rhs) TVMET_CXX_ALWAYS_INLINE;
narshu 1:cc2a9eb0bd55 214
narshu 1:cc2a9eb0bd55 215
narshu 1:cc2a9eb0bd55 216 template<class T, std::size_t Rows, std::size_t Cols>
narshu 1:cc2a9eb0bd55 217 XprMatrix<
narshu 1:cc2a9eb0bd55 218 XprBinOp<
narshu 1:cc2a9eb0bd55 219 Fcnl_polar<T, T>,
narshu 1:cc2a9eb0bd55 220 MatrixConstReference<T, Rows, Cols>,
narshu 1:cc2a9eb0bd55 221 XprLiteral<T>
narshu 1:cc2a9eb0bd55 222 >,
narshu 1:cc2a9eb0bd55 223 Rows, Cols
narshu 1:cc2a9eb0bd55 224 >
narshu 1:cc2a9eb0bd55 225 polar(const Matrix<T, Rows, Cols>& lhs,
narshu 1:cc2a9eb0bd55 226 const T& rhs) TVMET_CXX_ALWAYS_INLINE;
narshu 1:cc2a9eb0bd55 227
narshu 1:cc2a9eb0bd55 228 #endif // defined(TVMET_HAVE_COMPLEX) && defined(TVMET_HAVE_COMPLEX_MATH1)
narshu 1:cc2a9eb0bd55 229
narshu 1:cc2a9eb0bd55 230
narshu 1:cc2a9eb0bd55 231 #if defined(TVMET_HAVE_COMPLEX) && defined(TVMET_HAVE_COMPLEX_MATH2)
narshu 1:cc2a9eb0bd55 232 // to be written (atan2)
narshu 1:cc2a9eb0bd55 233 #endif // defined(TVMET_HAVE_COMPLEX) && defined(TVMET_HAVE_COMPLEX_MATH2)
narshu 1:cc2a9eb0bd55 234
narshu 1:cc2a9eb0bd55 235
narshu 1:cc2a9eb0bd55 236 /*********************************************************
narshu 1:cc2a9eb0bd55 237 * PART II: IMPLEMENTATION
narshu 1:cc2a9eb0bd55 238 *********************************************************/
narshu 1:cc2a9eb0bd55 239
narshu 1:cc2a9eb0bd55 240 /*
narshu 1:cc2a9eb0bd55 241 * binary_function(Matrix<T1, Rows, Cols>, Matrix<T2, Rows, Cols>)
narshu 1:cc2a9eb0bd55 242 * binary_function(Matrix<T1, Rows, Cols>, XprMatrix<E, Rows, Cols>)
narshu 1:cc2a9eb0bd55 243 * binary_function(XprMatrix<E, Rows, Cols>, Matrix<T, Rows, Cols>)
narshu 1:cc2a9eb0bd55 244 */
narshu 1:cc2a9eb0bd55 245 #define TVMET_IMPLEMENT_MACRO(NAME) \
narshu 1:cc2a9eb0bd55 246 template<class T1, class T2, std::size_t Rows, std::size_t Cols> \
narshu 1:cc2a9eb0bd55 247 inline \
narshu 1:cc2a9eb0bd55 248 XprMatrix< \
narshu 1:cc2a9eb0bd55 249 XprBinOp< \
narshu 1:cc2a9eb0bd55 250 Fcnl_##NAME<T1, T2>, \
narshu 1:cc2a9eb0bd55 251 MatrixConstReference<T1, Rows, Cols>, \
narshu 1:cc2a9eb0bd55 252 MatrixConstReference<T2, Rows, Cols> \
narshu 1:cc2a9eb0bd55 253 >, \
narshu 1:cc2a9eb0bd55 254 Rows, Cols \
narshu 1:cc2a9eb0bd55 255 > \
narshu 1:cc2a9eb0bd55 256 NAME(const Matrix<T1, Rows, Cols>& lhs, const Matrix<T2, Cols, Cols>& rhs) { \
narshu 1:cc2a9eb0bd55 257 typedef XprBinOp < \
narshu 1:cc2a9eb0bd55 258 Fcnl_##NAME<T1, T2>, \
narshu 1:cc2a9eb0bd55 259 MatrixConstReference<T1, Rows, Cols>, \
narshu 1:cc2a9eb0bd55 260 MatrixConstReference<T2, Rows, Cols> \
narshu 1:cc2a9eb0bd55 261 > expr_type; \
narshu 1:cc2a9eb0bd55 262 return XprMatrix<expr_type, Rows, Cols>( \
narshu 1:cc2a9eb0bd55 263 expr_type(lhs.const_ref(), rhs.const_ref())); \
narshu 1:cc2a9eb0bd55 264 } \
narshu 1:cc2a9eb0bd55 265 \
narshu 1:cc2a9eb0bd55 266 template<class E, class T, std::size_t Rows, std::size_t Cols> \
narshu 1:cc2a9eb0bd55 267 inline \
narshu 1:cc2a9eb0bd55 268 XprMatrix< \
narshu 1:cc2a9eb0bd55 269 XprBinOp< \
narshu 1:cc2a9eb0bd55 270 Fcnl_##NAME<typename E::value_type, T>, \
narshu 1:cc2a9eb0bd55 271 MatrixConstReference<T, Rows, Cols>, \
narshu 1:cc2a9eb0bd55 272 XprMatrix<E, Rows, Cols> \
narshu 1:cc2a9eb0bd55 273 >, \
narshu 1:cc2a9eb0bd55 274 Rows, Cols \
narshu 1:cc2a9eb0bd55 275 > \
narshu 1:cc2a9eb0bd55 276 NAME(const XprMatrix<E, Rows, Cols>& lhs, const Matrix<T, Rows, Cols>& rhs) { \
narshu 1:cc2a9eb0bd55 277 typedef XprBinOp< \
narshu 1:cc2a9eb0bd55 278 Fcnl_##NAME<typename E::value_type, T>, \
narshu 1:cc2a9eb0bd55 279 XprMatrix<E, Rows, Cols>, \
narshu 1:cc2a9eb0bd55 280 MatrixConstReference<T, Rows, Cols> \
narshu 1:cc2a9eb0bd55 281 > expr_type; \
narshu 1:cc2a9eb0bd55 282 return XprMatrix<expr_type, Rows, Cols>( \
narshu 1:cc2a9eb0bd55 283 expr_type(lhs, rhs.const_ref())); \
narshu 1:cc2a9eb0bd55 284 } \
narshu 1:cc2a9eb0bd55 285 \
narshu 1:cc2a9eb0bd55 286 template<class E, class T, std::size_t Rows, std::size_t Cols> \
narshu 1:cc2a9eb0bd55 287 inline \
narshu 1:cc2a9eb0bd55 288 XprMatrix< \
narshu 1:cc2a9eb0bd55 289 XprBinOp< \
narshu 1:cc2a9eb0bd55 290 Fcnl_##NAME<T, typename E::value_type>, \
narshu 1:cc2a9eb0bd55 291 MatrixConstReference<T, Rows, Cols>, \
narshu 1:cc2a9eb0bd55 292 XprMatrix<E, Rows, Cols> \
narshu 1:cc2a9eb0bd55 293 >, \
narshu 1:cc2a9eb0bd55 294 Rows, Cols \
narshu 1:cc2a9eb0bd55 295 > \
narshu 1:cc2a9eb0bd55 296 NAME(const Matrix<T, Rows, Cols>& lhs, const XprMatrix<E, Rows, Cols>& rhs) { \
narshu 1:cc2a9eb0bd55 297 typedef XprBinOp< \
narshu 1:cc2a9eb0bd55 298 Fcnl_##NAME<T, typename E::value_type>, \
narshu 1:cc2a9eb0bd55 299 MatrixConstReference<T, Rows, Cols>, \
narshu 1:cc2a9eb0bd55 300 XprMatrix<E, Rows, Cols> \
narshu 1:cc2a9eb0bd55 301 > expr_type; \
narshu 1:cc2a9eb0bd55 302 return XprMatrix<expr_type, Rows, Cols>( \
narshu 1:cc2a9eb0bd55 303 expr_type(lhs.const_ref(), rhs)); \
narshu 1:cc2a9eb0bd55 304 }
narshu 1:cc2a9eb0bd55 305
narshu 1:cc2a9eb0bd55 306 TVMET_IMPLEMENT_MACRO(atan2)
narshu 1:cc2a9eb0bd55 307 TVMET_IMPLEMENT_MACRO(drem)
narshu 1:cc2a9eb0bd55 308 TVMET_IMPLEMENT_MACRO(fmod)
narshu 1:cc2a9eb0bd55 309 TVMET_IMPLEMENT_MACRO(hypot)
narshu 1:cc2a9eb0bd55 310 TVMET_IMPLEMENT_MACRO(jn)
narshu 1:cc2a9eb0bd55 311 TVMET_IMPLEMENT_MACRO(yn)
narshu 1:cc2a9eb0bd55 312 TVMET_IMPLEMENT_MACRO(pow)
narshu 1:cc2a9eb0bd55 313 #if defined(TVMET_HAVE_COMPLEX)
narshu 1:cc2a9eb0bd55 314 TVMET_IMPLEMENT_MACRO(polar)
narshu 1:cc2a9eb0bd55 315 #endif
narshu 1:cc2a9eb0bd55 316 #undef TVMET_IMPLEMENT_MACRO
narshu 1:cc2a9eb0bd55 317
narshu 1:cc2a9eb0bd55 318
narshu 1:cc2a9eb0bd55 319 /*
narshu 1:cc2a9eb0bd55 320 * binary_function(Matrix<T, Rows, Cols>, POD)
narshu 1:cc2a9eb0bd55 321 */
narshu 1:cc2a9eb0bd55 322 #define TVMET_IMPLEMENT_MACRO(NAME, TP) \
narshu 1:cc2a9eb0bd55 323 template<class T, std::size_t Rows, std::size_t Cols> \
narshu 1:cc2a9eb0bd55 324 inline \
narshu 1:cc2a9eb0bd55 325 XprMatrix< \
narshu 1:cc2a9eb0bd55 326 XprBinOp< \
narshu 1:cc2a9eb0bd55 327 Fcnl_##NAME<T, TP >, \
narshu 1:cc2a9eb0bd55 328 MatrixConstReference<T, Rows, Cols>, \
narshu 1:cc2a9eb0bd55 329 XprLiteral< TP > \
narshu 1:cc2a9eb0bd55 330 >, \
narshu 1:cc2a9eb0bd55 331 Rows, Cols \
narshu 1:cc2a9eb0bd55 332 > \
narshu 1:cc2a9eb0bd55 333 NAME(const Matrix<T, Rows, Cols>& lhs, TP rhs) { \
narshu 1:cc2a9eb0bd55 334 typedef XprBinOp< \
narshu 1:cc2a9eb0bd55 335 Fcnl_##NAME<T, TP >, \
narshu 1:cc2a9eb0bd55 336 MatrixConstReference<T, Rows, Cols>, \
narshu 1:cc2a9eb0bd55 337 XprLiteral< TP > \
narshu 1:cc2a9eb0bd55 338 > expr_type; \
narshu 1:cc2a9eb0bd55 339 return XprMatrix<expr_type, Rows, Cols>( \
narshu 1:cc2a9eb0bd55 340 expr_type(lhs.const_ref(), XprLiteral< TP >(rhs))); \
narshu 1:cc2a9eb0bd55 341 }
narshu 1:cc2a9eb0bd55 342
narshu 1:cc2a9eb0bd55 343 TVMET_IMPLEMENT_MACRO(atan2, int)
narshu 1:cc2a9eb0bd55 344 TVMET_IMPLEMENT_MACRO(drem, int)
narshu 1:cc2a9eb0bd55 345 TVMET_IMPLEMENT_MACRO(fmod, int)
narshu 1:cc2a9eb0bd55 346 TVMET_IMPLEMENT_MACRO(hypot, int)
narshu 1:cc2a9eb0bd55 347 TVMET_IMPLEMENT_MACRO(jn, int)
narshu 1:cc2a9eb0bd55 348 TVMET_IMPLEMENT_MACRO(yn, int)
narshu 1:cc2a9eb0bd55 349 TVMET_IMPLEMENT_MACRO(pow, int)
narshu 1:cc2a9eb0bd55 350
narshu 1:cc2a9eb0bd55 351 #if defined(TVMET_HAVE_LONG_LONG)
narshu 1:cc2a9eb0bd55 352 TVMET_IMPLEMENT_MACRO(atan2, long long int)
narshu 1:cc2a9eb0bd55 353 TVMET_IMPLEMENT_MACRO(drem, long long int)
narshu 1:cc2a9eb0bd55 354 TVMET_IMPLEMENT_MACRO(fmod, long long int)
narshu 1:cc2a9eb0bd55 355 TVMET_IMPLEMENT_MACRO(hypot, long long int)
narshu 1:cc2a9eb0bd55 356 TVMET_IMPLEMENT_MACRO(jn, long long int)
narshu 1:cc2a9eb0bd55 357 TVMET_IMPLEMENT_MACRO(yn, long long int)
narshu 1:cc2a9eb0bd55 358 TVMET_IMPLEMENT_MACRO(pow, long long int)
narshu 1:cc2a9eb0bd55 359 #endif // defined(TVMET_HAVE_LONG_LONG)
narshu 1:cc2a9eb0bd55 360
narshu 1:cc2a9eb0bd55 361 TVMET_IMPLEMENT_MACRO(atan2, float)
narshu 1:cc2a9eb0bd55 362 TVMET_IMPLEMENT_MACRO(drem, float)
narshu 1:cc2a9eb0bd55 363 TVMET_IMPLEMENT_MACRO(fmod, float)
narshu 1:cc2a9eb0bd55 364 TVMET_IMPLEMENT_MACRO(hypot, float)
narshu 1:cc2a9eb0bd55 365 TVMET_IMPLEMENT_MACRO(jn, float)
narshu 1:cc2a9eb0bd55 366 TVMET_IMPLEMENT_MACRO(yn, float)
narshu 1:cc2a9eb0bd55 367 TVMET_IMPLEMENT_MACRO(pow, float)
narshu 1:cc2a9eb0bd55 368
narshu 1:cc2a9eb0bd55 369 TVMET_IMPLEMENT_MACRO(atan2, double)
narshu 1:cc2a9eb0bd55 370 TVMET_IMPLEMENT_MACRO(drem, double)
narshu 1:cc2a9eb0bd55 371 TVMET_IMPLEMENT_MACRO(fmod, double)
narshu 1:cc2a9eb0bd55 372 TVMET_IMPLEMENT_MACRO(hypot, double)
narshu 1:cc2a9eb0bd55 373 TVMET_IMPLEMENT_MACRO(jn, double)
narshu 1:cc2a9eb0bd55 374 TVMET_IMPLEMENT_MACRO(yn, double)
narshu 1:cc2a9eb0bd55 375 TVMET_IMPLEMENT_MACRO(pow, double)
narshu 1:cc2a9eb0bd55 376
narshu 1:cc2a9eb0bd55 377 #if defined(TVMET_HAVE_LONG_DOUBLE)
narshu 1:cc2a9eb0bd55 378 TVMET_IMPLEMENT_MACRO(atan2, long double)
narshu 1:cc2a9eb0bd55 379 TVMET_IMPLEMENT_MACRO(drem, long double)
narshu 1:cc2a9eb0bd55 380 TVMET_IMPLEMENT_MACRO(fmod, long double)
narshu 1:cc2a9eb0bd55 381 TVMET_IMPLEMENT_MACRO(hypot, long double)
narshu 1:cc2a9eb0bd55 382 TVMET_IMPLEMENT_MACRO(jn, long double)
narshu 1:cc2a9eb0bd55 383 TVMET_IMPLEMENT_MACRO(yn, long double)
narshu 1:cc2a9eb0bd55 384 TVMET_IMPLEMENT_MACRO(pow, long double)
narshu 1:cc2a9eb0bd55 385 #endif // defined(TVMET_HAVE_LONG_DOUBLE)
narshu 1:cc2a9eb0bd55 386
narshu 1:cc2a9eb0bd55 387 #undef TVMET_IMPLEMENT_MACRO
narshu 1:cc2a9eb0bd55 388
narshu 1:cc2a9eb0bd55 389
narshu 1:cc2a9eb0bd55 390 /*
narshu 1:cc2a9eb0bd55 391 * complex math
narshu 1:cc2a9eb0bd55 392 */
narshu 1:cc2a9eb0bd55 393
narshu 1:cc2a9eb0bd55 394 #if defined(TVMET_HAVE_COMPLEX) && defined(TVMET_HAVE_COMPLEX_MATH1)
narshu 1:cc2a9eb0bd55 395 /**
narshu 1:cc2a9eb0bd55 396 * \fn pow(const Matrix<T, Rows, Cols>& lhs, const std::complex<T>& rhs)
narshu 1:cc2a9eb0bd55 397 * \ingroup _binary_function
narshu 1:cc2a9eb0bd55 398 */
narshu 1:cc2a9eb0bd55 399 template<class T, std::size_t Rows, std::size_t Cols>
narshu 1:cc2a9eb0bd55 400 inline
narshu 1:cc2a9eb0bd55 401 XprMatrix<
narshu 1:cc2a9eb0bd55 402 XprBinOp<
narshu 1:cc2a9eb0bd55 403 Fcnl_pow<T, std::complex<T> >,
narshu 1:cc2a9eb0bd55 404 MatrixConstReference<T, Rows, Cols>,
narshu 1:cc2a9eb0bd55 405 XprLiteral< std::complex<T> >
narshu 1:cc2a9eb0bd55 406 >,
narshu 1:cc2a9eb0bd55 407 Rows, Cols
narshu 1:cc2a9eb0bd55 408 >
narshu 1:cc2a9eb0bd55 409 pow(const Matrix<T, Rows, Cols>& lhs, const std::complex<T>& rhs) {
narshu 1:cc2a9eb0bd55 410 typedef XprBinOp<
narshu 1:cc2a9eb0bd55 411 Fcnl_pow<T, std::complex<T> >,
narshu 1:cc2a9eb0bd55 412 MatrixConstReference<T, Rows, Cols>,
narshu 1:cc2a9eb0bd55 413 XprLiteral< std::complex<T> >
narshu 1:cc2a9eb0bd55 414 > expr_type;
narshu 1:cc2a9eb0bd55 415 return XprMatrix<expr_type, Rows, Cols>(
narshu 1:cc2a9eb0bd55 416 expr_type(lhs.const_ref(), XprLiteral< std::complex<T> >(rhs)));
narshu 1:cc2a9eb0bd55 417 }
narshu 1:cc2a9eb0bd55 418
narshu 1:cc2a9eb0bd55 419
narshu 1:cc2a9eb0bd55 420 /**
narshu 1:cc2a9eb0bd55 421 * \fn pow(const Matrix<std::complex<T>, Rows, Cols>& lhs, const std::complex<T>& rhs)
narshu 1:cc2a9eb0bd55 422 * \ingroup _binary_function
narshu 1:cc2a9eb0bd55 423 */
narshu 1:cc2a9eb0bd55 424 template<class T, std::size_t Rows, std::size_t Cols>
narshu 1:cc2a9eb0bd55 425 inline
narshu 1:cc2a9eb0bd55 426 XprMatrix<
narshu 1:cc2a9eb0bd55 427 XprBinOp<
narshu 1:cc2a9eb0bd55 428 Fcnl_pow< std::complex<T>, std::complex<T> >,
narshu 1:cc2a9eb0bd55 429 MatrixConstReference<std::complex<T>, Rows, Cols>,
narshu 1:cc2a9eb0bd55 430 XprLiteral< std::complex<T> >
narshu 1:cc2a9eb0bd55 431 >,
narshu 1:cc2a9eb0bd55 432 Rows, Cols
narshu 1:cc2a9eb0bd55 433 >
narshu 1:cc2a9eb0bd55 434 pow(const Matrix<std::complex<T>, Rows, Cols>& lhs, const std::complex<T>& rhs) {
narshu 1:cc2a9eb0bd55 435 typedef XprBinOp<
narshu 1:cc2a9eb0bd55 436 Fcnl_pow< std::complex<T>, std::complex<T> >,
narshu 1:cc2a9eb0bd55 437 MatrixConstReference<std::complex<T>, Rows, Cols>,
narshu 1:cc2a9eb0bd55 438 XprLiteral< std::complex<T> >
narshu 1:cc2a9eb0bd55 439 > expr_type;
narshu 1:cc2a9eb0bd55 440 return XprMatrix<expr_type, Rows, Cols>(
narshu 1:cc2a9eb0bd55 441 expr_type(lhs.const_ref(), XprLiteral< std::complex<T> >(rhs)));
narshu 1:cc2a9eb0bd55 442 }
narshu 1:cc2a9eb0bd55 443
narshu 1:cc2a9eb0bd55 444
narshu 1:cc2a9eb0bd55 445 /**
narshu 1:cc2a9eb0bd55 446 * \fn pow(const Matrix<std::complex<T>, Rows, Cols>& lhs, const T& rhs)
narshu 1:cc2a9eb0bd55 447 * \ingroup _binary_function
narshu 1:cc2a9eb0bd55 448 */
narshu 1:cc2a9eb0bd55 449 template<class T, std::size_t Rows, std::size_t Cols>
narshu 1:cc2a9eb0bd55 450 inline
narshu 1:cc2a9eb0bd55 451 XprMatrix<
narshu 1:cc2a9eb0bd55 452 XprBinOp<
narshu 1:cc2a9eb0bd55 453 Fcnl_pow<std::complex<T>, T>,
narshu 1:cc2a9eb0bd55 454 MatrixConstReference<std::complex<T>, Rows, Cols>,
narshu 1:cc2a9eb0bd55 455 XprLiteral<T>
narshu 1:cc2a9eb0bd55 456 >,
narshu 1:cc2a9eb0bd55 457 Rows, Cols
narshu 1:cc2a9eb0bd55 458 >
narshu 1:cc2a9eb0bd55 459 pow(const Matrix<std::complex<T>, Rows, Cols>& lhs, const T& rhs) {
narshu 1:cc2a9eb0bd55 460 typedef XprBinOp<
narshu 1:cc2a9eb0bd55 461 Fcnl_pow<std::complex<T>, T>,
narshu 1:cc2a9eb0bd55 462 MatrixConstReference<std::complex<T>, Rows, Cols>,
narshu 1:cc2a9eb0bd55 463 XprLiteral<T>
narshu 1:cc2a9eb0bd55 464 > expr_type;
narshu 1:cc2a9eb0bd55 465 return XprMatrix<expr_type, Rows, Cols>(
narshu 1:cc2a9eb0bd55 466 expr_type(lhs.const_ref(), XprLiteral<T>(rhs)));
narshu 1:cc2a9eb0bd55 467 }
narshu 1:cc2a9eb0bd55 468
narshu 1:cc2a9eb0bd55 469
narshu 1:cc2a9eb0bd55 470 /**
narshu 1:cc2a9eb0bd55 471 * \fn pow(const Matrix<std::complex<T>, Rows, Cols>& lhs, int rhs)
narshu 1:cc2a9eb0bd55 472 * \ingroup _binary_function
narshu 1:cc2a9eb0bd55 473 */
narshu 1:cc2a9eb0bd55 474 template<class T, std::size_t Rows, std::size_t Cols>
narshu 1:cc2a9eb0bd55 475 inline
narshu 1:cc2a9eb0bd55 476 XprMatrix<
narshu 1:cc2a9eb0bd55 477 XprBinOp<
narshu 1:cc2a9eb0bd55 478 Fcnl_pow<std::complex<T>, int>,
narshu 1:cc2a9eb0bd55 479 MatrixConstReference<std::complex<T>, Rows, Cols>,
narshu 1:cc2a9eb0bd55 480 XprLiteral<int>
narshu 1:cc2a9eb0bd55 481 >,
narshu 1:cc2a9eb0bd55 482 Rows, Cols
narshu 1:cc2a9eb0bd55 483 >
narshu 1:cc2a9eb0bd55 484 pow(const Matrix<std::complex<T>, Rows, Cols>& lhs, int rhs) {
narshu 1:cc2a9eb0bd55 485 typedef XprBinOp<
narshu 1:cc2a9eb0bd55 486 Fcnl_pow<std::complex<T>, int>,
narshu 1:cc2a9eb0bd55 487 MatrixConstReference<std::complex<T>, Rows, Cols>,
narshu 1:cc2a9eb0bd55 488 XprLiteral<int>
narshu 1:cc2a9eb0bd55 489 > expr_type;
narshu 1:cc2a9eb0bd55 490 return XprMatrix<expr_type, Rows, Cols>(
narshu 1:cc2a9eb0bd55 491 expr_type(lhs.const_ref(), XprLiteral<int>(rhs)));
narshu 1:cc2a9eb0bd55 492 }
narshu 1:cc2a9eb0bd55 493
narshu 1:cc2a9eb0bd55 494
narshu 1:cc2a9eb0bd55 495 /**
narshu 1:cc2a9eb0bd55 496 * \fn polar(const Matrix<T, Rows, Cols>& lhs, const T& rhs)
narshu 1:cc2a9eb0bd55 497 * \ingroup _binary_function
narshu 1:cc2a9eb0bd55 498 */
narshu 1:cc2a9eb0bd55 499 template<class T, std::size_t Rows, std::size_t Cols>
narshu 1:cc2a9eb0bd55 500 inline
narshu 1:cc2a9eb0bd55 501 XprMatrix<
narshu 1:cc2a9eb0bd55 502 XprBinOp<
narshu 1:cc2a9eb0bd55 503 Fcnl_polar<T, T>,
narshu 1:cc2a9eb0bd55 504 MatrixConstReference<T, Rows, Cols>,
narshu 1:cc2a9eb0bd55 505 XprLiteral<T>
narshu 1:cc2a9eb0bd55 506 >,
narshu 1:cc2a9eb0bd55 507 Rows, Cols
narshu 1:cc2a9eb0bd55 508 >
narshu 1:cc2a9eb0bd55 509 polar(const Matrix<T, Rows, Cols>& lhs, const T& rhs) {
narshu 1:cc2a9eb0bd55 510 typedef XprBinOp<
narshu 1:cc2a9eb0bd55 511 Fcnl_polar<T, T>,
narshu 1:cc2a9eb0bd55 512 MatrixConstReference<T, Rows, Cols>,
narshu 1:cc2a9eb0bd55 513 XprLiteral<T>
narshu 1:cc2a9eb0bd55 514 > expr_type;
narshu 1:cc2a9eb0bd55 515 return XprMatrix<expr_type, Rows, Cols>(
narshu 1:cc2a9eb0bd55 516 expr_type(lhs.const_ref(), XprLiteral<T>(rhs)));
narshu 1:cc2a9eb0bd55 517 }
narshu 1:cc2a9eb0bd55 518
narshu 1:cc2a9eb0bd55 519 #endif // defined(TVMET_HAVE_COMPLEX) && defined(TVMET_HAVE_COMPLEX_MATH1)
narshu 1:cc2a9eb0bd55 520
narshu 1:cc2a9eb0bd55 521 #if defined(TVMET_HAVE_COMPLEX) && defined(TVMET_HAVE_COMPLEX_MATH2)
narshu 1:cc2a9eb0bd55 522 // to be written (atan2)
narshu 1:cc2a9eb0bd55 523 #endif // defined(TVMET_HAVE_COMPLEX) && defined(TVMET_HAVE_COMPLEX_MATH2)
narshu 1:cc2a9eb0bd55 524
narshu 1:cc2a9eb0bd55 525
narshu 1:cc2a9eb0bd55 526 } // namespace tvmet
narshu 1:cc2a9eb0bd55 527
narshu 1:cc2a9eb0bd55 528 #endif // TVMET_MATRIX_BINARY_FUNCTIONS_H
narshu 1:cc2a9eb0bd55 529
narshu 1:cc2a9eb0bd55 530 // Local Variables:
narshu 1:cc2a9eb0bd55 531 // mode:C++
narshu 1:cc2a9eb0bd55 532 // tab-width:8
narshu 1:cc2a9eb0bd55 533 // End: