src/foam/primitives/VectorSpace/products.H

   1 /*---------------------------------------------------------------------------*\
   2   =========                 |
   3   \\      /  F ield         | foam-extend: Open Source CFD
   4    \\    /   O peration     | Version:     3.2
   5     \\  /    A nd           | Web:         http://www.foam-extend.org
   6      \\/     M anipulation  | For copyright notice see file Copyright
   7 -------------------------------------------------------------------------------
   8 License
   9     This file is part of foam-extend.
  10
  11     foam-extend 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 3 of the License, or (at your
  14     option) any later version.
  15
  16     foam-extend is distributed in the hope that it will be useful, but
  17     WITHOUT ANY WARRANTY; without even the implied warranty of
  18     MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  19     General Public License for more details.
  20
  21     You should have received a copy of the GNU General Public License
  22     along with foam-extend.  If not, see <http://www.gnu.org/licenses/>.
  23
  24 InNamespace
  25     Foam
  26
  27 Description
  28     Traits classes for inner and outer products of primitives.
  29
  30 \*---------------------------------------------------------------------------*/
  31
  32 #ifndef products_H
  33 #define products_H
  34
  35 #include "pTraits.H"
  36
  37 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  38
  39 namespace Foam
  40 {
  41
  42 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  43
  44 template<class Cmpt, int rank>
  45 class typeOfRank
  46 {};
  47
  48
  49 template<class Cmpt, int rank>
  50 class symmTypeOfRank
  51 {};
  52
  53
  54 template<class arg1, class arg2>
  55 class typeOfSum
  56 {
  57 public:
  58
  59     typedef arg1 type;
  60 };
  61
  62
  63 template<class arg1, class arg2>
  64 class outerProduct
  65 {
  66 public:
  67
  68     typedef typename typeOfRank
  69     <
  70         typename pTraits<arg1>::cmptType,
  71         int(pTraits<arg1>::rank) + int(pTraits<arg2>::rank)
  72     >::type type;
  73 };
  74
  75
  76 template<class arg1, class arg2>
  77 class crossProduct
  78 {
  79 public:
  80
  81     typedef typename typeOfRank
  82     <
  83         typename pTraits<arg2>::cmptType,
  84         int(pTraits<arg1>::rank) + int(pTraits<arg2>::rank) - 1
  85     >::type type;
  86 };
  87
  88 template<class arg1, class arg2>
  89 class innerProduct
  90 {
  91 public:
  92
  93     typedef typename typeOfRank
  94     <
  95         typename pTraits<arg1>::cmptType,
  96         int(pTraits<arg1>::rank) + int(pTraits<arg2>::rank) - 2
  97     >::type type;
  98 };
  99
 100 template<class arg1, class arg2>
 101 class scalarProduct
 102 {
 103 public:
 104
 105     typedef typename pTraits<arg1>::cmptType type;
 106 };
 107
 108
 109 template<class arg1, int arg2>
 110 class powProduct
 111 {
 112 public:
 113
 114     typedef typename symmTypeOfRank
 115     <
 116         typename pTraits<arg1>::cmptType,
 117         arg2*int(pTraits<arg1>::rank)
 118     >::type type;
 119 };
 120
 121
 122 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 123
 124 } // End namespace Foam
 125
 126 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 127
 128 #endif
 129
 130 // ************************************************************************* //