src/math/mathvector_structs.h

   1 // ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
   2 // +                                                                      +
   3 // + This file is part of enGrid.                                         +
   4 // +                                                                      +
   5 // + Copyright 2008-2014 enGits GmbH                                      +
   6 // +                                                                      +
   7 // + enGrid is free software: you can redistribute it and/or modify       +
   8 // + it under the terms of the GNU General Public License as published by +
   9 // + the Free Software Foundation, either version 3 of the License, or    +
  10 // + (at your option) any later version.                                  +
  11 // +                                                                      +
  12 // + enGrid is distributed in the hope that it will be useful,            +
  13 // + but WITHOUT ANY WARRANTY; without even the implied warranty of       +
  14 // + MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the        +
  15 // + GNU General Public License for more details.                         +
  16 // +                                                                      +
  17 // + You should have received a copy of the GNU General Public License    +
  18 // + along with enGrid. If not, see <http://www.gnu.org/licenses/>.       +
  19 // +                                                                      +
  20 // ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
  21 template <class T>
  22 struct mv_p
  23 {
  24   static T apply(const T &a, const T &b) { return a+b; };
  25 };
  26
  27 template <class T>
  28 struct mv_m
  29 {
  30   static T apply(const T &a, const T &b) { return a-b; };
  31 };
  32
  33 template <class T>
  34 struct mv_ml
  35 {
  36   static T apply(const T &a, const T &b) { return a*b; };
  37 };
  38
  39 template <class L, class O, class R>
  40 struct ParseNode
  41 {
  42   typedef typename R::value_type value_type;
  43   const L &l;
  44   const R &r;
  45   ParseNode(const L &a, const R &b) : l(a), r(b) {};
  46   value_type operator[](const uint_t &i) const { return O::apply(l[i], r[i]); };
  47   uint_t size() const { return r.size(); };
  48   value_type abs() const;
  49   value_type abs2() const;
  50 };
  51
  52 template <class O, class R>
  53 struct ParseNode<double, O, R>
  54 {
  55   typedef typename R::value_type value_type;
  56   const double l;
  57   const R &r;
  58   ParseNode(const double a, const R &b) : l(a), r(b) {};
  59   value_type operator[](const uint_t &i) const { return O::apply(l, r[i]); };
  60   uint_t size() const { return r.size(); };
  61   value_type abs() const;
  62   value_type abs2() const;
  63 };