Oskar Weigl / tvmet

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

AliasProxy.h@0:feb4117d16d8, 2012-03-28 (annotated)

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: AliasProxy.h,v 1.8 2007-06-23 15:58:58 opetzold Exp $





madcowswe

0:feb4117d16d8


    22

 */





madcowswe

0:feb4117d16d8


    23







madcowswe

0:feb4117d16d8


    24

#ifndef TVMET_ALIAS_PROXY_H





madcowswe

0:feb4117d16d8


    25

#define TVMET_ALIAS_PROXY_H





madcowswe

0:feb4117d16d8


    26







madcowswe

0:feb4117d16d8


    27

namespace tvmet {





madcowswe

0:feb4117d16d8


    28







madcowswe

0:feb4117d16d8


    29







madcowswe

0:feb4117d16d8


    30

/** forwards */





madcowswe

0:feb4117d16d8


    31

template<class E> class AliasProxy;





madcowswe

0:feb4117d16d8


    32







madcowswe

0:feb4117d16d8


    33







madcowswe

0:feb4117d16d8


    34

/**





madcowswe

0:feb4117d16d8


    35

 * \brief Simplify syntax for alias Matrices and Vectors,





madcowswe

0:feb4117d16d8


    36

 *        where aliasing left hand values appear in the





madcowswe

0:feb4117d16d8


    37

 *        expression.





madcowswe

0:feb4117d16d8


    38

 * \par Example:





madcowswe

0:feb4117d16d8


    39

 * \code





madcowswe

0:feb4117d16d8


    40

 * typedef tvmet::Matrix<double, 10, 10>	matrix_type;





madcowswe

0:feb4117d16d8


    41

 * matrix_type					m;





madcowswe

0:feb4117d16d8


    42

 * ...





madcowswe

0:feb4117d16d8


    43

 * alias(m) += trans(m);





madcowswe

0:feb4117d16d8


    44

 * \endcode





madcowswe

0:feb4117d16d8


    45

 * \sa AliasProxy





madcowswe

0:feb4117d16d8


    46

 * \sa Some Notes \ref alias





madcowswe

0:feb4117d16d8


    47

 */





madcowswe

0:feb4117d16d8


    48

template<class E>





madcowswe

0:feb4117d16d8


    49

AliasProxy<E> alias(E& expr) { return AliasProxy<E>(expr); }





madcowswe

0:feb4117d16d8


    50







madcowswe

0:feb4117d16d8


    51







madcowswe

0:feb4117d16d8


    52

/**





madcowswe

0:feb4117d16d8


    53

 * \class AliasProxy AliasProxy.h "tvmet/AliasProxy.h"





madcowswe

0:feb4117d16d8


    54

 * \brief Assign proxy for alias Matrices and Vectors.





madcowswe

0:feb4117d16d8


    55

 *





madcowswe

0:feb4117d16d8


    56

 *        A short lived object to provide simplified alias syntax.





madcowswe

0:feb4117d16d8


    57

 *        Only the friend function alias is allowed to create





madcowswe

0:feb4117d16d8


    58

 *        such a object. The proxy calls the appropriate member





madcowswe

0:feb4117d16d8


    59

 *        alias_xyz() which have to use temporaries to avoid





madcowswe

0:feb4117d16d8


    60

 *        overlapping memory regions.





madcowswe

0:feb4117d16d8


    61

 * \sa alias





madcowswe

0:feb4117d16d8


    62

 * \sa Some Notes \ref alias





madcowswe

0:feb4117d16d8


    63

 * \note Thanks to ublas-dev group, where the principle idea





madcowswe

0:feb4117d16d8


    64

 *       comes from.





madcowswe

0:feb4117d16d8


    65

 */





madcowswe

0:feb4117d16d8


    66

template<class E>





madcowswe

0:feb4117d16d8


    67

class AliasProxy





madcowswe

0:feb4117d16d8


    68

{





madcowswe

0:feb4117d16d8


    69

  AliasProxy(const AliasProxy&);





madcowswe

0:feb4117d16d8


    70

  AliasProxy& operator=(const AliasProxy&);





madcowswe

0:feb4117d16d8


    71







madcowswe

0:feb4117d16d8


    72

  friend AliasProxy<E> alias<>(E& expr);





madcowswe

0:feb4117d16d8


    73







madcowswe

0:feb4117d16d8


    74

public:





madcowswe

0:feb4117d16d8


    75

  AliasProxy(E& expr) : m_expr(expr) { }





madcowswe

0:feb4117d16d8


    76







madcowswe

0:feb4117d16d8


    77







madcowswe

0:feb4117d16d8


    78

  template<class E2>





madcowswe

0:feb4117d16d8


    79

  E& operator=(const E2& expr) {





madcowswe

0:feb4117d16d8


    80

    return m_expr.alias_assign(expr);





madcowswe

0:feb4117d16d8


    81

  }





madcowswe

0:feb4117d16d8


    82







madcowswe

0:feb4117d16d8


    83

  template<class E2>





madcowswe

0:feb4117d16d8


    84

  E& operator+=(const E2& expr) {





madcowswe

0:feb4117d16d8


    85

    return m_expr.alias_add_eq(expr);





madcowswe

0:feb4117d16d8


    86

  }





madcowswe

0:feb4117d16d8


    87







madcowswe

0:feb4117d16d8


    88

  template<class E2>





madcowswe

0:feb4117d16d8


    89

  E& operator-=(const E2& expr) {





madcowswe

0:feb4117d16d8


    90

    return m_expr.alias_sub_eq(expr);





madcowswe

0:feb4117d16d8


    91

  }





madcowswe

0:feb4117d16d8


    92







madcowswe

0:feb4117d16d8


    93

  template<class E2>





madcowswe

0:feb4117d16d8


    94

  E& operator*=(const E2& expr) {





madcowswe

0:feb4117d16d8


    95

    return m_expr.alias_mul_eq(expr);





madcowswe

0:feb4117d16d8


    96

  }





madcowswe

0:feb4117d16d8


    97







madcowswe

0:feb4117d16d8


    98

  template<class E2>





madcowswe

0:feb4117d16d8


    99

  E& operator/=(const E2& expr) {





madcowswe

0:feb4117d16d8


   100

    return m_expr.alias_div_eq(expr);





madcowswe

0:feb4117d16d8


   101

  }





madcowswe

0:feb4117d16d8


   102







madcowswe

0:feb4117d16d8


   103

private:





madcowswe

0:feb4117d16d8


   104

  E&						m_expr;





madcowswe

0:feb4117d16d8


   105

};





madcowswe

0:feb4117d16d8


   106







madcowswe

0:feb4117d16d8


   107







madcowswe

0:feb4117d16d8


   108

#if 0





madcowswe

0:feb4117d16d8


   109

namespace element_wise {





madcowswe

0:feb4117d16d8


   110

// \todo to write





madcowswe

0:feb4117d16d8


   111

template<class E, class E2>





madcowswe

0:feb4117d16d8


   112

E& operator/=(AliasProxy<E>& proxy, const E2& rhs) {





madcowswe

0:feb4117d16d8


   113

  return proxy.div_upd(rhs);





madcowswe

0:feb4117d16d8


   114

}





madcowswe

0:feb4117d16d8


   115







madcowswe

0:feb4117d16d8


   116

}





madcowswe

0:feb4117d16d8


   117

#endif





madcowswe

0:feb4117d16d8


   118







madcowswe

0:feb4117d16d8


   119







madcowswe

0:feb4117d16d8


   120

} // namespace tvmet





madcowswe

0:feb4117d16d8


   121







madcowswe

0:feb4117d16d8


   122







madcowswe

0:feb4117d16d8


   123

#endif /* TVMET_ALIAS_PROXY_H */





madcowswe

0:feb4117d16d8


   124







madcowswe

0:feb4117d16d8


   125

// Local Variables:





madcowswe

0:feb4117d16d8


   126

// mode:C++





madcowswe

0:feb4117d16d8


   127

// tab-width:8





madcowswe

0:feb4117d16d8


   128

// End: