This is the Tiny Vector Matrix Expression Templates library found at http://tvmet.sourceforge.net. It is the fastest and most compact matrix lib out there (for < 10x10 matricies). I have done some minor tweaks to make it compile for mbed. For examples and hints on how to use, see: http://tvmet.sourceforge.net/usage.html

Dependents:   Eurobot_2012_Secondary

Committer:
madcowswe
Date:
Wed Mar 28 15:53:45 2012 +0000
Revision:
0:feb4117d16d8

        

Who changed what in which revision?

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