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

VectorUnaryFunctions.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: VectorUnaryFunctions.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_VECTOR_UNARY_FUNCTIONS_H





madcowswe

0:feb4117d16d8


    25

#define TVMET_VECTOR_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(Vector<T, Sz>)





madcowswe

0:feb4117d16d8


    36

 */





madcowswe

0:feb4117d16d8


    37

#define TVMET_DECLARE_MACRO(NAME)            \





madcowswe

0:feb4117d16d8


    38

template<class T, std::size_t Sz>            \





madcowswe

0:feb4117d16d8


    39

XprVector<                        \





madcowswe

0:feb4117d16d8


    40

  XprUnOp<                        \





madcowswe

0:feb4117d16d8


    41

    Fcnl_##NAME<T>,                    \





madcowswe

0:feb4117d16d8


    42

    VectorConstReference<T, Sz>                \





madcowswe

0:feb4117d16d8


    43

  >,                            \





madcowswe

0:feb4117d16d8


    44

  Sz                            \





madcowswe

0:feb4117d16d8


    45

>                            \





madcowswe

0:feb4117d16d8


    46

NAME(const Vector<T, Sz>& rhs) TVMET_CXX_ALWAYS_INLINE;





madcowswe

0:feb4117d16d8


    47







madcowswe

0:feb4117d16d8


    48

TVMET_DECLARE_MACRO(abs)





madcowswe

0:feb4117d16d8


    49

TVMET_DECLARE_MACRO(cbrt)





madcowswe

0:feb4117d16d8


    50

TVMET_DECLARE_MACRO(ceil)





madcowswe

0:feb4117d16d8


    51

TVMET_DECLARE_MACRO(floor)





madcowswe

0:feb4117d16d8


    52

TVMET_DECLARE_MACRO(rint)





madcowswe

0:feb4117d16d8


    53

TVMET_DECLARE_MACRO(sin)





madcowswe

0:feb4117d16d8


    54

TVMET_DECLARE_MACRO(cos)





madcowswe

0:feb4117d16d8


    55

TVMET_DECLARE_MACRO(tan)





madcowswe

0:feb4117d16d8


    56

TVMET_DECLARE_MACRO(sinh)





madcowswe

0:feb4117d16d8


    57

TVMET_DECLARE_MACRO(cosh)





madcowswe

0:feb4117d16d8


    58

TVMET_DECLARE_MACRO(tanh)





madcowswe

0:feb4117d16d8


    59

TVMET_DECLARE_MACRO(asin)





madcowswe

0:feb4117d16d8


    60

TVMET_DECLARE_MACRO(acos)





madcowswe

0:feb4117d16d8


    61

TVMET_DECLARE_MACRO(atan)





madcowswe

0:feb4117d16d8


    62

TVMET_DECLARE_MACRO(exp)





madcowswe

0:feb4117d16d8


    63

TVMET_DECLARE_MACRO(log)





madcowswe

0:feb4117d16d8


    64

TVMET_DECLARE_MACRO(log10)





madcowswe

0:feb4117d16d8


    65

TVMET_DECLARE_MACRO(sqrt)





madcowswe

0:feb4117d16d8


    66







madcowswe

0:feb4117d16d8


    67

#if defined(TVMET_HAVE_IEEE_MATH)





madcowswe

0:feb4117d16d8


    68

TVMET_DECLARE_MACRO(asinh)





madcowswe

0:feb4117d16d8


    69

TVMET_DECLARE_MACRO(acosh)





madcowswe

0:feb4117d16d8


    70

TVMET_DECLARE_MACRO(atanh)





madcowswe

0:feb4117d16d8


    71

TVMET_DECLARE_MACRO(expm1)





madcowswe

0:feb4117d16d8


    72

TVMET_DECLARE_MACRO(log1p)





madcowswe

0:feb4117d16d8


    73

TVMET_DECLARE_MACRO(erf)





madcowswe

0:feb4117d16d8


    74

TVMET_DECLARE_MACRO(erfc)





madcowswe

0:feb4117d16d8


    75

TVMET_DECLARE_MACRO(j0)





madcowswe

0:feb4117d16d8


    76

TVMET_DECLARE_MACRO(j1)





madcowswe

0:feb4117d16d8


    77

TVMET_DECLARE_MACRO(y0)





madcowswe

0:feb4117d16d8


    78

TVMET_DECLARE_MACRO(y1)





madcowswe

0:feb4117d16d8


    79

TVMET_DECLARE_MACRO(lgamma)





madcowswe

0:feb4117d16d8


    80

/** \todo isnan etc. - default return is only an int! */





madcowswe

0:feb4117d16d8


    81







madcowswe

0:feb4117d16d8


    82

TVMET_DECLARE_MACRO(finite)





madcowswe

0:feb4117d16d8


    83

#endif // defined(TVMET_HAVE_IEEE_MATH)





madcowswe

0:feb4117d16d8


    84







madcowswe

0:feb4117d16d8


    85

#undef TVMET_DECLARE_MACRO





madcowswe

0:feb4117d16d8


    86







madcowswe

0:feb4117d16d8


    87







madcowswe

0:feb4117d16d8


    88

/*





madcowswe

0:feb4117d16d8


    89

 * unary_function(Vector<std::complex<T>, Sz>)





madcowswe

0:feb4117d16d8


    90

 */





madcowswe

0:feb4117d16d8


    91

#if defined(TVMET_HAVE_COMPLEX)





madcowswe

0:feb4117d16d8


    92

#define TVMET_DECLARE_MACRO(NAME)                    \





madcowswe

0:feb4117d16d8


    93

template<class T, std::size_t Sz>                    \





madcowswe

0:feb4117d16d8


    94

XprVector<                                \





madcowswe

0:feb4117d16d8


    95

  XprUnOp<                                \





madcowswe

0:feb4117d16d8


    96

    Fcnl_##NAME< std::complex<T> >,                    \





madcowswe

0:feb4117d16d8


    97

    VectorConstReference<std::complex<T>, Sz>                \





madcowswe

0:feb4117d16d8


    98

  >,                                    \





madcowswe

0:feb4117d16d8


    99

  Sz                                    \





madcowswe

0:feb4117d16d8


   100

>                                    \





madcowswe

0:feb4117d16d8


   101

NAME(const Vector<std::complex<T>, Sz>& rhs) TVMET_CXX_ALWAYS_INLINE;





madcowswe

0:feb4117d16d8


   102







madcowswe

0:feb4117d16d8


   103

TVMET_DECLARE_MACRO(real)





madcowswe

0:feb4117d16d8


   104

TVMET_DECLARE_MACRO(imag)





madcowswe

0:feb4117d16d8


   105

TVMET_DECLARE_MACRO(arg)





madcowswe

0:feb4117d16d8


   106

TVMET_DECLARE_MACRO(norm)





madcowswe

0:feb4117d16d8


   107

TVMET_DECLARE_MACRO(conj)





madcowswe

0:feb4117d16d8


   108







madcowswe

0:feb4117d16d8


   109

#undef TVMET_DECLARE_MACRO





madcowswe

0:feb4117d16d8


   110







madcowswe

0:feb4117d16d8


   111

#endif // defined(TVMET_HAVE_COMPLEX)





madcowswe

0:feb4117d16d8


   112







madcowswe

0:feb4117d16d8


   113







madcowswe

0:feb4117d16d8


   114

/*********************************************************





madcowswe

0:feb4117d16d8


   115

 * PART II: IMPLEMENTATION





madcowswe

0:feb4117d16d8


   116

 *********************************************************/





madcowswe

0:feb4117d16d8


   117







madcowswe

0:feb4117d16d8


   118







madcowswe

0:feb4117d16d8


   119

/*





madcowswe

0:feb4117d16d8


   120

 * unary_function(Vector<T, Sz>)





madcowswe

0:feb4117d16d8


   121

 */





madcowswe

0:feb4117d16d8


   122

#define TVMET_IMPLEMENT_MACRO(NAME)                    \





madcowswe

0:feb4117d16d8


   123

template<class T, std::size_t Sz>                    \





madcowswe

0:feb4117d16d8


   124

inline                                    \





madcowswe

0:feb4117d16d8


   125

XprVector<                                \





madcowswe

0:feb4117d16d8


   126

  XprUnOp<                                \





madcowswe

0:feb4117d16d8


   127

    Fcnl_##NAME<T>,                            \





madcowswe

0:feb4117d16d8


   128

    VectorConstReference<T, Sz>                        \





madcowswe

0:feb4117d16d8


   129

  >,                                    \





madcowswe

0:feb4117d16d8


   130

  Sz                                    \





madcowswe

0:feb4117d16d8


   131

>                                    \





madcowswe

0:feb4117d16d8


   132

NAME(const Vector<T, Sz>& rhs) {                    \





madcowswe

0:feb4117d16d8


   133

  typedef XprUnOp<                            \





madcowswe

0:feb4117d16d8


   134

      Fcnl_##NAME<T>,                            \





madcowswe

0:feb4117d16d8


   135

      VectorConstReference<T, Sz>                    \





madcowswe

0:feb4117d16d8


   136

    >                             expr_type;    \





madcowswe

0:feb4117d16d8


   137

    return XprVector<expr_type, Sz>(expr_type(rhs.const_ref()));    \





madcowswe

0:feb4117d16d8


   138

}





madcowswe

0:feb4117d16d8


   139







madcowswe

0:feb4117d16d8


   140

TVMET_IMPLEMENT_MACRO(abs)





madcowswe

0:feb4117d16d8


   141

TVMET_IMPLEMENT_MACRO(cbrt)





madcowswe

0:feb4117d16d8


   142

TVMET_IMPLEMENT_MACRO(ceil)





madcowswe

0:feb4117d16d8


   143

TVMET_IMPLEMENT_MACRO(floor)





madcowswe

0:feb4117d16d8


   144

TVMET_IMPLEMENT_MACRO(rint)





madcowswe

0:feb4117d16d8


   145

TVMET_IMPLEMENT_MACRO(sin)





madcowswe

0:feb4117d16d8


   146

TVMET_IMPLEMENT_MACRO(cos)





madcowswe

0:feb4117d16d8


   147

TVMET_IMPLEMENT_MACRO(tan)





madcowswe

0:feb4117d16d8


   148

TVMET_IMPLEMENT_MACRO(sinh)





madcowswe

0:feb4117d16d8


   149

TVMET_IMPLEMENT_MACRO(cosh)





madcowswe

0:feb4117d16d8


   150

TVMET_IMPLEMENT_MACRO(tanh)





madcowswe

0:feb4117d16d8


   151

TVMET_IMPLEMENT_MACRO(asin)





madcowswe

0:feb4117d16d8


   152

TVMET_IMPLEMENT_MACRO(acos)





madcowswe

0:feb4117d16d8


   153

TVMET_IMPLEMENT_MACRO(atan)





madcowswe

0:feb4117d16d8


   154

TVMET_IMPLEMENT_MACRO(exp)





madcowswe

0:feb4117d16d8


   155

TVMET_IMPLEMENT_MACRO(log)





madcowswe

0:feb4117d16d8


   156

TVMET_IMPLEMENT_MACRO(log10)





madcowswe

0:feb4117d16d8


   157

TVMET_IMPLEMENT_MACRO(sqrt)





madcowswe

0:feb4117d16d8


   158







madcowswe

0:feb4117d16d8


   159

#if defined(TVMET_HAVE_IEEE_MATH)





madcowswe

0:feb4117d16d8


   160

TVMET_IMPLEMENT_MACRO(asinh)





madcowswe

0:feb4117d16d8


   161

TVMET_IMPLEMENT_MACRO(acosh)





madcowswe

0:feb4117d16d8


   162

TVMET_IMPLEMENT_MACRO(atanh)





madcowswe

0:feb4117d16d8


   163

TVMET_IMPLEMENT_MACRO(expm1)





madcowswe

0:feb4117d16d8


   164

TVMET_IMPLEMENT_MACRO(log1p)





madcowswe

0:feb4117d16d8


   165

TVMET_IMPLEMENT_MACRO(erf)





madcowswe

0:feb4117d16d8


   166

TVMET_IMPLEMENT_MACRO(erfc)





madcowswe

0:feb4117d16d8


   167

TVMET_IMPLEMENT_MACRO(j0)





madcowswe

0:feb4117d16d8


   168

TVMET_IMPLEMENT_MACRO(j1)





madcowswe

0:feb4117d16d8


   169

TVMET_IMPLEMENT_MACRO(y0)





madcowswe

0:feb4117d16d8


   170

TVMET_IMPLEMENT_MACRO(y1)





madcowswe

0:feb4117d16d8


   171

TVMET_IMPLEMENT_MACRO(lgamma)





madcowswe

0:feb4117d16d8


   172

/** \todo isnan etc. - default return is only an int! */





madcowswe

0:feb4117d16d8


   173







madcowswe

0:feb4117d16d8


   174

TVMET_IMPLEMENT_MACRO(finite)





madcowswe

0:feb4117d16d8


   175

#endif // defined(TVMET_HAVE_IEEE_MATH)





madcowswe

0:feb4117d16d8


   176







madcowswe

0:feb4117d16d8


   177

#undef TVMET_IMPLEMENT_MACRO





madcowswe

0:feb4117d16d8


   178







madcowswe

0:feb4117d16d8


   179







madcowswe

0:feb4117d16d8


   180

/*





madcowswe

0:feb4117d16d8


   181

 * unary_function(Vector<std::complex<T>, Sz>)





madcowswe

0:feb4117d16d8


   182

 */





madcowswe

0:feb4117d16d8


   183

#if defined(TVMET_HAVE_COMPLEX)





madcowswe

0:feb4117d16d8


   184

#define TVMET_IMPLEMENT_MACRO(NAME)                    \





madcowswe

0:feb4117d16d8


   185

template<class T, std::size_t Sz>                    \





madcowswe

0:feb4117d16d8


   186

inline                                    \





madcowswe

0:feb4117d16d8


   187

XprVector<                                \





madcowswe

0:feb4117d16d8


   188

  XprUnOp<                                \





madcowswe

0:feb4117d16d8


   189

    Fcnl_##NAME< std::complex<T> >,                    \





madcowswe

0:feb4117d16d8


   190

    VectorConstReference<std::complex<T>, Sz>                \





madcowswe

0:feb4117d16d8


   191

  >,                                    \





madcowswe

0:feb4117d16d8


   192

  Sz                                    \





madcowswe

0:feb4117d16d8


   193

>                                    \





madcowswe

0:feb4117d16d8


   194

NAME(const Vector<std::complex<T>, Sz>& rhs) {                \





madcowswe

0:feb4117d16d8


   195

  typedef XprUnOp<                            \





madcowswe

0:feb4117d16d8


   196

      Fcnl_##NAME< std::complex<T> >,                    \





madcowswe

0:feb4117d16d8


   197

      VectorConstReference<std::complex<T>, Sz>                \





madcowswe

0:feb4117d16d8


   198

    >                             expr_type;    \





madcowswe

0:feb4117d16d8


   199

    return XprVector<expr_type, Sz>(expr_type(rhs.const_ref()));    \





madcowswe

0:feb4117d16d8


   200

}





madcowswe

0:feb4117d16d8


   201







madcowswe

0:feb4117d16d8


   202

TVMET_IMPLEMENT_MACRO(real)





madcowswe

0:feb4117d16d8


   203

TVMET_IMPLEMENT_MACRO(imag)





madcowswe

0:feb4117d16d8


   204

TVMET_IMPLEMENT_MACRO(arg)





madcowswe

0:feb4117d16d8


   205

TVMET_IMPLEMENT_MACRO(norm)





madcowswe

0:feb4117d16d8


   206

TVMET_IMPLEMENT_MACRO(conj)





madcowswe

0:feb4117d16d8


   207







madcowswe

0:feb4117d16d8


   208

#undef TVMET_IMPLEMENT_MACRO





madcowswe

0:feb4117d16d8


   209







madcowswe

0:feb4117d16d8


   210

#endif // defined(TVMET_HAVE_COMPLEX)





madcowswe

0:feb4117d16d8


   211







madcowswe

0:feb4117d16d8


   212







madcowswe

0:feb4117d16d8


   213

} // namespace tvmet





madcowswe

0:feb4117d16d8


   214







madcowswe

0:feb4117d16d8


   215

#endif // TVMET_VECTOR_UNARY_FUNCTIONS_H





madcowswe

0:feb4117d16d8


   216







madcowswe

0:feb4117d16d8


   217

// Local Variables:





madcowswe

0:feb4117d16d8


   218

// mode:C++





madcowswe

0:feb4117d16d8


   219

// tab-width:8





madcowswe

0:feb4117d16d8


   220

// End: