symmTensorField.C
Go to the documentation of this file.
1 /*---------------------------------------------------------------------------*\
2  ========= |
3  \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
4  \\ / O peration |
5  \\ / A nd | Copyright (C) 2011-2016 OpenFOAM Foundation
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
13  the Free Software Foundation, either version 3 of the License, or
14  (at your 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, see <http://www.gnu.org/licenses/>.
23 
24 \*---------------------------------------------------------------------------*/
25 
26 #include "symmTensorField.H"
27 #include "transformField.H"
28 
29 #define TEMPLATE
30 #include "FieldFunctionsM.C"
31 
32 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
33 
34 namespace Foam
35 {
36 
37 // * * * * * * * * * * * * * * * global functions * * * * * * * * * * * * * //
38 
41 
50 
51 void inv(Field<symmTensor>& tf, const UList<symmTensor>& tf1)
52 {
53  if (tf.empty())
54  {
55  return;
56  }
57 
58  scalar scale = magSqr(tf1[0]);
59  Vector<bool> removeCmpts
60  (
61  magSqr(tf1[0].xx())/scale < SMALL,
62  magSqr(tf1[0].yy())/scale < SMALL,
63  magSqr(tf1[0].zz())/scale < SMALL
64  );
65 
66  if (removeCmpts.x() || removeCmpts.y() || removeCmpts.z())
67  {
68  symmTensorField tf1Plus(tf1);
69 
70  if (removeCmpts.x())
71  {
72  tf1Plus += symmTensor(1,0,0,0,0,0);
73  }
74 
75  if (removeCmpts.y())
76  {
77  tf1Plus += symmTensor(0,0,0,1,0,0);
78  }
79 
80  if (removeCmpts.z())
81  {
82  tf1Plus += symmTensor(0,0,0,0,0,1);
83  }
84 
86 
87  if (removeCmpts.x())
88  {
89  tf -= symmTensor(1,0,0,0,0,0);
90  }
91 
92  if (removeCmpts.y())
93  {
94  tf -= symmTensor(0,0,0,1,0,0);
95  }
96 
97  if (removeCmpts.z())
98  {
99  tf -= symmTensor(0,0,0,0,0,1);
100  }
101  }
102  else
103  {
105  }
106 }
107 
109 {
110  tmp<symmTensorField> result(new symmTensorField(tf.size()));
111  inv(result.ref(), tf);
112  return result;
113 }
114 
116 {
117  tmp<symmTensorField> tRes = New(tf);
118  inv(tRes.ref(), tf());
119  tf.clear();
120  return tRes;
121 }
122 
123 
124 template<>
126 (
127  const tensorField& tf
128 )
129 {
130  return symm(tf);
131 }
132 
133 template<>
135 (
136  const tmp<tensorField>& ttf
137 )
138 {
140  ttf.clear();
141  return ret;
142 }
143 
144 
145 template<>
147 (
148  const symmTensorField& stf
149 )
150 {
151  return stf;
152 }
153 
154 template<>
156 (
157  const tmp<symmTensorField>& tstf
158 )
159 {
160  return tstf;
161 }
162 
163 
164 // * * * * * * * * * * * * * * * global operators * * * * * * * * * * * * * //
165 
167 
170 
171 
172 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
173 
174 } // End namespace Foam
175 
176 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
177 
178 #include "undefFieldFunctionsM.H"
179 
180 // ************************************************************************* //
#define BINARY_OPERATOR(ReturnType, Type1, Type2, Op, OpName, OpFunc)
Field< symmTensor > symmTensorField
Specialisation of Field<T> for symmTensor.
void clear() const
If object pointer points to valid object:
Definition: tmpI.H:230
#define UNARY_FUNCTION(ReturnType, Type1, Func, Dfunc)
T & ref() const
Return non-const reference or generate a fatal error.
Definition: tmpI.H:174
dimensionedSymmTensor sqr(const dimensionedVector &dv)
dimensionedSphericalTensor inv(const dimensionedSphericalTensor &dt)
const Cmpt & z() const
Definition: VectorI.H:87
Vector< scalar > vector
A scalar version of the templated Vector.
Definition: vector.H:49
dimensionedScalar det(const dimensionedSphericalTensor &dt)
#define TFOR_ALL_F_OP_FUNC_F(typeF1, f1, OP, FUNC, typeF2, f2)
Definition: FieldM.H:112
SphericalTensor< Cmpt > sph(const DiagTensor< Cmpt > &dt)
Return the spherical part of a diagonal tensor.
Definition: DiagTensorI.H:288
const tensorField & tf
dimensionedSymmTensor twoSymm(const dimensionedSymmTensor &dt)
const Cmpt & y() const
Definition: VectorI.H:81
void dot(FieldField< Field1, typename innerProduct< Type1, Type2 >::type > &f, const FieldField< Field1, Type1 > &f1, const FieldField< Field2, Type2 > &f2)
tmp< DimensionedField< TypeR, GeoMesh > > New(const tmp< DimensionedField< TypeR, GeoMesh >> &tdf1, const word &name, const dimensionSet &dimensions)
Spatial transformation functions for primitive fields.
dimensionedSymmTensor dev(const dimensionedSymmTensor &dt)
Pre-declare SubField and related Field type.
Definition: Field.H:57
SymmTensor< scalar > symmTensor
SymmTensor of scalars.
Definition: symmTensor.H:48
dimensionedScalar tr(const dimensionedSphericalTensor &dt)
A 1D vector of objects of type <T>, where the size of the vector is known and can be used for subscri...
Definition: HashTable.H:61
const Cmpt & x() const
Definition: VectorI.H:75
void hdual(pointPatchField< vector > &, const pointPatchField< tensor > &)
Templated 3D Vector derived from VectorSpace adding construction from 3 components, element access using x(), y() and z() member functions and the inner-product (dot-product) and cross product operators.
Definition: Vector.H:57
dimensioned< scalar > magSqr(const dimensioned< Type > &)
dimensionedSymmTensor innerSqr(const dimensionedSymmTensor &dt)
dimensionedSymmTensor dev2(const dimensionedSymmTensor &dt)
dimensionedSymmTensor symm(const dimensionedSymmTensor &dt)
#define UNARY_OPERATOR(ReturnType, Type1, Op, OpFunc, Dfunc)
SphericalTensor< scalar > sphericalTensor
SphericalTensor of scalars.
A class for managing temporary objects.
Definition: PtrList.H:53
label size() const
Return the number of elements in the UList.
Definition: UListI.H:299
#define BINARY_TYPE_OPERATOR(ReturnType, Type1, Type2, Op, OpName, OpFunc)
tmp< Field< symmTensor > > transformFieldMask< symmTensor >(const tensorField &tf)
dimensionedSymmTensor cof(const dimensionedSymmTensor &dt)
Namespace for OpenFOAM.