src/OpenFOAM/primitives/VectorSpace/VectorSpaceM.H

   1
   2 template<int N, int I>
   3 class VectorSpaceOps
   4 {
   5 public:
   6
   7     static const int endLoop = (I < N-1) ? 1 : 0;
   8
   9     template<class V, class S, class EqOp>
  10     static inline void eqOpS(V& vs, const S& s, EqOp eo)
  11     {
  12         eo(vs.v_[I], s);
  13         VectorSpaceOps<endLoop*N, endLoop*(I+1)>::eqOpS(vs, s, eo);
  14     }
  15
  16     template<class S, class V, class EqOp>
  17     static inline void SeqOp(S& s, const V& vs, EqOp eo)
  18     {
  19         eo(s, vs.v_[I]);
  20         VectorSpaceOps<endLoop*N, endLoop*(I+1)>::SeqOp(s, vs, eo);
  21     }
  22
  23     template<class V1, class V2, class EqOp>
  24     static inline void eqOp(V1& vs1, const V2& vs2, EqOp eo)
  25     {
  26         eo(vs1.v_[I], vs2.v_[I]);
  27         VectorSpaceOps<endLoop*N, endLoop*(I+1)>::eqOp(vs1, vs2, eo);
  28     }
  29
  30
  31     template<class V, class V1, class S, class Op>
  32     static inline void opVS(V& vs, const V1& vs1, const S& s, Op o)
  33     {
  34         vs.v_[I] = o(vs1.v_[I], s);
  35         VectorSpaceOps<endLoop*N, endLoop*(I+1)>::opVS(vs, vs1, s, o);
  36     }
  37
  38     template<class V, class S, class V1, class Op>
  39     static inline void opSV(V& vs, const S& s,  const V1& vs1, Op o)
  40     {
  41         vs.v_[I] = o(s, vs1.v_[I]);
  42         VectorSpaceOps<endLoop*N, endLoop*(I+1)>::opSV(vs, s, vs1, o);
  43     }
  44
  45     template<class V, class V1, class Op>
  46     static inline void op(V& vs, const V1& vs1, const V1& vs2, Op o)
  47     {
  48         vs.v_[I] = o(vs1.v_[I], vs2.v_[I]);
  49         VectorSpaceOps<endLoop*N, endLoop*(I+1)>::op(vs, vs1, vs2, o);
  50     }
  51 };
  52
  53
  54 template<>
  55 class VectorSpaceOps<0, 0>
  56 {
  57 public:
  58
  59     template<class V, class S, class EqOp>
  60     static inline void eqOpS(V&, const S&, EqOp)
  61     {}
  62
  63     template<class S, class V, class EqOp>
  64     static inline void SeqOp(S&, const V&, EqOp)
  65     {}
  66
  67     template<class V1, class V2, class EqOp>
  68     static inline void eqOp(V1&, const V2&, EqOp)
  69     {}
  70
  71
  72     template<class V, class V1, class S, class Op>
  73     static inline void opVS(V& vs, const V1&, const S&, Op)
  74     {}
  75
  76     template<class V, class S, class V1, class Op>
  77     static inline void opSV(V& vs, const S&, const V1&, Op)
  78     {}
  79
  80     template<class V, class V1, class Op>
  81     static inline void op(V& vs, const V1&, const V1&, Op)
  82     {}
  83 };