src/VectorN/OpenFOAM/Fields/TensorNFields.C

   1 /*---------------------------------------------------------------------------*\
   2   =========                 |
   3   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   4    \\    /   O peration     |
   5     \\  /    A nd           | Copyright held by original author
   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 \*---------------------------------------------------------------------------*/
  26
  27 #include "TensorNFields.H"
  28 #include "transformField.H"
  29
  30 #define TEMPLATE
  31 #include "FieldFunctionsM.C"
  32 #include "ExpandTensorN.H"
  33
  34 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  35
  36 #define TensorN_FieldFunctions(tensorType,diagTensorType,                    \
  37                                sphericalTensorType,vectorType,CmptType,      \
  38                                args...)                                      \
  39                                                                              \
  40 UNARY_FUNCTION(tensorType, tensorType, inv)                                  \
  41 UNARY_FUNCTION(diagTensorType, tensorType, diag)                             \
  42 UNARY_FUNCTION(tensorType, tensorType, negSumDiag)                           \
  43                                                                              \
  44 BINARY_OPERATOR(tensorType, CmptType, tensorType, /, divide)                 \
  45 BINARY_TYPE_OPERATOR(tensorType, CmptType, tensorType, /, divide)            \
  46                                                                              \
  47 BINARY_OPERATOR(vectorType, vectorType, tensorType, /, divide)               \
  48 BINARY_TYPE_OPERATOR(vectorType, vectorType, tensorType, /, divide)          \
  49                                                                              \
  50 BINARY_OPERATOR(tensorType, tensorType, tensorType, /, divide)               \
  51 BINARY_TYPE_OPERATOR(tensorType, tensorType, tensorType, /, divide)          \
  52                                                                              \
  53 BINARY_OPERATOR(tensorType, tensorType, diagTensorType, /, divide)           \
  54 BINARY_TYPE_OPERATOR(tensorType, tensorType, diagTensorType, /, divide)      \
  55                                                                              \
  56 BINARY_OPERATOR(tensorType, diagTensorType, tensorType, /, divide)           \
  57 BINARY_TYPE_OPERATOR(tensorType, diagTensorType, tensorType, /, divide)      \
  58                                                                              \
  59 BINARY_OPERATOR(tensorType, sphericalTensorType, tensorType, /, divide)      \
  60 BINARY_TYPE_OPERATOR(tensorType, sphericalTensorType, tensorType, /, divide) \
  61                                                                              \
  62 BINARY_OPERATOR(tensorType, tensorType, sphericalTensorType, /, divide)      \
  63 BINARY_TYPE_OPERATOR(tensorType, tensorType, sphericalTensorType, /, divide) \
  64                                                                              \
  65 BINARY_OPERATOR(tensorType, tensorType, tensorType, +, add)                  \
  66 BINARY_OPERATOR(tensorType, tensorType, tensorType, -, subtract)             \
  67                                                                              \
  68 BINARY_TYPE_OPERATOR(tensorType, tensorType, tensorType, +, add)             \
  69 BINARY_TYPE_OPERATOR(tensorType, tensorType, tensorType, -, subtract)        \
  70                                                                              \
  71 BINARY_OPERATOR(tensorType, diagTensorType, tensorType, +, add)              \
  72 BINARY_OPERATOR(tensorType, diagTensorType, tensorType, -, subtract)         \
  73                                                                              \
  74 BINARY_TYPE_OPERATOR(tensorType, diagTensorType, tensorType, +, add)         \
  75 BINARY_TYPE_OPERATOR(tensorType, diagTensorType, tensorType, -, subtract)    \
  76                                                                              \
  77 BINARY_OPERATOR(tensorType, tensorType, diagTensorType, +, add)              \
  78 BINARY_OPERATOR(tensorType, tensorType, diagTensorType, -, subtract)         \
  79                                                                              \
  80 BINARY_TYPE_OPERATOR(tensorType, tensorType, diagTensorType, +, add)         \
  81 BINARY_TYPE_OPERATOR(tensorType, tensorType, diagTensorType, -, subtract)    \
  82                                                                              \
  83 BINARY_OPERATOR(tensorType, sphericalTensorType, tensorType, +, add)         \
  84 BINARY_OPERATOR(tensorType, sphericalTensorType, tensorType, -, subtract)    \
  85                                                                              \
  86 BINARY_TYPE_OPERATOR(tensorType, sphericalTensorType, tensorType, +, add)    \
  87 BINARY_TYPE_OPERATOR(tensorType, sphericalTensorType, tensorType, -, subtract) \
  88                                                                              \
  89 BINARY_OPERATOR(tensorType, tensorType, sphericalTensorType, +, add)         \
  90 BINARY_OPERATOR(tensorType, tensorType, sphericalTensorType, -, subtract)    \
  91                                                                              \
  92 BINARY_TYPE_OPERATOR(tensorType, tensorType, sphericalTensorType, +, add)    \
  93 BINARY_TYPE_OPERATOR(tensorType, tensorType, sphericalTensorType, -, subtract) \
  94                                                                              \
  95 template<>                                                                   \
  96 tmp<Field<tensorType> > transformFieldMask<tensorType>                       \
  97 (                                                                            \
  98     const Field<diagTensorType>& dtf                                         \
  99 )                                                                            \
 100 {                                                                            \
 101     tmp<Field<tensorType> > tRes(new Field<tensorType>(dtf.size()));         \
 102     Field<tensorType>& res = tRes();                                         \
 103     TFOR_ALL_F_OP_F(tensorType, res, =, diagTensorType, dtf)                 \
 104     return tRes;                                                             \
 105 }                                                                            \
 106                                                                              \
 107 template<>                                                                   \
 108 tmp<Field<tensorType> > transformFieldMask<tensorType>                       \
 109 (                                                                            \
 110     const Field<sphericalTensorType>& stf                                    \
 111 )                                                                            \
 112 {                                                                            \
 113     tmp<Field<tensorType> > tRes(new Field<tensorType>(stf.size()));         \
 114     Field<tensorType>& res = tRes();                                         \
 115     TFOR_ALL_F_OP_F(tensorType, res, =, sphericalTensorType, stf)            \
 116     return tRes;                                                             \
 117 }                                                                            \
 118
 119
 120 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 121
 122 namespace Foam
 123 {
 124
 125 forAllVectorTensorNTypes(TensorN_FieldFunctions)
 126
 127 }
 128
 129 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 130
 131 #undef TensorN_FieldFunctions
 132
 133 #include "undefFieldFunctionsM.H"
 134
 135 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 136
 137 // ************************************************************************* //