src/foam/primitives/VectorN/vector2/diagTensor2I.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 Type
  25     diagTensor2
  26
  27 Description
  28     DiagTensorN of 2 scalars.
  29
  30 \*---------------------------------------------------------------------------*/
  31
  32
  33 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  34
  35 namespace Foam
  36 {
  37
  38 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  39
  40 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
  41
  42 //- Assign to a SphericalTensorN
  43 template<>
  44 inline void diagTensor2::operator=(const sphericalTensor2& st)
  45 {
  46     this->v_[0] = st[0];
  47     this->v_[1] = st[0];
  48 }
  49
  50 // * * * * * * * * * * * * * * * Global Operators  * * * * * * * * * * * * * //
  51
  52 //- Inner-product between two diagonal tensors
  53 inline diagTensor2
  54 operator&(const diagTensor2& dt1, const diagTensor2& dt2)
  55 {
  56     diagTensor2 result;
  57
  58     result[0] = dt1[0]*dt2[0];
  59     result[1] = dt1[1]*dt2[1];
  60
  61     return result;
  62 }
  63
  64
  65 //- Inner-product between a spherical tensor and a diagonal tensor
  66 inline diagTensor2
  67 operator&(const sphericalTensor2& st1, const diagTensor2& dt2)
  68 {
  69     diagTensor2 result;
  70
  71     result[0] = st1[0]*dt2[0];
  72     result[1] = st1[0]*dt2[1];
  73
  74     return result;
  75 }
  76
  77 //- Inner-product between a diagonal tensor and a spherical tensor
  78 inline diagTensor2
  79 operator&(const diagTensor2& dt1, const sphericalTensor2& st2)
  80 {
  81     diagTensor2 result;
  82
  83     result[0] = dt1[0]*st2[0];
  84     result[1] = dt1[1]*st2[0];
  85
  86     return result;
  87 }
  88
  89
  90 //- Inner-product between a diagonal tensor and a vector
  91 inline vector2
  92 operator&(const diagTensor2& dt, const vector2& v)
  93 {
  94     vector2 result;
  95
  96     result[0] = dt[0]*v[0];
  97     result[1] = dt[1]*v[1];
  98
  99     return result;
 100 }
 101
 102
 103 //- Inner-product between a vector and a diagonal tensor
 104 inline vector2
 105 operator&(const vector2& v, const diagTensor2& dt)
 106 {
 107     vector2 result;
 108
 109     result[0] = v[0]*dt[0];
 110     result[1] = v[1]*dt[1];
 111
 112     return result;
 113 }
 114
 115
 116 //- Division of a scalar by a diagTensor2
 117 inline diagTensor2 operator/(const scalar s, const diagTensor2& dt)
 118 {
 119     diagTensor2 result;
 120
 121     result[0] = s/dt[0];
 122     result[1] = s/dt[1];
 123
 124     return result;
 125 }
 126
 127
 128 //- Return the inverse of a diagTensor2
 129 inline diagTensor2 inv(const diagTensor2& dt)
 130 {
 131     diagTensor2 result;
 132
 133     result[0] = 1.0/dt[0];
 134     result[1] = 1.0/dt[1];
 135
 136     return result;
 137 }
 138
 139
 140 //- Return the component sum
 141 inline scalar sum(const diagTensor2& dt)
 142 {
 143     return dt[0] + dt[1];
 144 }
 145
 146 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 147
 148 } // End namespace Foam
 149
 150 // ************************************************************************* //