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: CommaInitializer.h,v 1.18 2007-06-23 15:58:58 opetzold Exp $
madcowswe 0:feb4117d16d8 22 */
madcowswe 0:feb4117d16d8 23
madcowswe 0:feb4117d16d8 24 #ifndef TVMET_COMMA_INITIALIZER_H
madcowswe 0:feb4117d16d8 25 #define TVMET_COMMA_INITIALIZER_H
madcowswe 0:feb4117d16d8 26
madcowswe 0:feb4117d16d8 27 #include <tvmet/CompileTimeError.h>
madcowswe 0:feb4117d16d8 28
madcowswe 0:feb4117d16d8 29 namespace tvmet {
madcowswe 0:feb4117d16d8 30
madcowswe 0:feb4117d16d8 31
madcowswe 0:feb4117d16d8 32 /**
madcowswe 0:feb4117d16d8 33 * \class CommaInitializer CommaInitializer.h "tvmet/CommaInitializer.h"
madcowswe 0:feb4117d16d8 34 * \brief Initialize classes using a comma separated lists.
madcowswe 0:feb4117d16d8 35 *
madcowswe 0:feb4117d16d8 36 * The comma operator is called when it appears next to an object of
madcowswe 0:feb4117d16d8 37 * the type the comma is defined for. However, "operator," is not called
madcowswe 0:feb4117d16d8 38 * for function argument lists, only for objects that are out in the open,
madcowswe 0:feb4117d16d8 39 * separated by commas (Thinking C++
madcowswe 0:feb4117d16d8 40 * <a href=http://www.ida.liu.se/~TDDA14/online/v1ticpp/Chapter12.html>
madcowswe 0:feb4117d16d8 41 * Ch.12: Operator comma</a>).
madcowswe 0:feb4117d16d8 42 *
madcowswe 0:feb4117d16d8 43 * This implementation uses the same technique as described in Todd Veldhuizen
madcowswe 0:feb4117d16d8 44 * Techniques for Scientific C++
madcowswe 0:feb4117d16d8 45 * <a href=http://extreme.indiana.edu/~tveldhui/papers/techniques/techniques01.html#l43>
madcowswe 0:feb4117d16d8 46 * chapter 1.11 Comma overloading</a>.
madcowswe 0:feb4117d16d8 47 *
madcowswe 0:feb4117d16d8 48 * The initializer list is avaible after instanciation of the object,
madcowswe 0:feb4117d16d8 49 * therefore use it like:
madcowswe 0:feb4117d16d8 50 * \code
madcowswe 0:feb4117d16d8 51 * vector3d t;
madcowswe 0:feb4117d16d8 52 * t = 1.0, 2.0, 3.0;
madcowswe 0:feb4117d16d8 53 * \endcode
madcowswe 0:feb4117d16d8 54 * It's evaluated to (((t = 1.0), 2.0), 3.0)
madcowswe 0:feb4117d16d8 55 *
madcowswe 0:feb4117d16d8 56 * For matrizes the initilization is done row wise.
madcowswe 0:feb4117d16d8 57 *
madcowswe 0:feb4117d16d8 58 * If the comma separted list of values longer then the size of the vector
madcowswe 0:feb4117d16d8 59 * or matrix a compile time error will occour. Otherwise the pending values
madcowswe 0:feb4117d16d8 60 * will be written random into the memory.
madcowswe 0:feb4117d16d8 61 *
madcowswe 0:feb4117d16d8 62 */
madcowswe 0:feb4117d16d8 63 template<class Obj, std::size_t LEN>
madcowswe 0:feb4117d16d8 64 class CommaInitializer
madcowswe 0:feb4117d16d8 65 {
madcowswe 0:feb4117d16d8 66 CommaInitializer();
madcowswe 0:feb4117d16d8 67 CommaInitializer& operator=(const CommaInitializer&);
madcowswe 0:feb4117d16d8 68
madcowswe 0:feb4117d16d8 69 private:
madcowswe 0:feb4117d16d8 70 /**
madcowswe 0:feb4117d16d8 71 * \class Initializer
madcowswe 0:feb4117d16d8 72 * \brief Helper fo recursive overloaded comma operator.
madcowswe 0:feb4117d16d8 73 */
madcowswe 0:feb4117d16d8 74 template<class T, std::size_t N> class Initializer
madcowswe 0:feb4117d16d8 75 {
madcowswe 0:feb4117d16d8 76 Initializer();
madcowswe 0:feb4117d16d8 77 Initializer& operator=(const Initializer&);
madcowswe 0:feb4117d16d8 78
madcowswe 0:feb4117d16d8 79 public:
madcowswe 0:feb4117d16d8 80 typedef T value_type;
madcowswe 0:feb4117d16d8 81 typedef T* iterator;
madcowswe 0:feb4117d16d8 82
madcowswe 0:feb4117d16d8 83 public:
madcowswe 0:feb4117d16d8 84 Initializer(iterator iter) : m_iter(iter) { }
madcowswe 0:feb4117d16d8 85
madcowswe 0:feb4117d16d8 86 /** Overloads the comma operator for recursive assign values from comma
madcowswe 0:feb4117d16d8 87 separated list. */
madcowswe 0:feb4117d16d8 88 Initializer<value_type, N+1> operator,(value_type rhs)
madcowswe 0:feb4117d16d8 89 {
madcowswe 0:feb4117d16d8 90 TVMET_CT_CONDITION(N < LEN, CommaInitializerList_is_too_long)
madcowswe 0:feb4117d16d8 91 *m_iter = rhs;
madcowswe 0:feb4117d16d8 92 return Initializer<value_type, N+1>(m_iter + 1);
madcowswe 0:feb4117d16d8 93 }
madcowswe 0:feb4117d16d8 94
madcowswe 0:feb4117d16d8 95 private:
madcowswe 0:feb4117d16d8 96 iterator m_iter;
madcowswe 0:feb4117d16d8 97 };
madcowswe 0:feb4117d16d8 98
madcowswe 0:feb4117d16d8 99 public:
madcowswe 0:feb4117d16d8 100 typedef typename Obj::value_type value_type;
madcowswe 0:feb4117d16d8 101 typedef value_type* iterator;
madcowswe 0:feb4117d16d8 102
madcowswe 0:feb4117d16d8 103 public:
madcowswe 0:feb4117d16d8 104 CommaInitializer(const CommaInitializer& rhs)
madcowswe 0:feb4117d16d8 105 : m_object(rhs.m_object),
madcowswe 0:feb4117d16d8 106 m_data(rhs.m_data),
madcowswe 0:feb4117d16d8 107 m_wipeout_on_destruct(true)
madcowswe 0:feb4117d16d8 108 {
madcowswe 0:feb4117d16d8 109 rhs.disable();
madcowswe 0:feb4117d16d8 110 }
madcowswe 0:feb4117d16d8 111
madcowswe 0:feb4117d16d8 112 /** Constructor used by Vector or Matrix operator(value_type rhs) */
madcowswe 0:feb4117d16d8 113 CommaInitializer(Obj& obj, value_type x)
madcowswe 0:feb4117d16d8 114 : m_object(obj),
madcowswe 0:feb4117d16d8 115 m_data(x),
madcowswe 0:feb4117d16d8 116 m_wipeout_on_destruct(true)
madcowswe 0:feb4117d16d8 117 { }
madcowswe 0:feb4117d16d8 118
madcowswe 0:feb4117d16d8 119 /** Destructs and assigns the comma separated value. */
madcowswe 0:feb4117d16d8 120 ~CommaInitializer() {
madcowswe 0:feb4117d16d8 121 if(m_wipeout_on_destruct) m_object.assign_value(m_data);
madcowswe 0:feb4117d16d8 122 }
madcowswe 0:feb4117d16d8 123
madcowswe 0:feb4117d16d8 124 /** Overloaded comma operator, called only once for the first occoured comma. This
madcowswe 0:feb4117d16d8 125 means the first value is assigned by %operator=() and the 2nd value after the
madcowswe 0:feb4117d16d8 126 comma. Therfore we call the %Initializer::operator,() for the list starting
madcowswe 0:feb4117d16d8 127 after the 2nd. */
madcowswe 0:feb4117d16d8 128 Initializer<value_type, 2> operator,(value_type rhs);
madcowswe 0:feb4117d16d8 129
madcowswe 0:feb4117d16d8 130 void disable() const { m_wipeout_on_destruct = false; }
madcowswe 0:feb4117d16d8 131
madcowswe 0:feb4117d16d8 132 private:
madcowswe 0:feb4117d16d8 133 Obj& m_object;
madcowswe 0:feb4117d16d8 134 value_type m_data;
madcowswe 0:feb4117d16d8 135 mutable bool m_wipeout_on_destruct;
madcowswe 0:feb4117d16d8 136 };
madcowswe 0:feb4117d16d8 137
madcowswe 0:feb4117d16d8 138
madcowswe 0:feb4117d16d8 139 /*
madcowswe 0:feb4117d16d8 140 * Implementation
madcowswe 0:feb4117d16d8 141 */
madcowswe 0:feb4117d16d8 142 template<class Obj, std::size_t LEN>
madcowswe 0:feb4117d16d8 143 typename CommaInitializer<Obj, LEN>::template Initializer<typename Obj::value_type, 2>
madcowswe 0:feb4117d16d8 144 CommaInitializer<Obj, LEN>::operator,(typename Obj::value_type rhs)
madcowswe 0:feb4117d16d8 145 {
madcowswe 0:feb4117d16d8 146 m_wipeout_on_destruct = false;
madcowswe 0:feb4117d16d8 147 iterator iter1 = m_object.data();
madcowswe 0:feb4117d16d8 148 *iter1 = m_data;
madcowswe 0:feb4117d16d8 149 iterator iter2 = iter1 + 1;
madcowswe 0:feb4117d16d8 150 *iter2 = rhs;
madcowswe 0:feb4117d16d8 151 return Initializer<value_type, 2>(iter2 + 1);
madcowswe 0:feb4117d16d8 152 }
madcowswe 0:feb4117d16d8 153
madcowswe 0:feb4117d16d8 154
madcowswe 0:feb4117d16d8 155
madcowswe 0:feb4117d16d8 156 } // namespace tvmet
madcowswe 0:feb4117d16d8 157
madcowswe 0:feb4117d16d8 158
madcowswe 0:feb4117d16d8 159 #endif // TVMET_COMMA_INITIALIZER_H
madcowswe 0:feb4117d16d8 160
madcowswe 0:feb4117d16d8 161 // Local Variables:
madcowswe 0:feb4117d16d8 162 // mode:C++
madcowswe 0:feb4117d16d8 163 // tab-width:8
madcowswe 0:feb4117d16d8 164 // End: