src/OpenFOAM/primitives/VectorSpace/products.H

   1 /*---------------------------------------------------------------------------*\
   2   =========                 |
   3   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   4    \\    /   O peration     |
   5     \\  /    A nd           | Copyright (C) 1991-2008 OpenCFD Ltd.
   6      \\/     M anipulation  |
   7 -------------------------------------------------------------------------------
   8 License
   9     This file is part of OpenFOAM.
  10
  11     OpenFOAM is free software; you can redistribute it and/or modify it
  12     under the terms of the GNU General Public License as published by the
  13     Free Software Foundation; either version 2 of the License, or (at your
  14     option) any later version.
  15
  16     OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
  17     ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
  18     FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License
  19     for more details.
  20
  21     You should have received a copy of the GNU General Public License
  22     along with OpenFOAM; if not, write to the Free Software Foundation,
  23     Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA
  24
  25 InNamespace
  26     Foam
  27
  28 Description
  29     Traits classes for inner and outer products of primitives.
  30
  31 \*---------------------------------------------------------------------------*/
  32
  33 #ifndef products_H
  34 #define products_H
  35
  36 #include "pTraits.H"
  37
  38 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  39
  40 namespace Foam
  41 {
  42
  43 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  44
  45 template<class Cmpt, int rank>
  46 class typeOfRank
  47 {};
  48
  49
  50 template<class Cmpt, int rank>
  51 class symmTypeOfRank
  52 {};
  53
  54
  55 template<class arg1, class arg2>
  56 class typeOfSum
  57 {
  58 public:
  59
  60     typedef arg1 type;
  61 };
  62
  63
  64 #ifndef __CINT__
  65
  66 template<class arg1, class arg2>
  67 class outerProduct
  68 {
  69 public:
  70
  71     typedef typename typeOfRank
  72     <
  73         typename pTraits<arg1>::cmptType,
  74         int(pTraits<arg1>::rank) + int(pTraits<arg2>::rank)
  75     >::type type;
  76 };
  77
  78
  79 template<class arg1, class arg2>
  80 class crossProduct
  81 {
  82 public:
  83
  84     typedef typename typeOfRank
  85     <
  86         typename pTraits<arg2>::cmptType,
  87         int(pTraits<arg1>::rank) + int(pTraits<arg2>::rank) - 1
  88     >::type type;
  89 };
  90
  91 template<class arg1, class arg2>
  92 class innerProduct
  93 {
  94 public:
  95
  96     typedef typename typeOfRank
  97     <
  98         typename pTraits<arg1>::cmptType,
  99         int(pTraits<arg1>::rank) + int(pTraits<arg2>::rank) - 2
 100     >::type type;
 101 };
 102
 103 template<class arg1, class arg2>
 104 class scalarProduct
 105 {
 106 public:
 107
 108     typedef typename pTraits<arg1>::cmptType type;
 109 };
 110
 111
 112 template<class arg1, int arg2>
 113 class powProduct
 114 {
 115 public:
 116
 117     typedef typename symmTypeOfRank
 118     <
 119         typename pTraits<arg1>::cmptType,
 120         arg2*int(pTraits<arg1>::rank)
 121     >::type type;
 122 };
 123
 124 #else
 125
 126 template<class arg1, class arg2>
 127 class outerProduct
 128 {
 129 public:
 130
 131     typedef arg1 type;
 132 };
 133
 134
 135 template<class arg1, class arg2>
 136 class crossProduct
 137 {
 138 public:
 139
 140     typedef arg1 type;
 141 };
 142
 143 template<class arg1, class arg2>
 144 class innerProduct
 145 {
 146 public:
 147
 148     typedef arg1 type;
 149 };
 150
 151 template<class arg1, class arg2>
 152 class scalarProduct
 153 {
 154 public:
 155
 156     typedef typename pTraits<arg1>::cmptType type;
 157 };
 158
 159
 160 template<class arg1, int arg2>
 161 class powProduct
 162 {
 163 public:
 164
 165     typedef arg1 type;
 166 };
 167
 168 #endif
 169
 170
 171 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 172
 173 } // End namespace Foam
 174
 175 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 176
 177 #endif
 178
 179 // ************************************************************************* //