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: MatrixUnaryFunctions.h,v 1.13 2007-06-23 15:58:58 opetzold Exp $
madcowswe 0:feb4117d16d8 22 */
madcowswe 0:feb4117d16d8 23
madcowswe 0:feb4117d16d8 24 #ifndef TVMET_MATRIX_UNARY_FUNCTIONS_H
madcowswe 0:feb4117d16d8 25 #define TVMET_MATRIX_UNARY_FUNCTIONS_H
madcowswe 0:feb4117d16d8 26
madcowswe 0:feb4117d16d8 27 namespace tvmet {
madcowswe 0:feb4117d16d8 28
madcowswe 0:feb4117d16d8 29
madcowswe 0:feb4117d16d8 30 /*********************************************************
madcowswe 0:feb4117d16d8 31 * PART I: DECLARATION
madcowswe 0:feb4117d16d8 32 *********************************************************/
madcowswe 0:feb4117d16d8 33
madcowswe 0:feb4117d16d8 34 /*
madcowswe 0:feb4117d16d8 35 * unary_function(Matrix<T, Rows, Cols>)
madcowswe 0:feb4117d16d8 36 */
madcowswe 0:feb4117d16d8 37 #define TVMET_DECLARE_MACRO(NAME) \
madcowswe 0:feb4117d16d8 38 template<class T, std::size_t Rows, std::size_t Cols> \
madcowswe 0:feb4117d16d8 39 inline \
madcowswe 0:feb4117d16d8 40 XprMatrix< \
madcowswe 0:feb4117d16d8 41 XprUnOp< \
madcowswe 0:feb4117d16d8 42 Fcnl_##NAME<T>, \
madcowswe 0:feb4117d16d8 43 MatrixConstReference<T, Rows, Cols> \
madcowswe 0:feb4117d16d8 44 >, \
madcowswe 0:feb4117d16d8 45 Rows, Cols \
madcowswe 0:feb4117d16d8 46 > \
madcowswe 0:feb4117d16d8 47 NAME(const Matrix<T, Rows, Cols>& rhs) TVMET_CXX_ALWAYS_INLINE;
madcowswe 0:feb4117d16d8 48
madcowswe 0:feb4117d16d8 49 TVMET_DECLARE_MACRO(abs)
madcowswe 0:feb4117d16d8 50 TVMET_DECLARE_MACRO(cbrt)
madcowswe 0:feb4117d16d8 51 TVMET_DECLARE_MACRO(ceil)
madcowswe 0:feb4117d16d8 52 TVMET_DECLARE_MACRO(floor)
madcowswe 0:feb4117d16d8 53 TVMET_DECLARE_MACRO(rint)
madcowswe 0:feb4117d16d8 54 TVMET_DECLARE_MACRO(sin)
madcowswe 0:feb4117d16d8 55 TVMET_DECLARE_MACRO(cos)
madcowswe 0:feb4117d16d8 56 TVMET_DECLARE_MACRO(tan)
madcowswe 0:feb4117d16d8 57 TVMET_DECLARE_MACRO(sinh)
madcowswe 0:feb4117d16d8 58 TVMET_DECLARE_MACRO(cosh)
madcowswe 0:feb4117d16d8 59 TVMET_DECLARE_MACRO(tanh)
madcowswe 0:feb4117d16d8 60 TVMET_DECLARE_MACRO(asin)
madcowswe 0:feb4117d16d8 61 TVMET_DECLARE_MACRO(acos)
madcowswe 0:feb4117d16d8 62 TVMET_DECLARE_MACRO(atan)
madcowswe 0:feb4117d16d8 63 TVMET_DECLARE_MACRO(exp)
madcowswe 0:feb4117d16d8 64 TVMET_DECLARE_MACRO(log)
madcowswe 0:feb4117d16d8 65 TVMET_DECLARE_MACRO(log10)
madcowswe 0:feb4117d16d8 66 TVMET_DECLARE_MACRO(sqrt)
madcowswe 0:feb4117d16d8 67
madcowswe 0:feb4117d16d8 68 #if defined(TVMET_HAVE_IEEE_MATH)
madcowswe 0:feb4117d16d8 69 TVMET_DECLARE_MACRO(asinh)
madcowswe 0:feb4117d16d8 70 TVMET_DECLARE_MACRO(acosh)
madcowswe 0:feb4117d16d8 71 TVMET_DECLARE_MACRO(atanh)
madcowswe 0:feb4117d16d8 72 TVMET_DECLARE_MACRO(expm1)
madcowswe 0:feb4117d16d8 73 TVMET_DECLARE_MACRO(log1p)
madcowswe 0:feb4117d16d8 74 TVMET_DECLARE_MACRO(erf)
madcowswe 0:feb4117d16d8 75 TVMET_DECLARE_MACRO(erfc)
madcowswe 0:feb4117d16d8 76 TVMET_DECLARE_MACRO(j0)
madcowswe 0:feb4117d16d8 77 TVMET_DECLARE_MACRO(j1)
madcowswe 0:feb4117d16d8 78 TVMET_DECLARE_MACRO(y0)
madcowswe 0:feb4117d16d8 79 TVMET_DECLARE_MACRO(y1)
madcowswe 0:feb4117d16d8 80 TVMET_DECLARE_MACRO(lgamma)
madcowswe 0:feb4117d16d8 81 /** \todo isnan etc. - default return is only an int! */
madcowswe 0:feb4117d16d8 82
madcowswe 0:feb4117d16d8 83 TVMET_DECLARE_MACRO(finite)
madcowswe 0:feb4117d16d8 84 #endif // defined(TVMET_HAVE_IEEE_MATH)
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 * unary_function(Matrix<std::complex<T>, Rows, Cols>)
madcowswe 0:feb4117d16d8 91 */
madcowswe 0:feb4117d16d8 92 #if defined(TVMET_HAVE_COMPLEX)
madcowswe 0:feb4117d16d8 93 #define TVMET_DECLARE_MACRO(NAME) \
madcowswe 0:feb4117d16d8 94 template<class T, std::size_t Rows, std::size_t Cols> \
madcowswe 0:feb4117d16d8 95 inline \
madcowswe 0:feb4117d16d8 96 XprMatrix< \
madcowswe 0:feb4117d16d8 97 XprUnOp< \
madcowswe 0:feb4117d16d8 98 Fcnl_##NAME< std::complex<T> >, \
madcowswe 0:feb4117d16d8 99 MatrixConstReference<std::complex<T>, Rows, Cols> \
madcowswe 0:feb4117d16d8 100 >, \
madcowswe 0:feb4117d16d8 101 Rows, Cols \
madcowswe 0:feb4117d16d8 102 > \
madcowswe 0:feb4117d16d8 103 NAME(const Matrix<std::complex<T>, Rows, Cols>& rhs) TVMET_CXX_ALWAYS_INLINE;
madcowswe 0:feb4117d16d8 104
madcowswe 0:feb4117d16d8 105 TVMET_DECLARE_MACRO(real)
madcowswe 0:feb4117d16d8 106 TVMET_DECLARE_MACRO(imag)
madcowswe 0:feb4117d16d8 107 TVMET_DECLARE_MACRO(arg)
madcowswe 0:feb4117d16d8 108 TVMET_DECLARE_MACRO(norm)
madcowswe 0:feb4117d16d8 109 TVMET_DECLARE_MACRO(conj)
madcowswe 0:feb4117d16d8 110
madcowswe 0:feb4117d16d8 111 #undef TVMET_DECLARE_MACRO
madcowswe 0:feb4117d16d8 112
madcowswe 0:feb4117d16d8 113 #endif // defined(TVMET_HAVE_COMPLEX)
madcowswe 0:feb4117d16d8 114
madcowswe 0:feb4117d16d8 115
madcowswe 0:feb4117d16d8 116 /*********************************************************
madcowswe 0:feb4117d16d8 117 * PART II: IMPLEMENTATION
madcowswe 0:feb4117d16d8 118 *********************************************************/
madcowswe 0:feb4117d16d8 119
madcowswe 0:feb4117d16d8 120
madcowswe 0:feb4117d16d8 121 /*
madcowswe 0:feb4117d16d8 122 * unary_function(Matrix<T, Rows, Cols>)
madcowswe 0:feb4117d16d8 123 */
madcowswe 0:feb4117d16d8 124 #define TVMET_IMPLEMENT_MACRO(NAME) \
madcowswe 0:feb4117d16d8 125 template<class T, std::size_t Rows, std::size_t Cols> \
madcowswe 0:feb4117d16d8 126 inline \
madcowswe 0:feb4117d16d8 127 XprMatrix< \
madcowswe 0:feb4117d16d8 128 XprUnOp< \
madcowswe 0:feb4117d16d8 129 Fcnl_##NAME<T>, \
madcowswe 0:feb4117d16d8 130 MatrixConstReference<T, Rows, Cols> \
madcowswe 0:feb4117d16d8 131 >, \
madcowswe 0:feb4117d16d8 132 Rows, Cols \
madcowswe 0:feb4117d16d8 133 > \
madcowswe 0:feb4117d16d8 134 NAME(const Matrix<T, Rows, Cols>& rhs) { \
madcowswe 0:feb4117d16d8 135 typedef XprUnOp< \
madcowswe 0:feb4117d16d8 136 Fcnl_##NAME<T>, \
madcowswe 0:feb4117d16d8 137 MatrixConstReference<T, Rows, Cols> \
madcowswe 0:feb4117d16d8 138 > expr_type; \
madcowswe 0:feb4117d16d8 139 return XprMatrix<expr_type, Rows, Cols>(expr_type(rhs.const_ref())); \
madcowswe 0:feb4117d16d8 140 }
madcowswe 0:feb4117d16d8 141
madcowswe 0:feb4117d16d8 142 TVMET_IMPLEMENT_MACRO(abs)
madcowswe 0:feb4117d16d8 143 TVMET_IMPLEMENT_MACRO(cbrt)
madcowswe 0:feb4117d16d8 144 TVMET_IMPLEMENT_MACRO(ceil)
madcowswe 0:feb4117d16d8 145 TVMET_IMPLEMENT_MACRO(floor)
madcowswe 0:feb4117d16d8 146 TVMET_IMPLEMENT_MACRO(rint)
madcowswe 0:feb4117d16d8 147 TVMET_IMPLEMENT_MACRO(sin)
madcowswe 0:feb4117d16d8 148 TVMET_IMPLEMENT_MACRO(cos)
madcowswe 0:feb4117d16d8 149 TVMET_IMPLEMENT_MACRO(tan)
madcowswe 0:feb4117d16d8 150 TVMET_IMPLEMENT_MACRO(sinh)
madcowswe 0:feb4117d16d8 151 TVMET_IMPLEMENT_MACRO(cosh)
madcowswe 0:feb4117d16d8 152 TVMET_IMPLEMENT_MACRO(tanh)
madcowswe 0:feb4117d16d8 153 TVMET_IMPLEMENT_MACRO(asin)
madcowswe 0:feb4117d16d8 154 TVMET_IMPLEMENT_MACRO(acos)
madcowswe 0:feb4117d16d8 155 TVMET_IMPLEMENT_MACRO(atan)
madcowswe 0:feb4117d16d8 156 TVMET_IMPLEMENT_MACRO(exp)
madcowswe 0:feb4117d16d8 157 TVMET_IMPLEMENT_MACRO(log)
madcowswe 0:feb4117d16d8 158 TVMET_IMPLEMENT_MACRO(log10)
madcowswe 0:feb4117d16d8 159 TVMET_IMPLEMENT_MACRO(sqrt)
madcowswe 0:feb4117d16d8 160
madcowswe 0:feb4117d16d8 161 #if defined(TVMET_HAVE_IEEE_MATH)
madcowswe 0:feb4117d16d8 162 TVMET_IMPLEMENT_MACRO(asinh)
madcowswe 0:feb4117d16d8 163 TVMET_IMPLEMENT_MACRO(acosh)
madcowswe 0:feb4117d16d8 164 TVMET_IMPLEMENT_MACRO(atanh)
madcowswe 0:feb4117d16d8 165 TVMET_IMPLEMENT_MACRO(expm1)
madcowswe 0:feb4117d16d8 166 TVMET_IMPLEMENT_MACRO(log1p)
madcowswe 0:feb4117d16d8 167 TVMET_IMPLEMENT_MACRO(erf)
madcowswe 0:feb4117d16d8 168 TVMET_IMPLEMENT_MACRO(erfc)
madcowswe 0:feb4117d16d8 169 TVMET_IMPLEMENT_MACRO(j0)
madcowswe 0:feb4117d16d8 170 TVMET_IMPLEMENT_MACRO(j1)
madcowswe 0:feb4117d16d8 171 TVMET_IMPLEMENT_MACRO(y0)
madcowswe 0:feb4117d16d8 172 TVMET_IMPLEMENT_MACRO(y1)
madcowswe 0:feb4117d16d8 173 TVMET_IMPLEMENT_MACRO(lgamma)
madcowswe 0:feb4117d16d8 174 /** \todo isnan etc. - default return is only an int! */
madcowswe 0:feb4117d16d8 175
madcowswe 0:feb4117d16d8 176 TVMET_IMPLEMENT_MACRO(finite)
madcowswe 0:feb4117d16d8 177 #endif // defined(TVMET_HAVE_IEEE_MATH)
madcowswe 0:feb4117d16d8 178
madcowswe 0:feb4117d16d8 179 #undef TVMET_IMPLEMENT_MACRO
madcowswe 0:feb4117d16d8 180
madcowswe 0:feb4117d16d8 181
madcowswe 0:feb4117d16d8 182 /*
madcowswe 0:feb4117d16d8 183 * unary_function(Matrix<std::complex<T>, Rows, Cols>)
madcowswe 0:feb4117d16d8 184 */
madcowswe 0:feb4117d16d8 185 #if defined(TVMET_HAVE_COMPLEX)
madcowswe 0:feb4117d16d8 186 #define TVMET_IMPLEMENT_MACRO(NAME) \
madcowswe 0:feb4117d16d8 187 template<class T, std::size_t Rows, std::size_t Cols> \
madcowswe 0:feb4117d16d8 188 inline \
madcowswe 0:feb4117d16d8 189 XprMatrix< \
madcowswe 0:feb4117d16d8 190 XprUnOp< \
madcowswe 0:feb4117d16d8 191 Fcnl_##NAME< std::complex<T> >, \
madcowswe 0:feb4117d16d8 192 MatrixConstReference<std::complex<T>, Rows, Cols> \
madcowswe 0:feb4117d16d8 193 >, \
madcowswe 0:feb4117d16d8 194 Rows, Cols \
madcowswe 0:feb4117d16d8 195 > \
madcowswe 0:feb4117d16d8 196 NAME(const Matrix<std::complex<T>, Rows, Cols>& rhs) { \
madcowswe 0:feb4117d16d8 197 typedef XprUnOp< \
madcowswe 0:feb4117d16d8 198 Fcnl_##NAME< std::complex<T> >, \
madcowswe 0:feb4117d16d8 199 MatrixConstReference<std::complex<T>, Rows, Cols> \
madcowswe 0:feb4117d16d8 200 > expr_type; \
madcowswe 0:feb4117d16d8 201 return XprMatrix<expr_type, Rows, Cols>(expr_type(rhs.const_ref())); \
madcowswe 0:feb4117d16d8 202 }
madcowswe 0:feb4117d16d8 203
madcowswe 0:feb4117d16d8 204 TVMET_IMPLEMENT_MACRO(real)
madcowswe 0:feb4117d16d8 205 TVMET_IMPLEMENT_MACRO(imag)
madcowswe 0:feb4117d16d8 206 TVMET_IMPLEMENT_MACRO(arg)
madcowswe 0:feb4117d16d8 207 TVMET_IMPLEMENT_MACRO(norm)
madcowswe 0:feb4117d16d8 208 TVMET_IMPLEMENT_MACRO(conj)
madcowswe 0:feb4117d16d8 209
madcowswe 0:feb4117d16d8 210 #undef TVMET_IMPLEMENT_MACRO
madcowswe 0:feb4117d16d8 211
madcowswe 0:feb4117d16d8 212 #endif // defined(TVMET_HAVE_COMPLEX)
madcowswe 0:feb4117d16d8 213
madcowswe 0:feb4117d16d8 214
madcowswe 0:feb4117d16d8 215 } // namespace tvmet
madcowswe 0:feb4117d16d8 216
madcowswe 0:feb4117d16d8 217 #endif // TVMET_MATRIX_UNARY_FUNCTIONS_H
madcowswe 0:feb4117d16d8 218
madcowswe 0:feb4117d16d8 219 // Local Variables:
madcowswe 0:feb4117d16d8 220 // mode:C++
madcowswe 0:feb4117d16d8 221 // tab-width:8
madcowswe 0:feb4117d16d8 222 // End: