src/foam/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     template<class V, class V1, class V2, class Op>
  53     static inline void opVV(V& vs, const V1& vs1, const V2& vs2, Op o)
  54     {
  55         vs.v_[I] = o(vs1.v_[I], vs2.v_[I]);
  56         VectorSpaceOps<endLoop*N, endLoop*(I+1)>::opVV(vs, vs1, vs2, o);
  57     }
  58
  59     template<class S, class V1, class V2, class EqOp, class Op>
  60     static inline void SopEqOpVV(S& s, const V1& vs1, const V2& vs2, EqOp eo, Op o)
  61     {
  62         eo(s, o(vs1.v_[I], vs2.v_[I]));
  63         VectorSpaceOps<endLoop*N, endLoop*(I+1)>::SopEqOpVV(s, vs1, vs2, eo, o);
  64     }
  65 };
  66
  67
  68 template<>
  69 class VectorSpaceOps<0, 0>
  70 {
  71 public:
  72
  73     template<class V, class S, class EqOp>
  74     static inline void eqOpS(V&, const S&, EqOp)
  75     {}
  76
  77     template<class S, class V, class EqOp>
  78     static inline void SeqOp(S&, const V&, EqOp)
  79     {}
  80
  81     template<class V1, class V2, class EqOp>
  82     static inline void eqOp(V1&, const V2&, EqOp)
  83     {}
  84
  85
  86     template<class V, class V1, class S, class Op>
  87     static inline void opVS(V& vs, const V1&, const S&, Op)
  88     {}
  89
  90     template<class V, class S, class V1, class Op>
  91     static inline void opSV(V& vs, const S&, const V1&, Op)
  92     {}
  93
  94     template<class V, class V1, class Op>
  95     static inline void op(V& vs, const V1&, const V1&, Op)
  96     {}
  97
  98     template<class V, class V1, class V2, class Op>
  99     static inline void opVV(V& vs, const V1& vs1, const V2& vs2, Op o)
 100     {}
 101
 102     template<class S, class V1, class V2, class EqOp, class Op>
 103     static inline void SopEqOpVV(S& s, const V1& vs1, const V2& vs2, EqOp eo, Op o)
 104     {}
 105 };
 106